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Chapter 1 
Introduction 


A uniform geometrical optics (UGO) and an extended uniform geometrical theory 
of diffraction (EUTD) are developed for evaluating high frequency electromagnetic 
(EM) fields within transition regions associated with smooth caustics of reflected 
rays and composite shadow boundaries. Although such transition regions usually 
occupy a small portion of the coordinate space surrounding an opaque scattering 
object, they contain some very complicated and yet most interesting wave phenomena. 
Composite shadow boundaries occur when a smooth caustic formed by an envelope of 
geometrical optics (GO) [1] reflected rays is terminated due to an edge in the reflecting 
boundary, or when the smooth caustic is in the immediate vicinity of the reflection 
shadow boundary (RSB). Figures 1.1 and 1.2 depict radiation/scattering problems 
that represent typical situations of interest for this analysis. In the concave surface 
example of Figure 1.1 the composite shadow boundary is formed by the termination 
of the smooth caustic occurring at T. This termination also generates a complex ray 
shadow boundary (CRSB) on the dark side of the smooth caustic. In the case of 
the concave- convex boundary of Figure 1.2 the smooth caustic is not terminated and 
therefore a CRSB does not exist; however, a transition region is formed when the 
edge at Q a is near the point of inflection, Q p , and the RSB is in close proximity of 
the smooth caustic. 

The UGO is a uniform version of the classic GO. It retains the simple ray optical 
expressions of classic GO and employs a new set of uniform high frequency reflection 
coefficients that provide for a finite reflected field near a smooth caustic on the lit 
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Reflection Shadow Boundary (RSB) 


Figure 1.1: Scattering and diffraction from a concave 


boundary containing an edge. 
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Figure 1.2: Scattering and diffraction from a concave- convex boundary containing an 
edge. 
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side. The UGO also encompasses a uniform version of the complex GO (CGO) ray 
field of Ikuno and Felsen [2, 3] that exists on the dark side of the smooth caustic. The 
EUTD is meant as an extension of the classic uniform geometrical theory of diffraction 
(UTD) of Kouyoumjian and Pathak [4] which is regarded as a uniform version of 
Keller’s geometrical theory of diffraction (GTD) [5]. The UTD provides for a finite 
field associated with the diffracted rays and also compensates the discontinuities of 
the classic GO fields across the optical boundaries. However, in the immediate vicinity 
of caustics the incident and/or reflected fields do not exhibit the simple ray optical 
behavior predicted by classic GO and therefore the UTD fails to provide the necessary 
corrections near a composite shadow boundary formed by the confluence of a smooth 
caustic and a ray shadow boundary. The EUTD developed in this report accounts for 
the non-ray optical behavior of the UGO reflected field near caustics by using a more 
sophisticated transition function in the expression for the edge diffraction coefficient. 
It also uniformly recovers the classic UTD behavior of the diffracted field outside the 
composite shadow boundary transition region. 

The approach employed in this report to construct the UGO/EUTD solution for 
the high frequency EM fields near smooth caustics and composite shadow boundaries 
is based on a spatial domain radiation integral representation for the fields scattered 
by a perfectly conducting smoothly indented surface containing an edge. The radi- 
ation integral is then evaluated using uniform asymptotic procedures to recover the 
uniform reflected, complex ray, and edge diffracted fields. The asymptotic estimate 
for the surface current in the radiation integral representation is found using the phys- 
ical optics (PO) [6] approximation and therefore the resulting solution only includes 
first order scattering mechanisms that remain valid away from regions where grazing 
fields exist. In addition, the resulting expressions for the edge diffracted field are 
heuristically modified to satisfy reciprocity and the local boundary condition on the 
scattering surface. Two methods are employed for the asymptotic reduction of the PO 
radiation integral. The first method employs the Chester et al. expansion [7] which 
results in a rigorous asymptotic reduction of the PO radiation integral. Although the 
resulting solution is not in a UGO/EUTD format, that is the expression for the fields 
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is not in terms of the conventional ray optical GO/GTD description multiplied by an 
appropriate transition function, it provides for a systematic approach for obtaining 
the higher order terms in the uniform asymptotic expansion for the fields in cases 
where slope scattering effects are present. The second method for the asymptotic 
reduction of the PO radiation integral employs the method of steepest descent [8]. 
The resulting expressions in this case are indeed in a UGO/EUTD format and are in 
agreement with the results of the first method to the first order term. Both solutions 
reduce to the conventional GO/UTD and CGO expressions outside the transition 
regions associated with the smooth caustic and composite shadow boundary. 

The subject of high frequency field behavior near a terminated caustic of a com- 
posite shadow boundary was investigated by Levey and Felsen [9] using a boundary 
layer expansion. However, they did not provide a useful engineering solution since 
they only concentrated on the mathematical aspects of the problem and specifically 
on the asymptotic reduction of the relevant diffraction integral and the asymptotic 
properties of the incomplete Airy functions that serve as canonical functions in the 
uniform asymptotic representation for the fields. On the other hand, the subject of 
evaluating fields in regions of smooth caustics but away from caustic terminations has 
been investigated extensively in the past. The solution obtained by Rahnavard and 
Rusch [10] is strictly valid in the caustic transition region, where Albertsen et al. [11] 
heuristically modified their solution near caustics of the diffracted rays from a curved 
edge to reduce to the non-uniform result on the deep lit side of the caustic. Other 
solutions which use Maslov’s method [12, 13] and uniformly recover the non-uniform 
fields outside the caustic transition regions are not very useful for engineering pur- 
poses since the relevant parameters involved in the solutions depend on the details of 
the caustic geometry. A uniform asymptotic high frequency solution which remains 
valid across smooth caustics of reflected rays, uniformly recovers the non-uniform ex- 
pressions outside the caustic transition region, and is useful for engineering purposes 
was recently developed by Pathak and Liang [14]. Their solution, however, is not 
cast in a UGO format, that is the scattered field contributions do not retain the GO 
ray optical expressions multiplied by an appropriate transition function, and also the 
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effects of the caustic termination of a composite shadow boundary are not accounted 
for. 

The ultimate goal of the work presented in this report is to develop a useful 
and efficient engineering solution for treating high frequency radiation and scattering 
problems involving higher order geometric primitives such as polynomial or spline 
defined surfaces that are commonly used in computer aided design (CAD) systems 
for the modeling of complex structures such as aircraft, ships, reflectors, etc. Such 
higher order surfaces often exhibit rapid curvature variation and discontinuities in sur- 
face derivatives which in turn create caustics and phase catastrophes that cannot be 
handled by the classic GO/UTD techniques. Numerical techniques for treating these 
problems such as the method of moments (MM) or PO with numerical integration still 
require vast amounts of computer resources for even intermediate size objects. Uni- 
form asymptotic closed form solutions on the other hand are inherently very efficient 
provided the canonical functions involved are efficiently computed. They also provide 
significant physical insight into the various scattering and radiation processes, a fact 
that proves useful in radar cross section (RCS) as well as antenna pattern analysis 
and design. 

The outline of this report is as follows: In Chapter 2, the near-zone problem for 
the scattered EM fields near a two dimensional smooth caustic and composite shadow 
boundary is formulated in terms of a spatial domain radiation integral which is then 
reduced using two different uniform asymptotic procedures. The first method yields 
a rigorous uniform asymptotic solution that provides for a systematic approach in ob- 
taining the higher order terms in the uniform asymptotic expansion for the fields. The 
second method yields the desired UGO/EUTD expressions for the reflected, complex 
ray, and edge diffracted fields which are in agreement to the first order term with 
the results of the first method. The same procedure is repeated in Chapter 3 for 
the near-zone problem for the scattered EM fields near a three dimensional smooth 
caustic and composite shadow boundary. In Chapter 4, the f&r-zone RCS problem 
for the plane wave scattering from two and three dimensional polynomial defined 
surfaces is formulated using the procedures developed in Chapters 2 and 3. In the 
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far-zone problem the two methods introduced in Chapters 2 and 3 produce the same 
results for the reflection and zero-curvature diffraction mechanisms with the associ- 
ated field expressions cast in the UGO/EUTD format. However, they result in slightly 
different expressions for the EUTD edge diffraction coefficient. Although the zero- 
curvature diffraction contribution in the far-zone problem is equivalent to the complex 
ray mechanism in the near-zone problem, it is easier to use since it eliminates the 
need for a complex extension of the scattering surface. In Chapter 5, some indicative 
numerical results for the scattering from cubic and fourth order polynomial strips 
are presented. The UGO/EUTD results are validated by comparison with results ob- 
tained using an independent method of moments (MM) solution. The UGO/EUTD 
results are also compared with classic GO/UTD results to demonstrate the need for 
the new solution. Chapter 6 summarizes the accomplishments of this report, offers 
some concluding remarks, and briefly discusses topics for future research. All the 
new UGO/EUTD reflection and diffraction coefficients are also tabulated in Chap- 
ter 6. Several appendices are also included in this report. Appendices A-E provide 
important mathematical details on the asymptotic properties and efficient numeri- 
cal evaluation of the canonical functions involved in the UGO/EUTD reflection and 
diffraction coefficients. Appendices F and G provide some background derivations for 
the uniform asymptotic reduction of the spatial domain PO radiation integral. 

All fields are time harmonic with an assumed e JUJt time dependance that is sup- 
pressed throughout this report. Also the wavenumber, A, of free space is assumed to 
have a small negative imaginary part to ensure that the radiation condition at infinity 
is satisfied. 
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Chapter 2 

Near- Zone Uniform Asymptotic 
Analysis and UGO/EUTD 
Formulation for the Scattered 
Fields near a Two Dimensional 
Smooth Caustic and Composite 
Shadow Boundary 


In this chapter, the behavior of near-zone high frequency scattered EM fields near a 
two dimensional smooth caustic and composite shadow boundary is investigated and 
a UGO/EUTD solution is developed for the description of the scattering mechanisms. 
A composite shadow boundary occurs when a smooth caustic formed by an envelope of 
GO reflected rays is terminated due to an edge in the reflecting boundary, or when the 
smooth caustic is in the immediate vicinity of the RSB as shown in Figures 1.1 and 1.2. 
The caustic termination also generates a complex ray shadow boundary (CRSB) on 
the dark side of the smooth caustic. When the RSB is not in the immediate vicinity of 
the smooth caustic, the conventional UTD edge diffraction solution [4] which involves 
the Fresnel integral as a canonical function can be used effectively together with the 
solution in [14] to describe the fields away from the composite shadow boundary. 
However, near the composite shadow boundary where a confluence of both reflected 
and caustic type shadowing occurs, the Fresnel integral is not adequate to describe 
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the diffracted field behavior and it must be appropriately replaced by the incomplete 
Airy functions [9] via a uniform asymptotic procedure. 

The present analysis is based on a spatial domain radiation integral representa- 
tion for the fields scattered by an analytically described smoothly indented perfectly 
conducting boundary containing an edge. The asymptotic estimate for the surface 
current induced by the incident field from a source whose vector pattern function 
can also be described analytically is found using the PO approximation. The PO 
radiation integral is subsequently reduced using two uniform asymptotic methods to 
recover the uniform expressions for the scattering mechanisms. The first method em- 
ploys the Chester et al. expansion [7] and results in a rigorous uniform asymptotic 
reduction of the PO radiation integral. Although the resulting solution is not in the 
desired UGO /EUTD format, that is the expression for the fields is not in terms of the 
conventional ray optical GO/GTD description multiplied by an appropriate transition 
function, it provides for a systematic approach for obtaining the higher order terms 
in the uniform asymptotic expansion for the fields in cases where slope scattering 
effects are present. The second method for the uniform asymptotic reduction of the 
PO radiation integral employs the method of steepest descent [8] and the resulting 
expressions in this case are indeed in a UGO/EUTD format. 

The PO integral for finite concave or concave-convex shaped boundaries is char- 
acterized by two stationary phase points that are arbitrarily close to one another or 
to the integration endpoint associated with the edge. The stationary points are real 
and correspond to the two points of reflection when the observation point is on the 
lit side of the caustic. For observation points near the smooth caustic itself, the two 
stationary points coalesce to form a second order stationary point that moves near the 
integration endpoint for observation points near the caustic termination or composite 
shadow boundary. For observation points on the dark side of the smooth caustic, the 
two stationary points become a pair of complex conjugate stationary phase points 
with only one of them satisfying the radiation condition and eventually contributing 
to the fields on the dark side. Away from the caustic transition region the contri- 
bution from the complex stationary point was interpreted by Ikuno and Felsen [2, 3] 
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as a complex GO ray field reflected from the complex extension of the boundary. 
This complex GO ray field obeys a generalized form of Fermat’s principle with all 
the parameters involved computed in complex space through analytic continuation 
and in precisely the same manner as for real GO rays as shown in [14]. In the case 
of a caustic termination as in Figure 1.1, the complex ray field only exists in the 
illuminated side of the CRSB, and therefore the edge diffracted field has to provide 
continuity in the total field across this boundary as well as the optical boundaries 
associated with the real GO rays on the lit side of the caustic. 

Although higher order mechanisms such as edge excited creeping waves and whis- 
pering gallery modes may also exist and contribute significantly to the scattered field 
from the geometries considered, they will not be addressed in the analysis that follows. 
The problem formulation and the derivation of a uniform asymptotic solution using 
the Chester expansion for the asymptotic reduction of the PO radiation integral are 
presented in Section 2.1. In Section 2.2 the method of steepest descent is employed for 
the asymptotic reduction of the PO radiation integral and the UGO/EUTD expres- 
sions for the reflection, complex ray, and edge diffraction mechanisms are formulated. 

2.1 Derivation of a uniform asymptotic solution 
using the Chester expansion 

In this section, the spatial domain PO radiation integral is formulated and then 
asymptotically evaluated using a rigorous uniform asymptotic procedure that employs 
the Chester expansion. The canonical geometry for the analysis that follows is a semi- 
infinite two dimensional concave or concave- convex boundary as shown in Figures 2.1a 
and 2.1b, respectively. The scattered electric field at any point P exterior to the 
surface boundary C can be expressed in terms of the usual radiation integral over the 
electric current J induced on C by the source at P' as follows: 

2?'(P) » _» ^ [8 x R X J{Q')] H™{kR)dl' (2.1) 

where 

R = QP, J2= |J?|, R = R/R, 
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Figure 2.1: Canonical geometries for the uniform asymptotic analysis: (a) semi- 
infinite concave boundary, and (b) semi-infinite concave-convex boundary. 
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J(Q') 

H^(kR) 

Z 0 

k 

dt 


value of J at any point Q' on the surface, 

cylindrical Hankel function of the second kind of order zero, 

impedance of free space, 

wavenumber of free space, and 

line integration element. 


Also, the observation point P cannot be too close to the scattering boundary for 
Equation (2.1) to be valid. When the scattering boundary is electrically large and 
well illuminated, one can use the PO approximation to (2.1) in which the induced 
current is assumed to be given by GO as follows: 


J(Q') 


< 


2to' x H'(Q ') , on the lit portion of the boundary 
0 , on the shadowed portion of the boundary 


( 2 . 2 ) 


in which n' is the outward unit normal to the surface at Q' y and H'(Q ') is the direct 
incident magnetic field from the source whose phase center is at P* and may expressed 
ray optically as follows: 


H‘m ~ I S' A = 0, 

V s’ 

= Z;'S‘xE‘(Q'), 

E‘(Q') = A(Q') ^rr, 

Vs’ 


(2.3) 

(2.4) 

(2.5) 


where A is the vector source pattern, s' — P'Q', s' = |J* ‘|, and s' = s' I s'. Using the 
assumed current in (2.2) and the large argument form of Hankel function given by 

H?\kR) ~ e~ ihR , for kR > 1, (2.6) 

—* 

E‘(P) is given by the following stationary phase integral: 


where 


E‘{P) w j F{i-)e- ik ^W 

la 


F{1) = 

m = 


I Jk R(t') x R(t') x [n(<') x ?•(<') x !(<■)] 
2 tt 


y/sWW) 


, and 


*•(/') + R(l ') . 


(2.7) 
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(2.8) 

(2.9) 



The phase function contains two real or complex stationary phase points l ia satisfying 
the condition 4>{l, a ) = 0, and are arbitrarily close to one another or to the integration 
endpoint l a . Before proceeding with the asymptotic reduction of the integral in (2.7), 
we choose the reference point on the scattering surface to be Q 0 , which is a first order 
stationary point of the derivative of the phase function, that is <j >(£ 0 ) = 0 and must 
necessarily lie between the two stationary phase points on the concave part of C. For 
the moment, we assume that this point is real without losing the generality of the 
analysis. Then, the quantities R(t') and s'(i') are given by 

R(t) = \Ro- W| > and ( 2 - 10 ) 

= |aj+f(0|, ( 2 - n ) 


where H = Q 0 P, « j = PQo, and r(t) = Q 0 Q' is the natural representation of the 


scattering boundary. 

The first step in the asymptotic reduction of the stationary phase integral in 
Equation (2.7) is to expand the phase function in a Taylor series around £ 0 , that is 


^ 4>{£o) + {£' ~ toWo) + 


(l'-*o) 2 

2 ! 


m + 


V - *»? 

3! 


3 _ 


Hi 0 ). 


( 2 . 12 ) 


The derivatives of the phase function can be expressed in terms of the geometric 
parameters of the surface and using (F.9), (F.14), and (F.19) they are given by 


4>{lo) 

m 

?(«.) 


— (^ o — Ro) • to j 

= n 0 • (Ro ~ 3 0 )K g (Q 0 ) + : 1- „ - u > 




Ro 


= n 0 • (Ro — a * ) k' 9 (Q 0 ) + t 0 -{Ro- r 0 ) k)(Q 0 ) 


+ 


+ 


3(?; •£„)(?' - to) 




3(f2o * w 0 )(J?o ■ to) 

Ro 


, r \ \ a o '-0 

*»(<? 0) — “ 

*.m + 


n, 

8 

Ti ■ n„ 


n, 


= m (Qo) , 


(2.13) 

(2.14) 


(2.15) 


where t 0 and n 0 are the unit tangent and normal vectors to C at Q 0 1 respectively, 
and /c e ((2o) = ]■! pg{Qo) is the normal surface curvature at Q 0 . From Equation (2.14) 
one can observe that the point l 0 is real for a concave or concave- convex boundary, 
with k„(Q 0 ) < 0. On the other hand, for a strictly convex boundary, the point £ 0 , 
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if it exists, must be complex. In the special case of plane wave incidence with the 
observer in the Fraunhofer region of the boundary, that is s ' 0 , Ro -> oo, the point Q 0 
is a zero-curvature (inflection) point that satisfies k 9 (Q 0 ) = 0. For this case, l 0 may 
be real for strictly convex boundaries as well. 

Next, we make the following linear transformation: 



(2.16) 


where m(Q 0 ) is given in (2.15), and Equation (2.9) becomes 



<f>(l ) = t(s) = a + (3s + «s 3 /3 , 

(2.17) 

— 

where 



“ = m) = \W>) + <Kl>)\ , or 

(2.18) 


= «<...) ±|(-r, or 

(2.19) 


= «<.)- 0s. - s’J3 , 

P = m L<w] ’ ” 

(2.20) 

(2.21) 


= - W*) - ^(A)]J , or 

(2.22) 


= *{f !*&)-*(*,,,)]} ‘ , and 

(2.23) 

— 


(2.24) 

— 

= ±[^)-M)rf-4rl‘; j + “ Co3C 

L m (^o)J [ - if Qo G C. 

(2.25) 

— 

The proper branches for /? and s a depend on the sign of m(Q 0 ). Thus, 
Equation (2.7) transforms to the following expression: 

using (2.17) 

— 

"oo” 

E‘(P) w e~ ika 1 G(s)e- jk(0 ‘ + ’ i/3) ds 

•a 

(2.26) 


where 


— 

Co 

II 

8-1% 

(2.27) 

— 
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dt 

ds 


m 


ml 


, and 


“oo” — oo exp 


(2.28) 

(2.29) 


(£)] • 

The proper branch for ^ also depends on the sign of m(Q 0 ). 

Next we employ the Chester et al. expansion [7] for the amplitude function in 
Equation (2.26), that is 

GW = f; [«.(•' +/>r] 

m — 0 

= a 0 + sb 0 + (s 2 + fi)g («) , (2.30) 


where 


a. = i [<?(.,) + <?(.,)] . 

= i [<?(«,)- <?(«,)] . 


K 

m = t, N* 1 + /’)■" + + «”"] - 


>1,2 
G{s 1>2 ) 


dt 

ds 


= ±(-/3) 1/2 ; rKO = 0, 

dt I 




ds 


* 1,2 


= (-/») 


1/4 


* 1,2 


±2 


V(a, 2 ) 


(2.31) 

(2.32) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 


and thus Equation (2.26) becomes 


E‘(P) — e 


-jka 


“oo" OO 

d 0 I e- jk{e,+ ’ S/3) ds + b 0 j se- ik{p ‘ +,S/3) ds 


“oo" 

+ / (s'+P)g(s)e-w+‘ S ' 3 'ds 


(2.37) 


The third integral on the right hand side of Equation (2.37) is regular over the entire 
path of integration, thus employing integration by parts for this integral we get 


E‘(P) ~ e 


-jka 


'oo oo 

do I e - 2 W S/ » ds + b 0 J se~ Mti,+t ' l3) ds 

• a 

+ 0(k-’). 


JK 


(2.38) 
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Now, from Equation (2.30) we have that 


-* / \ ®(^«) <*0 Sab 0 

9M= K + fi 

and omitting the 0(k~ 7 ) terms, Equation (2.38) becomes 


(2.39) 


E‘(P) 


~ e 




+ € 


~~j kr (*a) 


00 "oo 

a 0 I e -w+‘ s /») da + b 0 J s ds 

1 

G[sg) a 0 + s o b 0 


[jfcr'(s a ) jk(/3 + 3|) 

Next, we let s = k~ l/3 t, 0 = k~ 2/3 7, and s a = *r 1/3 £ 0 in (2.40), thus 


(2.40) 


E‘{P) 


D -jka 


y/k 


"00 

a° J 


e-^ t+t3 ' 3 'dt + b 0 j te~^ t+,3/i) dt 

Ca C a 

G(*a) j(flQ + Ca^o) 

i*7-'(s a ) v £( 7 + £) 


(2.41) 


Finally, depending on the sign of m{Q 0 ) the scattered electric field E’(P) is given by 


E'(P) 


\/k 


d 0 Ai'( 7 ,C a )+i6o J^Ai*(7,C.) 


F(l a ) 3 (no + Cbo) 


E'(P) 


o~J ka 


Vic [ 


[mi a ) v^(7+oJ 

3 

“oAi(7, Ca) - jbo ^Ai(7, Co) 


, if m(Q 0 ) > 0 (2.42) 


F(t a ) 3 (a 0 + Cab 0 ) 


, if m(Q 0 ) < 0, (2.43) 


[mi.) vS(7+a 

where (*) denotes the complex conjugate operator, 7 = k 2/3 0, ( a = k 1/3 s a , and a, a 0 , 
and b 0 are given by the following expressions: 


a 0 

b 0 


a = W a )-6k~'( 7C. + C/3), or 
= *(<,.,) ±tt- |(— r)*" , 

S = sgn[m((J 0 )) , 

( — l) ±tJt 




\ 


IM)' 


± F(Q 


i 


lto)l 


(2.44) 

(2.45) 

(2.46) 

(2.47) 
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(2.48) 


The function Ai(7,£) is the incomplete Airy integral defined by 

oo 

Ai( 7 ,0 = / e> ( -" +zi ' 3) dz. 

( 

The second term in both Equations (2.42) and (2.43) is an edge diffraction con- 
tribution where the first term in (2.42) and (2.43) appears to be a contribution from 
the two real or complex stationary points. If left in this form, the solution would 
clearly be in a UAT format [15], since the necessary corrections to the fields near 
the optical boundaries are performed on the reflected field as opposed to the edge 
diffracted field in a UTD type solution [4]. In fact, the first term in Equations (2.42) 
and (2.43) besides containing reflection and complex ray mechanisms, it also contains 
edge diffraction contributions implicitly embedded in the incomplete Airy integral. 
Therefore, in order to obtain a UTD type solution it is necessary to separate the edge 
diffraction contributions in (2.42) and (2.43) by using the properties and complex 
plane topology of the incomplete Airy integral found in Appendix C. In this form, 
the solution would be applicable to more general problems and also provide important 
physical insight into the individual scattering processes. 


2.1.1 Scattered field solution on the lit side of the caustic 


On the lit side of the caustic, the two stationary phase points i 1<2 are real and cor- 
respond to real reflection points Q l 2 on C. Therefore, the parameter 7 in Equa- 
tions (2.42) and (2.43) must be negative. Using (C.9)-(C.ll), the scattered field 
solution on the lit side of the smooth caustic is given by 
1. m((g 0 ) > 0 : 


EW ~ ivf e_i ‘“ { [‘WM + + A W .'W] } 


. e -J»«*a) f E(lg) On 

1 jkW.) Vk 1 




, • fee 

3 Vk [ 


37 M7 ’ W +c 


}■ 


(2.49) 


if (.<-(-7)-", 
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+ e 


-jk<p(ta) 


f F(Q 3. 
\jki(l.) A 


e «.*d/> )g ; ( ., X) + 


+ 3 


. b 0 


y/k 


J(l(a + <J/»)A e 'fT, t \ _L . A a ... - l 

a 1 g - (1,c * ) + 7 + cJj ’ 


if -(- 7) 1 ' 1 < C. < (- 7 )-" , 


£'(/>!> ~ e -M‘.) | f i <J A + °° 


^ .K_ 

3 A 


if C. > (- 7 ) 


jk<j>(l a ) y/k L 

e i(-rt.+d/»iA-*f~ t \ I ^ 

87 g,l7 ' W 7 + «j 


e ««.*<S/.> g;( 7 i(a) + _l_ 


}' 


- 1/2 


(2.50) 


(2.51) 


2. m(Qp) < 0 : 

EW ~ -isj\ { [«oW,( 7 ) - j 6 „w;( 7 )] - fcw.ft) - i*„w;( 7 )] } 

| g— 


- 3 


. 60 


Vk 

if C<-(- 7)- 1/2 


J EUl il _L A- g-i( 7 <a + d/S) e f\ i 

a [ 84 J 7 +ej 


}■ 


E-(P l ) ~ i v /|e- i *”[S.W 1 (7)-,iw;(7)] 


+ e 




f F(Q . aq 

1^4) V* 


g XKa-K*/») gj ( 7 ^ a ) _ _ 


. bo 


g->(7(a + d/S)j?_- M + 

5 7 gs(7 ’ U + 7 + Cj 


Vfc 1 

if ~(- 7) 1/2 < C. < (~ 7)- 1/2 


}• 


E-{P l ) 


~ e 


- 3 


-}k*(ta) 


/ ~^ a ) 4. -^2- fg-iCKa+Cj/*)- ^ / \ 2 

\ + vlfe l gl(7,Caj 7 + 0 


d 


y/k [ 

if C>(- 7)-‘ /2 


p-i(nrc« + cj/®) ^ _L 

»r g,l7 ’ w + 7 + C 


_6_n 

-r + cj J ’ 


(2.52) 


(2.53) 


(2.54) 


where W 12 ( 7 ) are the Fock-type Airy functions [16], and g,(7,C<.)? * = 1,2,3 are 
the incomplete Airy functions [9]. From (2.8), (2.9), (F.9), and (F.28) we have the 
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AW - i/s 


following relations on the lit side of the caustic: 

j jk S 1,2 X 3 1,2 X [ n i»2 X 5 1,2 X -^(Ql.a)] 

\j S 1,2 ^1,2 

^(^ 1,0 = S l ,2 + 5 i ,2 1 

*w = (^+^). 


and 


0(/o) = («J - &) *<o = -(cosy’o + C°sy7 0 ), 


(2.55) 

(2.56) 

(2.57) 

(2.58) 


where 0\ 7 is the usual angle of incidence at Q l>2 , y> 0 are defined in Figure 2.1, 

and p\ 2 is the radius of curvature of the reflected wavefront from Q 12 given in (F.29). 
Also, the vector expression in the numerator of (2.55) may be written as follows: 


^1,2 X ^ 1,2 


ft'l ,2 


X 3 


X 


A(Q 1i2 )] = cos 0[ 2 A(Q^)- R (<?i, 2 ), 


where R (Q 12 ) is the dyadic GO reflection coefficient and is given by 


R (<?i, 2 ) = flj.s »i,a — I — R* e±e± + Rh e J|(^ I<2 )e jj(Qi, 2 ) > 
h 1 ; for a perfectly conducting boundary, 


(2.59) 


(2.60) 

(2.61) 


/ is the unit dyad, and the vectors e x and e | ( ’ r are defined in Figure F.2. Thus, 
using (2.21)— (2.23), (2.45), (2.47), and (2.55)-(2.59), the parameters a, 7, a 0 , and b 0 
for the lit side of the caustic are given by 


a 


= 3\ t + s[ 7 ±6k , 


I 2 I* 

7 = -fc 2/s |cos¥>; + cosyjol — 77^-d , or 

= - |* !(«! + *;) -(«! + *:)!]* , or 


fl-0 

60 


I* 


(•So + Ro) — («!,, + ^,2 

MQYi RjM I ir 

\ pi + *l 


Ik (-7) ±1/4 
7T 




I el 

p; + 


(2.62) 

(2.63) 

(2.64) 

(2.65) 


(2.66) 
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Figure 2.2: Scattered field regions associated with a concave boundary containing an 
edge. 

Before proceeding with the scattered field solution on the dark side of the caustic, 
a few observations on Equations (2.49)-(2.54) are in order. First, Equations (2.49) 
and (2.52) apply on the illuminated side of the reflection shadow boundary (RSB) 
designated as region I in Figures 2.2 and 2.3, where two reflected ray contributions 
exist and are represented by the first term in (2.49) and (2.52) involving the Fock-type 
Airy functions and their derivatives. The second term in (2.49) and (2.52) represents 
the diffracted field contribution from the edge. Equations (2.50) and (2.53) apply 
on the shadow side of the RSB designated as region II in Figures 2.2 and 2.3, where 
only a single reflected ray contribution exists and is represented by the first term 
in (2.50) and (2.53). The second term in (2.50) and (2.53) represents the diffracted 
field contribution from the edge. Equations (2.51) and (2.54) apply in regions III and 
VI in Figures 2.2 and 2.3 where only edge diffracted rays exist. 
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Figure 2.3: Scattered field regions associated with a concave-convex boundary con- 
taining an edge. 
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2.1.2 Scattered field solution on the dark side of the caustic 


On the dark side of the caustic, the two stationary phase points £ 12 are complex 
and correspond to complex reflection points Q\ 2 which are interpreted to lie on the 
complex extension of C . Therefore, the parameter 7 in Equations (2.42) and (2.43) 
must be positive. Using (C. 12 ) and (C.13), the scattered field solution on the dark 
side of the smooth caustic is given by 


1. m(Q q) > 0 : 

E-( p s ) ~ ^ e-*" [s.Ai( 7 ) + j6 0 Ai'(f)] 

/ jfyk) i °° 

\jk4>(l a ) y/k [ 


+ e' 


e i(,c a+ c 2 / 3)g . (7jCa) + _J_ 


1 -A_ 

+ 3 Vk 


0 >(TCa + C /* \ , 


if Ca< 0, 

E\p s ) ~ £->**(<<■> J 

VS [ 


)• 


^(^Ca + d/ 3 )^* 


g*(7,Ca) + 


7 + C 


+ j 


b 0 


2 . 7n(Q 0 ) < 0 : 


\/k 

if Ca>0 

27 r 


[ »7 8,l7,c ‘ j 7+dr 


E'(Ps) ~ Tje-»"[«.Ai( 7 )-,iAi'( 7 )] 


+ e 


-ik<ne a ) / EiEl _j_ °° 




e -i(,C a + d/3) g2(7tC) _ 


7 + CoJ 


~ J 


. b 0 


Vk[ 

if Ca < 0 , 


+ (~ (■ \ I l 

a 7 g s (7,C») + 7 + CJ j | , 




~ e 


w<„) + v*r g ' l7 ’ c ■ , 7+c 


. K 
3 Vk 


« ( ,,(.) + _k_ 

^7 7 + Ca 


)' 


if Co > 0 > 


(2.67) 


(2.68) 


(2.69) 


(2.70) 
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where Ai(7) is the ordinary or Miller-type Airy function [17]. In the dark side of the 
caustic we have the following relations: 


, fjk cos e^AiQl,) ■ R‘(Q[,) 

* 1*1,2; — Vo^ Tz — 7 T 7 

v 27r yjs s[% 

1, 2 ) = , and 
fa*) = cos 2 (-^ + -M , 

\P 1,2 *1 ,2 / 

where R C {Q\ 2 ) is the dyadic complex GO reflection coefficient and is given by 


(2.71) 

(2.72) 

(2.73) 


r c (Q'J = k,k,~ i = ( 2 - 74 ) 


h = rpl ; for a perfectly conducting boundary. 


(2.75) 


Notice that all the quantities in (2.71)— (2.75) with the superscript c are associated 
with the complex reflection points <2, 2 and are analytic continuations along the com- 
plex extension of C [2] of the quantities associated with the real reflection points Q t 2 
in the lit side of the caustic. They are evaluated in precisely the same fashion as the 
quantities associated with the real reflection points. Thus, using (2.21) (2.23}, (2.45), 
(2.47), and (2.71)-(2.73), the parameters a, 7, a 0 , and bo for the dark side of the caus- 
tic are given by 


a = s[ c 2 + 3 r ^TjSk- l -Y /7 , 


7 = k i/3 |cos<^>' + cos v? 0 | 


m 


(<M 


IW + «;•) - (< + «r 


1 £ 


or 


or 






it 

— 

[2* « + 

R 0 )- 

(-Hi + '"*) 

5 

_ jy/k-y 1 '* 

o* 

1 

R c (Ql) 

/ />r 

a 0 — ■ 

2y/7r 

V 

G? \ 

p[ c + 

+ 

Aqd-r 

e {Ql ) 

1 PV 



i 

v 

5 

- y/k') 1/4 


•R c (Qt)\ 

/ 

^0 — 

2v^F 

! |\ 

1 p\ c + *r e 


(2.76) 

(2.77) 

(2.78) 

(2.79) 


(2.80) 
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(2.81) 


MQl)-R c {Ql) 


p; e + 


A few observations on the nature of the scattered field solution on the dark side of 
the caustic are as follows: First, Equations (2.67) and (2.69) apply on the illu min ated 
side of the CRSB designated as region IV of the dark side in Figures 2.2 and 2.3, 
where both complex ray and edge diffraction contributions exist. The complex ray 
contribution is represented by the first term in Equations (2.67) and (2.69) involving 
the ordinary Airy function and its derivative. Although the values associated with 
both complex stationary points are present in the expression for the complex ray 
contribution, as will be shown later in this section, in its asymptotic limit it contains 
only the reflected field from one of the complex stationary points. The second term 
in (2.67) and (2.69) involving the incomplete Airy function and its derivative repre- 
sents the edge diffraction contribution in region IV. Equations (2.68) and (2.70) apply 
in shadow side of the CRSB designated as region V of the dark side in Figure 2.2, 
where only edge diffracted field contribution exists. 

2.1.3 The reflected field and its limiting forms on the lit 
side of the caustic 

From Equations (2.49)-(2.54) and using the fact that W^) = W*( 7 ) the reflected 
field on the lit side of the caustic may be written as follows: 

e-'f 


E r (P L ) ~ [^ 7 y*w 1 (-7 / )-ifl/77 l/4 w;(- 7< )] 

- u 2 [At 7 ; /4 w,(-7<) - jSa7 l,t w;(- 7< )] } 

if m(Q 0 ) > 0 , 

E’(Pl) ~ - V {“I [W w;(- 7l ) + ]B, w;-(- 7 ,)] 


(2.82) 


- [AnY' W;(-7,) + W;-(-7,)] } 

if m(Q 0 ) < 0 , 

«.* = «(-c-±7n, 


(2.83) 


(2.84) 
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where u(i) is the Heaviside unit step function and from (2.62), (2.64), and (2.66) the 
parameters 7 1 , a 12 , A/, and Bt are given by 

7 , = [j*i(.i+o-M+*;)i] ,, ’> ( 2 - 85 ) 

( 2 . 86 ) 


<*1,2 = S ’l,2 + J l,2 i 1 

3 7 ' . 



f At ) A{Q 1 )-r{Q 1 ) 

/ Pl 

A{Q 2 )-r{Q,) 

J « 

\ B t j ^ 

V p\ + 


V P2 + *5 


(2.87) 


with 8 given in (2.46). The first term in (2.82) and (2.83) is non-zero in regions I and 
II of Figures 2.2 and 2.3 and represents the uniform reflected field contribution from 
Q,. The second term in (2.82) and (2.83) is non-zero only in region I of Figures 2.2 
and 2.3 and represents the uniform reflected field contribution from Q 2 . Both terms 
in (2.82) and (2.83) are zero in the remaining regions of Figures 2.2 and 2.3 where 
only edge diffracted and/or complex ray fields exist. Also, it should be noted from 
Equations (2.82) and (2.83) that the uniform contributions from each reflection point 
are not completely independent of each other since the values from both reflection 
points are explicitly contained in the parameters A/, and Bt- 

Next, let us investigate the behavior of the reflected field by examining two limiting 


cases. 


Limiting case 1 : Observation point on the deep lit side of the caustic. 

In this case, the parameter 7, > 1, and using the large negative argument forms 
of the Fock-type Airy functions and their derivatives given in (A.21)-(A.24) Equa- 
tions (2.82) and (2.83) become 

E r {P L ) ~ \ e-MW) + Bt) u (-C + 7/ /2 ) 

z 

+ i e -i*< •;+•;> [A t - Bt) u(~c - 7 \' 7 ) ; if m(Q 0 ) > 0, (2.88) 

z 

E t (Pl) ~ J (At + Bt) u (-C. + 7i /2 ) 

+ - e _jM * 2+ ' 5) (At - Bt) u(-Ca - 7* l/2 ) ; m (Q ° ) < ( 2 - 89 ) 

z 

From (2.5) and (2.87), we have that 

\e~ ik ‘ I* (At ± Bt) = J ^(Q,. a ) • R (<?i. 2 ) 


P\a 


P\a + 


(2.90) 
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and the reflected field on the deep lit side of the smooth caustic is given by the 
following expression: 


E’(P l ) ~ u, #(<?,)• R (Q,)J 

V Pi + a l 


+ u 7 E i (Q 3 )> R(Q,) . 


Pi 


Pi + *1 


e~ ik ’l 


(2.91) 


with the understanding that 

P\ 


pi + 

pi 

Pl +3 


7<0 for t/i(<? o )< 0 , 


->0 for m(Q o )< 0 , and 


1 ^ -3 

L P' 

if ^ 

/ P T + 3 r 

V V + * r 

P r + 3 r 


< 0 . 

Equation (2.91) shows that the reflected field solution on the deep lit side of the 
smooth caustic reduces to the classic GO expression as expected. 

Limiting case 2 : Observation point lies on the smooth caustic of reflected rays. 

In this case, 7* ► 0 in (2.82) and (2.83), where the spread factors in (2.87) 

become singular for observation points at the caustic, that is p [ 2 -> -s, 2 . From (2.28) 
and (2.36) the reflected ray spread factor near the caustic assumes the following form: 

cos 0' o (fc/2) 1/a 7 < ' 1/4 


l-K- 

Pl + 31 


V^MQo)| ,/S 


(2.92) 


and thus using (2.87) and (2.92) we have 


l\ h A t = 

7 « 


2 cos fl* (fc/2) 1/6 

o)| 1/s 

0. 


A(Qo)- R(Qo) 


and 


(2.93) 

(2.94) 


Also, the Fock-type Airy functions, W 12 , assume the following small argument form: 

w, , ( 0) = - : ’'‘^3 /3) (3'« ± jVa) , (2.95) 

and thus the uniform asymptotic expression for the reflected field close to the smooth 
caustic is as follows: 


E r (Pc) ~ «(-C-) 


cos (fe/6) 1/6 T(l/3) 


V* l m (Qo)|‘ /S 


E'(Q 0 )-R(Qo) 


r> k ‘ 0 

vW 


(2.96) 
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2.1.4 The complex ray field and its limiting form on the 
dark side of the caustic 

From Equations (2.67)-(2.70) the complex ray field on the dark side of the caustic is 
given by the following expression: 

E‘{P S ) ~ iVJe- 1 **!.* [A, 7 y‘Ai(i,)TB,7; l,, %)] 1l H)' (2-97) 

for m(Qo) >0 , 

where from (2.76), (2.78), (2.80), and (2.81) the parameters 7., <„ A„ and B. are 
given by 

(2.98) 

(2.99) 


7. = 


,2 

A. 

B, 


-k |(*“ + «;*)- 
L4 

(*; e + *r)l 

s\ C 7 + 2 -F 2 7a ? 

a(QD-r{QD 

/ pv 

V* 1 

\J p v + «r 

A(Ql)- R(Qi) 

1 p\ c 

V*? 

V P\ C + 


, 2/3 


( 2 . 100 ) 


The uniform complex ray held solution in (2.97) is non-zero only in region IV of 
Figures 2.2 and 2.3. Although the values associated with both complex rejection 
points are present in the expression for the complex ray held through the parameters 
A,, and B„ on the deep dark side of the caustic only the contribution from one of 
the complex reflection points should be present in the solution since the contribution 
from its complex conjugate counterpart violates the radiation condition and produces 
an unbounded result. On the deep dark side of the caustic the parameter 7 . » 1. 
and using the large argument form of the ordinary Airy function and its derivative 
given in (A.10) and (A.ll) the complex ray field in the deep dark side of the smooth 
caustic is given by the following expression: 

bips)~«(-c.) £'(<?:..) “«w> 0 ’ (2101) 

with the understanding that p" / (p™, + s”) < 0 on the dark side of the caustic. 
Equation (2.101) indeed shows that the complex ray field solution in the deep dark 
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side of the smooth caustic reduces to the non-uniform complex GO expression of 

Ikuno and Felsen [2, 3] obtained independently via the complex stationary phase 
method [18]. 


2.1.5 The edge diffracted field and its limiting forms 

The edge diffracted field may be written as follows: 

E d (P) = E' u (P) + E?(P), (2.102) 

where E nu (P ) and e c (P) are the non-uniform and curvature correction terms, respec- 

lively. From Equations (2.49)-(2.54) and (2.67)-(2.70) they are given by the following 
expressions: 


Ku(P) 

Eli P) 


IMA 




e i(^c a+ d/3 ) g ;( 7a)Ca ) + _i_^ 

7 a + Cj 


+ jB d Je j(7a<a+ ^ /3 )-^-ff*( 7 ] | 

i 5 7 0 g,l7a,Caj+ 7o +dr 

if m(Q 0 ) > 0 , 

E d c {P) ~ | e -i(T«Crr-Kj/ S ) g _( 7a _ j ^ j 

- jB d [ e -^°ca+c*/j)_^_ / / ) i 1 1 

[ ^ g,(7a,CJ+ 7o +dr 


% — 


if m(Q 0 ) < 0 , 

1 in regions III, V, and VI 

2 in region I and IV 

3 in region II 

where 7a , A d , and B d are given by 


J 

la = — Sk> (cos tp' 0 -(- cos <p 0 ) 


m{Q 0 ) 


i 


Ca — u 0 k^ |(cos ip' 0 -f- cos ^ 0 ) — (cos <p' a -f- cos 


m{Qo) 


u n = 


+1 if Qo B C 
-l if Q 0 eC, 


(2.103) 


(2.104) 


(2.105) 


(2.106) 

(2.107) 

(2.108) 
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U 1 

1 _ 

7a) ±1/4 

At 

. B ‘ J 

2y/n 

B t 


= — < 
2y^ 

f 3 A, 

1 B, 

' , if 7a 


V * ) 


, if 7a < 0 

(2.109) 

>0, 

(2.110) 


where A t , B t and A„ B, are given in (2.87) and (2.100), respectively. The non- 
uniform term in the diffracted field expression of Equation (2.103) can be recognized 
as the non-uniform PO solution for the edge diffracted field from a half- plane. The 
curvature correction term given by Equations (2.104)-(2.105) cancels the singular 
non-uniform term near the real ray optical boundaries, namely the RSB and ISB. 
It also compensates the discontinuities of the incident, reflected, and complex ray 
field components at the ISB, RSB, and CRSB, respectively, thus ensuring continu- 
ity in the total field. As will be shown in more detail later in this subsection, the 
curvature correction term reduces to zero away from all optical boundaries. Before 
proceeding with a more detailed examination of the two terms comprising the edge 
diffracted field, some disadvantages of the formulation in Equation (2.102) should be 
mentioned. First, the singularity near the real ray optical boundaries is subtracted 
by the curvature correction term. For practical implementations this could result in 
numerical instabilities due to the subtraction of two large numbers that tend to in- 
finity at the aforementioned boundaries, and therefore it might be necessary to take 
special limiting procedures near these boundaries. Second, the expression for the 
curvature correction term involves information about the reflection points explicitly 
contained in the vector quantities A d and B d , and at least from a philosophical point 
of view, the locality ansatz of GTD/UTD postulating that the edge diffracted field 
only depends on the geometrical properties of the surface near the edge is clearly 
violated. Nevertheless, this non-local nature of the diffracted field is a direct result 
of the rigorous Chester et al. expansion [7] employed at the outset of the asymptotic 
analysis which provides for a systematic method for obtaining the higher order terms 
in the asymptotic expansion for the fields in the vicinity of the smooth caustic. 
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Next, let us proceed with a more in depth examination of the two terms com- 
prising the edge diffracted field. The non-uniform term will be examined first. From 
Equations (2.8) and (F.9) we have 


fjk 8 a * S a * [ Wa * S a * ^(^«)] 


F(Q = 

${l a ) = - cos <p'a - cos <p a , 


2 * 


( 2 . 111 ) 

( 2 . 112 ) 


where the vector quantity in the numerator of (2.111) can be written as follows: 

3 d x s d x [n„ x x A(<2„)] = - A(Q a ) • [sin ip' a e ± e ± — sin tp a e j,? Jfj , (2.113) 

with <p' a and ip a being the angles of incidence and observation, respectively, measured 
from the edge-fixed coordinate system, and e x , e]j < ' are the unit vectors fixed in the 
incident and diffracted rays as shown in Figure 2.4. Thus, using (2.111)— (2.113) the 
non-uniform term of the diffracted field is given by 


E d n SP)~E\Q a )-D n » (Q a ) 




(2.114) 


where E'(Q a ) is the incident electric field at the edge and D„ u is the non-uniform 
dyadic PO diffraction coefficient given by 


D n u(Qa) = D.(<p' a ,<p a )exe± + D h (<p' a ,(fi a )e' n e d , 

sin 


(2.115) 




sin <p a 


DM, V*) ± (cos ^ + cos • ( 2 * 116 ) 

It is well known that the PO diffraction coefficient does not satisfy reciprocity nor 
does it satisfy the local boundary condition on the surface for the soft acoustic or TM 
polarization case. These shortcomings are expected to produce erroneous results for 
observation points away from the optical boundaries where the total diffracted field 
would equal the non-uniform PO part. A rigorous method of providing the necessary 
corrections to the PO approximation is to employ the PTD technique [19] at the out- 
set of the analysis. PTD introduces a non-uniform induced current near the edge in 
addition to the GO current. However, this addition most often makes the asymptotic 
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reduction of the resulting radiation integral quite difficult if not impossible, thus ne- 
cessitating the use of time consuming numerical integrations. A rather heuristic but 
indeed simpler approach for correcting the edge diffracted field contribution is to im- 
prove the PO diffraction coefficient by introducing a pair of multiplicative correction 
terms (one for each polarization) derived by James [20], that is 


= sin(p/2) sec(<p'/2 ) , and (2.117) 

= cos(v?'/2) csc(ip/2 ) , (2.118) 


which enforce reciprocity and the local boundary condition for the soft acoustic or 
TM polarization case so that the PO diffraction coefficient yields the exact diffrac- 
tion coefficient outside the optical boundaries. A plot of the PO correction factors 
appears in Figure 2.5. The correction factors are approximately unity near the real 
ray optical boundaries (ip = 7r ± <p'), and therefore they do not affect the cancellation 
of singularities or the compensation of discontinuities in the incident, reflected, and 
complex ray field components provided by the curvature correction term in the edge 
diffracted field. 

Next, we will investigate the behavior of the curvature correction term of the 
edge diffracted field by examining certain limiting cases. The curvature correction 
term is expected to cancel the singular non-uniform term and enforce total field 
continuity near the optical and composite shadow boundaries, while reducing to zero 
for observation points far removed from such regions. 

Limiting case la : Observation point on the deep lit side of the smooth caustic of 
reflected rays. 

In this case 'y a < -1 in Equations (2.104) and (2.105), and using the asymptotic 
formulas for the incomplete Airy functions and their derivatives in Appendix B, the 
curvature correction term assumes the following form: 


Et(p L ) 
Ei(P i) 


f"S-' 






; gcg) - 1 ~ 

. 2 i*K-7a) l/ \ 
' 1 - F’tf) ' 
2j»K-7a) l/4 . 


if m(Q 0 ) > 0, 
if m(Q 0 ) < 0, 


(2.119) 
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( 2 . 120 ) 

( 2 . 121 ) 



PO Correction Factors 



Figure 2.5: Plot of the PO correction factors with <p‘ — 45°. 
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where 


V = ± 


7.C. + C/3 + j(-7.r 


1 i 


C < ( — 7«) 1/2 S g n ( Ca ) > 


( 2 . 122 ) 


and .F(x) is the UTD transition function [4]. Note that when the observation point is 
far removed from the optical boundaries, that is |i/| 1, F(tj 2 ) — ► 1 and the curvature 

correction term reduces to zero as expected. 

Limiting case lb : Observation point on the deep lit side of the smooth caustic of 
reflected rays and close to the RSB. 

In this case, ( a « — ( — y a ) l/2 with 7 0 <C —1 in Equations (2.104) and (2.105) or tj « 0 
in (2.119) and (2.120). Again using the asymptotic formulas for the incomplete Airy 
functions and their derivatives in Appendix B, the curvature correction term assumes 
the following form: 


Et(P L ) ~ [A d — (— 7„) 1/2 2?d] | 

if m(Qo) > 0 , 

E d {P L ) ~ e- , * ( *“ +, “ ) [A d -(-7 a ) 1/2 i d ]| 
if m(Q 0 ) < 0 . 


J 

7a + (l 


~ J 

7a + C 


± 


± 


2(-7a) l/4 J ’ 

VFFe-^f 1 

2(-7a) 1/4 J ’ 


(2.123) 

(2.124) 


Now, from Equations (2.28), (2.36), (2.87), (2.106)— (2.109), (2.115), and (2.116) we 


have that 

A d -(- la y /2 B d = 


(~7a) 1/4 

'a(Q.)R(Q.) 

1 p: 


Vk J 

VPl + K 


(Qa) 


**«* |$j| 

sin ip' a R ( Q a ) 


(2.125) 
, (2.126) 
(2.127) 


\/2irA:(cos 9?^ + cos y> a ) 
and using the fact that s T a = s* at the RSB, the curvature correction term is given by 

Et(P L ) ~ -E^Qa) ■ Dnu (Qa) l;? ' 










(2.128) 


with the ± sign for the shadowed and illuminated sides of the RSB, respectively. 
Therefore, the singular non-uniform term is cancelled, and the discontinuity in the 
reflected field is appropriately compensated to provide continuity in the total field. 
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Limiting case lc : Observation point is close to the ISB. 

In this case, £ a ss (— 7 a ) 1/2 with 7 0 <C —1 in (2.104) and (2.105) or 7/ w 0 in (2.119) 
and (2.120). Again using the asymptotic formulas for the incomplete Airy functions 
and their derivatives in Appendix B, the curvature correction term assumes the fol- 
lowing form: 


e!(p l ) 


~ e 


\X I ( <y y/ 2 B 1 ! - 4 ^ 1 

[a,+( 7 .) B 4\ la+Cl ± 2(- Ia r<j' 

if m(Q 0 ) > 0 , 

^ e -J*(*k+*a) [ T _L y/ 2 R 1 / zi 4- \ 

1 a + ( 7>) d| l7. + C 2(-7.)'"/’ 

if m(<2o) < o . 


(2.129) 


(2.130) 


Now, from (2.87), (2.109), (2.111), (2.115), and (2.116) we have that 

e‘H- 7.)'" [!«?.)•/' 


A, + (— 7.)' ,, ft = 


V* 


Vk + 


Dnu ( Qa ) 


7.+C - -sVk-r.r ZM. 

V sm <Pa / V K+ S a 

e~ j * sin <p' a I 


jsb \/27rfc(cos ^ -f cos <p a ) ’ 

,-gd 

ir ii 


/ = e^ej. -f ef.e J, 


(2.131) 

(2.132) 

(2.133) 

(2.134) 


and using the fact that 4 -s* = \P'P L \ = s’ p and A(Q a ) = A(P L ) at the ISB, the 
curvature correction term is given by 


o~i k, a 


1 r. 


Et{P L ) ~ -E\Q a )- D„AQ*)-7^±7;A{P,) 


e~ ik ’p 


3 r 


= -E^(P,)±\e‘(P l ), 


(2.135) 

(2.136) 


with the 4 sign for the shadowed and illuminated sides of the ISB, respectively. 
Therefore, the singular non-uniform term is cancelled, and the discontinuity in the 
incident field is appropriately compensated to produce a continuous total field. 
Limiting case 2a : Observation point lies on the smooth caustic of reflected rays. 
In this case Q l 2 — > Q 0 thus, -y a = 0 with ( a < 0. Also, B d « 0 and using the 
small argument form of the incomplete Airy functions in Appendix B, the curvature 
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3 -i*(**a + »a) 


— e'<S'V*S-*' , r(l/3,iC/3) + i, , 
- e -«i»e->f3- , 'T(l/3,-;C/3) - j; 


(2.137) 


(2.138) 


correction term is given by 

E 4 e (P) ~ A d e 

if m(Q 0 ) > 0 , 

E d e (P) ~ 

if m(Q 0 ) < 0 , 

where T(x,y) is the incomplete Gamma function [21]. When £„ <C — 1 or the obser- 
vation point is far removed from the RSB, the incomplete Gamma function assumes 
the following large argument form: 

(2.139) 

and letting x = 1/3 and y = jCl/3 one can show that E*(P) — ► 0 . 

Limiting case 2 b : Observation point coincides with the caustic termination point 


In this case Q 12 = Q 0 -+ Q a , thus 7 * = 0 , ( a = 0 % and (2.137), (2.138) become 




_e*f 3-*'* r (i/ 3 ) + ^ 

»a . 


^4 


E*(T~) 

E d (T + ) 

Now, from (2.92), (2.93), and (2. 107)— (2. 109) we have 

A d = 


e~ i6 t 3~ 2/a T(l/3) + ^ 

’a J 


e^(fc/ 2) 1/6 cos 0 ; 


|m(Q 0 | l/a 

. 

V< *1 

J 

= — 6 fc J/s (cos + cos y> a ) 

2 

m(Qo) 

1 

> 


D„u (Qa) 


sin <p' a R ( Q a ) 


(2.140) 

(2.141) 

(2.142) 

(2.143) 

(2.144) 


It %/27t k (cos <p' a + cos tp a ) ’ 

and using the fact that s*;* = s' Q ' r and cos 0' o = sin <p' a near T, the curvature correction 

term assumes the following form: 

'cosfl‘ (fc/6) 1/6 r(l/3)e^' 


E*{T*) ~ -E\Q a ) -D nu (Q a )^^± 


\/3 


v^HQo)| ,/J 


x E\Q 0 )- R (Q 0 ) 


e~ ik ‘o 




(2.145) 

(2.146) 
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where E r (T ) is given in (2.96). The curvature correction term cancels the singular 
non-uniform term and also compensates the discontinuity of the reflected field to 
ensure continuity in the total field. 

Limiting case 3a : Observation point on the deep dark side of the smooth caustic 
of reflected rays. 

In this last case, 7 a > 1 and using the asymptotic formulas for the incomplete Airy 
functions in Appendix B, the curvature correction term assumes the following form: 

Et(P s ) ~ «g»(C.) (A + C.A) | -J, ~ } ■ (2.147) 


if m(Q 0 ) > 0 , 

K{P S ) ~ sgn 

if m(Q 0 ) < 0 , 


Vo/ v , 

( C)e - J .ci..i, (i <+ (.A) | • ( 2148 > 


where 


( = sgn(Ca) 


7aCa + C/3 + Jg7a /2 


1 i 


(2.149) 


and F(C is the Fresnel transition function of complex argument defined in Ap- 
pendix E. When the observation point is far removed from the CRSB on either side, 
3i(£e W4 ) > 0 with \(\ > 1 and F(() — ► 1, thus the curvature correction term reduces 
to zero. 

Limiting case 3b : Observation point near the CRSB. 

In this case, = 0*, &(£e W4 ) « 0, F(£) assumes the following form: 

F(t) = e i(1 +1, (2.150) 

and the curvature correction term in Equations (2.147) and (2.148) becomes 

(2.151) 


E d (P s ) ~ sgn(Ca )^ A d 

3/ 0^1/4 d 


2 71' 

Now, from Equations (2.44) and (2.76) with C, a = 0, (2.100), and (2.110) we have 

2 


•*. + < = + 7i 

A d ■■ 


S/2 


hi" 

'AQU)-R(QU) 

I Pt,2 


yp5 

V pia + 


(2.152) 

(2.153) 
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and the curvature correction term near the CRSB is given by 


E‘.(Ps) ~ 


jsgn((. )E‘(Q;,)-i(Qi, 
isgn(C„)£‘(P s ), 



c--> k ‘u 


m (Qo) >o , 


(2.154) 

(2.155) 


where ^‘(Ps) is the complex GO ray field. Since ( a < 0 in the illuminated side of 
the CRSB and £„ > 0 in the shadow side of the CRSB, the curvature correction term 
indeed compensates the discontinuity in the complex ray field across the CRSB. 


2.2 Derivation of the UGO/EUTD solution using 
the method of steepest descent 

In this section, the UGO/EUTD solution for the reflected, complex ray, and edge 
diffracted fields is derived using the method of steepest descent [8] for the asymptotic 
reduction of the PO radiation integral representation for the scattered field formu- 
lated in the previous section. We begin with the expression for the scattered field in 
Equation (2.26), that is 

E‘(P) « e~ iha 1(0, s Q ; k ) , (2.156) 

where a is given in (2.18)-(2.20) and l(0,s a ',k ) is a stationary phase integral in the 
s-plane given by 

oc 

1(0, s a] k) = j G(s) e-' w * ,/s) ds , (2.157) 

•a 

where 0, s a , and G(s) are given in (2.21 )— (2.23), (2.25) and (2.27), respectively. Only 
the case of m(Q 0 ) > 0 will be considered here since the results for m(Q 0 ) < 0 can be 
easily obtained by letting 0 = 0exp(j2ir/3 ) and s a = s a exp(jir/3) in (2.157). The 
complex plane topology of the stationary phase integral in (2.157) for the lit side 
of the caustic (0 < 0) is shown in Figure 2.6. The original path of integration, P 0 , 
running from s a to oo is deformed into the appropriate steepest descent paths through 
the endpoint and the two real saddle points s 12 = ±(— /?) 1/2 , and Equation (2.157) 
may be written as follows: 

1(0, Salk ) = u(s 3 -s a )[ G(s)e-w+- 3 '*ds 

21 
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Figure 2.6: Complex plane topology of the stationary phase integral I((3,s a ;k) for 
observation points in the lit side of the smooth caustic (/ 3 < 0). 
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Re (s) 



Figure 2.7: Complex plane topology of the stationary phase integral I((3,s a ]k) for 
observation points in the dark side of the smooth caustic (/ 3 > 0). 

+ u( Sl -s a )f G(s)e- ik ^- S/ ^d3 

^13 

+ f G(s)e--’* (s,+,5/s) ds, (2.158) 

JPn 

) 1 if <S 2 < 3 a < Si 

2 if s a < s 2 (2.159) 

3 if s a > 3, . 

In a similar fashion, in the dark side of the caustic (/? > 0), the original path of 
integration P 0 is deformed into the steepest descent paths through the endpoint and 
the complex saddle point s 2 as shown in Figure 2.7 and Equation (2.157) may be 
written as follows: 


/(£,*.;*) = «(-*.)/ G(3)e~W+‘ 3 '*d3 

-'J'ZS 

+ f G(3)e- j W a+ti Mds, 

JPttx 


(2.160) 
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m 


< 


2 if s a < 0 

3 if s a > 0. 


(2.161) 


Next, it is assumed that the vector amplitude function, G(s), is regular and slowly 
varying along each of the paths of integration in Equations (2.158) and (2.160), and 
since the principal contribution from each integral would come from the neighborhood 
of either a saddle point or an endpoint, (?(s) can be expanded in Taylor series near the 
pertinent critical point for each integral. Thus, using the definitions for the ordinary, 
Fock-type, and incomplete Airy functions in Appendices A and B, (2.158) and (2.160) 
may be written as follows: 


where 


I(0>S.ik) ~ G(s 2 )W 2 (k 2 ' 3 p)u(s 2 - 3 .) 

+ k- 1 ' 3 G(s a )g;(k'' i l3,k l ' i s a ) + 0(k- 1 ) if /3 < 0, 


1 if s a > s, 

2 if s a < s 2 

3 otherwise, 


l(0,s a ;k) ~ 2ivk- 1/3 G{s 2 )Ai(k 7 ' 3 p)u{-s a ) 


+ k- t ' t dMg;(k ,ft 0,V''s a ) + O(k-') if /3 > 0, 

1 if s a > 0 

2 if a a < 0, 


(2.162) 


(2.163) 


G(s ul ) 

= F ^ r {klv 

(2.164) 

GM 

= HQ (0+ ■]' 

(2.165) 


[jk cos ei ,A(Q lt7 )- R ($,.,) 

V sJs\A 

(2.166) 

^(^1,2) 

= cos,e; ’(^ + ^)’ 

(2.167) 

HQ 

HQ 

jkA(Q a )‘ Dn u (< 2 a) 

(2.168) 
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R (C?i,j) is the dyadic GO reflection coefficient given in Equations (2.60) and (2.61), 
and D nu ( Q a ) is the dyadic non-uniform PO diffraction coefficient given in (2.115) 
and (2.116). The first two terms in Equation (2.162) are the uniform contributions 
from the two real stationary phase points associated with the reflected field in the 
lit side of the caustic. The first term in (2.163) is the uniform contribution from 
the complex stationary phase point associated with the complex ray field in the dark 
side of the caustic. Finally, the third term in (2.162) and the second term in (2.163) 
are the uniform endpoint contributions associated with the edge diffracted field. The 
main advantage of (2.162) and (2.163) over the formulation of the previous section is 
that the expression for each term contributing to the integral is independent of the 
others, and also the correction to the non-uniform edge diffraction coefficient near 
the optical boundaries is multiplicative rather than additive. These very desirable 
features, however, come at the expense of a solution that does not contain the higher 
order terms involving the derivatives of the canonical functions. The derivative terms 
in the expressions of the previous section implicitly include slope scattering that oc- 
curs when the source pattern has a null in the direction corresponding to observation 
points near the caustic or the composite shadow boundary. It is believed, however, 
that when slope scattering effects are not important the higher order terms are not 
significant in terms of solution accuracy because they are only necessary for recov- 
ering the well known non-uniform expressions for observation points away from the 
transition regions. As will be shown later in this section, the UGO/EUTD expressions 
for the reflected, complex ray, and edge diffracted field also recover the non-uniform 
expressions outside the transition regions even though the derivative terms are ab- 
sent. Furthermore, the results of Section 2.1 may be used to derive slope reflection 
and diffraction coefficients to be used in conjunction with the regular UGO/EUTD 
expressions. 
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2.2.1 UGO expression for the reflected field 

The UGO expression for the reflected field is as follows: 


V Pi + 3 l 

V P5 + *2 


(2.169) 


where = 1 in regions I and II of Figures 2.2 and 2.3 and zero elsewhere, and 
u 2 = 1 only in region I and zero elsewhere. The quantity % (Q 1>2 ) is the dyadic UGO 
reflection coefficient and is given by 


n(Qui) 

= R'.iQ 1,2) e±e± + R-h{Qia) * Jj , 

(2.170) 

R^t,h(Q 1,2) 

= R, h e-'f W^-7^) , 

(2.171) 

R.,h 

= =pl , for a PEC boundary. 

(2.172) 


A plot of the scalar UGO reflection coefficients, appears in Figure 2.8. Equa- 
tion (2.169) shows that the complex conjugate of the UGO reflection coefficient should 
be used when the spread factor p r /{p r + «s r ) < 0. The parameter 7 12 depends on the 
relative phase distance between the reflection point and the point Q 0 associated with 
the smooth caustic and is given by 


7l,2 — 


§‘ 


(s’ + Ro) - (sj j + s[. 


(2.173) 


The caustic point Q 0 is found from the following relation: 

^R(Qo) - s‘(Qo)] *n(Q 0 ) [?’(Q 0 ) -n(Qo)] 2 [^(<2o) • «(Qo)] . . 

P t (Qo) + s<(Q 0 ) + R(Q 0 ) ~ (2 - 174 ^ 

which in general must be solved using a computer search procedure. Notice that in 
the case of plane wave incidence with the observer in the far-zone of the scattering 
boundary, Q 0 is simply an inflection or zero-curvature point. 

When the observation point is in the deep lit side of the smooth caustic, 7, 2 1 

and using the large negative argument form of the Fock-type Airy function, R, A — ► 
R., h = ^1 which is the usual GO reflection coefficient for a PEC boundary. When the 
observation point lies on the smooth caustic, the reflected ray spread factors in (2.169) 
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0.1 

7 


Figure 2.8: Plot of the scalar UGO real ray reflection coefficients. 
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become singular where the parameter 7 12 — > 0 in (2.171), and using Equations (2.28), 
(2.31), and (2.36) 71, h assumes the following form: 


tt. lh (Qo) = kiR^e-t* W,(0) |m r (Q„)|-* cos0;y^^ , (2.175) 


<n\ fl i >/ n \ 3 cos 20 o sin ^o 

m r {Q 0 ) = cos 6 0 K a (Q 0 ) 


1 


3 cos 2 8^ sin 8 r 0 

27: 


1 

2^ L?„(Qo)cos0’ 3’J 

1 


1 

+ — 


(2.176) 


[p a {Q°) cos 6i 

Hence, the singular spread factors are cancelled and the total reflected field remains 
finite and agrees with the expression in (2.96) derived in the previous section. 

2.2.2 UGO expression for the complex ray field 

The UGO expression for the complex ray field in the dark side of the caustic desig- 
nated as region IV of Figures 2.2 and 2.3 is as follows: 


E'{P s )~E'(Q\,)-n c {QU)< 


PZ 


p[% + 


e~ ik ‘Zi ; m(Qo) > 0 , 


(2.177) 


where p" 2 /(p" 2 + s" 2 ) < 0 in this case and H‘(QU) is the UGO dyadic complex ray 
reflection coefficient given by 


n c {Qu) 

KAQU) 


Kh = 


K:{Ql 2 )Z ± e ± + K‘ h (Ql 7 )e J-ep, 

e (2/3b * /2 Ai(7.) , 

=pl , for a PEC boundary. 


(2.178) 

(2.179) 

(2.180) 


A plot of the scalar UGO complex ray reflection coefficients, 7V t h , appears in Fig- 
ure 2.9. The parameter 7, depends on the relative phase distance between the complex 
reflection point Q\ or Q\ and the point Q 0 associated with the smooth caustic and is 
given by 

p o n f- 

(2.181) 


7. = 


j* 


(»; + &)- (»« + 


- 21 
a 


When the observation point is in the deep dark side of the smooth caustic, 7, 1 

and using the large positive argument form of the ordinary Airy function, 7t c t h — * 
R\ h = =Fl which is the complex GO reflection coefficient for a PEC boundary. 
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2.2.3 EUTD expression for the edge diffracted field 


The EUTD expression for the edge diffracted field is as follows: 


o~i k ’a 


E*(P)~E'(Q a )-V(Q a ) 


where V (Q a ) is the EUTD dyadic edge diffraction coefficient given by 


T>{Q a ) = V,(v' a ,<p a )e ± e ± + V h (<p' a ,'p a )e\.e* 




-e-*t 

2\/2Trk 


cos 


(*4*) cos(^)J 




^ r (7 a ,Ca) is the EUTD transition function given by 

E{h() - i(7 + Og;(7i() ej(lC+<>/1) , 

1 if C > (-t) 1/2 

2 if C < -(~ 7 ) l/2 

3 otherwise 


i = < 


I 


; for 7 < 0, or 


; for 7 > 0. 


1 if ( > 0 

2 if C < 0 

The parameters 7 0 and ( a are given by the following expressions: 


(cos^ + cospo) 

IMV , a>y J a)| l/S 

u a = sgn [L d (<p' a ,<p a )}, 

/3 yj\ cos v >' 0 + cos <p 0 - cos y>[ [ - cos p a | 

S a ~ tin re 


ti ( 


ti 0 = < 


+1 if Qo ^ C 
-1 i{Qo£C, 

where L d (<p' a ,ip a ) is the EUTD distance parameter and is given by 

t / , \ 1 , . , • v . , _ . COS V' + COS 'Pa 

«) = -( Sln <fi a + Sln <Pa)K„(Qa) + 2 ?(Q~) 


+ 

+ 


3 sin <p' a cos ip' a 

3 sin <f a cos <p a 

2 ? 


1 + sin ¥ > l 


lPg{Qa) 8' a 

1 + sin <p a 


Pl(Qa) 


(2.182) 

(2.183) 

(2.184) 

(2.185) 


(2.186) 

(2.187) 

(2.188) 
(2.189) 


(2.190) 
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Figure 2.11: Plots of the real and imaginary parts of the EUTD transition function 
in the second quadrant of the 7 ^-plane. 
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It is understood that when L d (<p' a ,<p a ) < 0, the EUTD transition function in (2.184) 
should be replaced by its complex conjugate version. Plots of the real and imaginary 
parts of the EUTD transition function in the four quadrants of the 7^-plane appear 
in Figures 2.10- 2.13. One can observe that the EUTD transition function exhibits 
three lines of discontinuity. The first discontinuity occurs along the £ = 0 line which is 
the boundary between the first and fourth quadrants of the 7^— plane, and physically 
corresponds to the CRSB. Figures 2.10 and 2.13 show that only the imaginary part of 
the transition function is discontinuous across the CRSB while the real part remains 
continuous. Also, as 7 > 1 the discontinuity diminishes. This transition function 
behavior makes sense physically since 7 1 corresponds to the deep dark side of the 

smooth caustic where the complex ray field becomes vanishingly small and accordingly 
its discontinuity across the CRSB should also diminish. The second discontinuity of 
the transition function occurs along the ( = (—7 ) 1/2 line in the second quadrant of the 
7^-plane as shown in Figure 2.11. This line physically corresponds to the ISB. Finally, 
the third discontinuity of the transition function occurs along the £ = — (— 7) 1/2 line 
in the third quadrant of the 7^-plane as shown in Figure 2.12. This line physically 
corresponds to the RSB. 

Next, the behavior of the EUTD transition function for observation points in the 
deep lit side of the smooth caustic will be examined, so that a comparison between 
the EUTD edge diffraction coefficient and the original UTD expression [4] can be 
made. In this case the EUTD transition function assumes the following form: 

~ (-t. + C) + 1 ■ ( 2 - 191 ) 

where the parameter rj is given in Equation (2.122). When the observation point is 
away from the optical boundaries, |7j| > 1, the UTD transition function F(rj 2 ) — ► 1, 
and „,C<.) 1- Hence, the EUTD edge diffraction coefficient assumes the ex- 

act expression outside the transition regions associated with the optical boundaries. 
When the observation point is near an optical boundary, we have that 7 a + c » 
2 tj ( — y a ) I/< and therefore T{ 7 a ,C<>) ~ F(t] 2 ), with 77 reducing to the following expres- 
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sion: 


V 




1 / 2k 

--(cos^ + cosv,.)^^ , 

sin <p' a + sin y? a sin 2 sin 2 y? a 

P.(Q.) < *2 


(2.192) 

(2.193) 


Near the RSB, <p a « 7r - ^ and thus 


cos’ ( V> 2 —) . “i 

’’’ K 


(2.194) 

(2.195) 


Near the ISB, <p a « 7r + ^ and thus 

*(*>>.) * (jr + ji) «*’ (^) . < 2 - 196 ) 

t/ 3 w 2k ( cos 1 ( ^° 2^° ) = — p a ) , (2.197) 

where L i r are the UTD distance parameters. Both L' T and a(y>'±y>) are defined in [4]. 
The above results indeed show that near the optical boundaries the EUTD transition 
function is identical to the UTD transition function. Therefore the two diffraction 
coefficients should give equivalent results despite the fact that the correction near the 
optical boundaries is applied differently in the two methods. In the UTD case, two 
transition functions with different arguments are used, one for each of the two terms 
in the non-uniform diffraction coefficient. In the EUTD case, only one transition 
function is used and operates on the entire non-uniform diffraction coefficient. 
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Chapter 3 

Near-Zone Uniform Asymptotic 
Analysis and UGO/EUTD 
Formulation for the Scattered 
Fields near a Three Dimensional 
Smooth Caustic and Composite 
Shadow Boundary 


In this chapter, the behavior of near-zone high frequency scattered EM fields near 
a three dimensional smooth caustic and composite shadow boundary is investigated, 
and a UGO/EUTD solution is developed for the description of the scattering mech- 
anisms. The analysis for the 3-D case is very similar to the 2-D case of Chapter 2. 
The three dimensional surfaces considered in this chapter are assumed to be rotation- 
ally symmetric or slowly varying in the transverse direction, with the 2-D boundaries 
of Figures 1.1 and 1.2 being the generator curves representing one of the principal 
coordinate curves of the 3-D surface. 

In order to further simplify the asymptotic analysis, it is also assumed that the 
plane of incidence defined by the incident ray direction and the normal to the surface 
at a real reflection point coincides with one of the principal coordinate curves of the 
surface. Although this is a very special case in a general scattering problem, the 
derivation of key expressions in the asymptotic analysis such as the reflected ray 
spread factors will be kept as general as possible. Because of the highly local nature 
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of high-frequency fields, the results for this special case will remain valid for general 
3-D boundaries that are not necessarily rotationally symmetric with the plane of 
incidence not necessarily coinciding with one of the principal coordinate curves of the 
surface. 

The PO approximation is also employed for the 3-D case. The problem formu- 
lation and the derivation of a uniform asymptotic solution that remains valid near 
a three dimensional smooth caustic and a composite shadow boundary is derived in 
Section 3.1 using the Chester expansion for the asymptotic reduction of the PO ra- 
diation integral. In Section 3.2 the method of steepest descent is employed for the 
asymptotic reduction of the PO radiation integral and the UGO/EUTD expressions 
for the reflected, complex ray, and edge diffracted fields are derived. 

3.1 Derivation of a uniform asymptotic solution 
using the Chester expansion 

The electric field E s which is scattered by a 3-D surface containing an edge illumi- 
nated by a ray optical field due to a point source at P ; and observed at a point P 
exterior to the surface as shown in Figure 3.1 can be expressed in terms of the usual 
radiation integral over the induced current J on the surface S as follows: 

E‘(P) « ^ J/[RxRx /(<?')] ^ is 1 , (3.1) 

where Q' is any point on 5, R = Q'P, R = |i2|, R = R/ R, and P cannot be too 
close to S for (3.1) to be valid. As in the 2-D case, the PO approximation to (3.1) 
is employed using the GO currents of Equation (2.2), where n' is the outward unit 
normal vector to S at Q'. The incident magnetic field at Q' due to the point source 
at P' is expressed ray optically as follows: 

H'(Q') ~ Z~ l (?' x A) ; r-A = 0, 

= zprxEW), 

= A{Q') —jr~ 
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(3.2) 

(3.3) 

(3.4) 
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where A is the vector source pattern, s' = P'Q', s' = and s' = s' /s'. Thus, 
using the PO approximation Equation (3.1) is given by the following stationary phase 
integral: 

oo oo 

E‘{P) « J J F(t, r') e - ik « t,y) dl'dr ' , (3.5) 

la ~ oc 

where 


F{1',t') = 


jk R x R x [n' x a' x A(l',r ') ] 
27t s' R 


+ r 1 ), 


(3.6) 

(3.7) 


and t, t 1 are the principal coordinate curves of the 3-D surface S as shown in Fig- 
ure 3.1. It is also assumed in (3.5) that the edge is located at and therefore 

it extends along the principal coordinate curve t' which is always orthogonal to t 
and assumed to have a slowly varying radius of curvature. In order to simplify the 
analysis, it is further assumed that the plane of incidence defined by the incident 
ray direction and the normal to the surface at the reflection point coincides with 
the principal coordinate curve l' on the surface. Although this situation occurs only 
in special cases, it is chosen because it greatly simplifies the asymptotic reduction 
of (3.5). Nevertheless, due to the highly local nature of high frequency fields, the 
results for this special case can also be applied to general 3-D boundaries that are 
not necessarily rotationally symmetric and the plane of incidence does not necessarily 
coincide with one of the principal coordinate curves of the surface. 

The phase function in (3.7) contains two real or complex stationary phase points 
corresponding to real or complex reflection points Q y and Q 3 on S that are arbitrarily 
close together or to an edge diffraction point Q a as in the 2-D case of Chapter 2. The 
stationary phase condition is given by 

V • t') = 0 and T ' . t') = 0 , (3.8) 

where l \ t' are the principal surface directions, and both equations in (3.8) must 
be satisfied simultaneously at the stationary points. With the assumption that the 
plane of incidence coincides with the principal coordinate curve l' , Equation (3.8) 
yields stationary phase points for which r' = r, along with t = i x for Q x and t — l? 
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for Qi. Therefore, the coordinates of the stationary phase points are = (£, , r.) and 
Q 3 = with the edge diffraction point coordinate being Q a = (^„,t,). Then, 

using a local quadratic approximation, the phase function in (3.7) may be written as 
follows: 

W, r') ~ Ht, t.) + (r' - T.)4,.(t,T.) + i(r' - r.)V„(<',r.) , (3.9) 

where the notation ^ = <f>zv{z\y') was used in (3.9). Upon inserting (3.9) in 

Equation (3.5) we get 

oo oo 

E‘(P) w J dt e • J dr' F{1‘, r ') , (3.10) 


la 


where 


,,,, n - Mt> rU 

i T >) — i T >) -7 


24> Tr (t',T.) 

= r.- 


(3.11) 

(3.12) 

(3.13) 


2* Tr (/',r.) 

When the plane of incidence coincides with the principal coordinate curve t ', <£ T (£', r.) = 
0, V and thus rp(i',T.) = and (( 1',t ,) = t„ V However, we will keep 

all expressions as general as possible since we would like our results to remain valid 
for an arbitrary plane of incidence, in which case <j> r (l',T.) = 0 only when t = t 1>3 
from the stationary phase condition. It is also assumed in (3.11) and (3.13) that 
<j> TT (£',T.) ^ 0 VI'. An additional necessary condition is the matching of the phase 
function when one or both reflection points coincide with the edge diffraction point, 
that is 

V'(Cr.) = W. 2 ,r.). (3.14) 

Prom (3.11) Equation (3.14) is satisfied only if j> T (l a ,T.) = 0, and using (G.81) this 
edge condition can be expressed as follows: 


(?:-«!)•?. = o, 


(3.15) 


where ?„ is the unit edge tangent vector at Q a . Equation (3.15) can be used to find 
Q a and is recognized as Keller’s law of edge diffraction. An additional consequence 
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of Equation (3.15) used in (3.13) is that 


t{la,T.) = T,. 


(3.16) 


Next, the inner integral in (3.10) denoted as 1,(1* ,r,) can be evaluated using the 
stationary phase method [8] which yields 


fa F(f,{) 

~ J-Tj- , 


(3.17) 


Hence, employing (3.17) in (3.10) the 3-D analysis effectively reduces to a 2-D situ- 
ation with E*(P ) given by 


oo 

E‘{P) * j d£' G(£\t.)< 

la 


,r s ) 


where 


C(i> T \- [jk Rx[n' xs i x A(4,0] 

,T ' V 2ff 


(3.18) 


(3.19) 


and the phase function 0(£',r,) is given by Equation (3.11) and has stationary phase 
points at £' = A, 2 corresponding to Q , i2 , that is 0(4,2, 7-.) = 0. Therefore, the asymp- 
totic reduction of the integral in (3.18) would follow the same steps as the asymptotic 
reduction of the integral in (2.7) for the 2-D case. 

The phase function V'( 4,r ,) is then expanded in a Taylor series around the point £ 0 
which is a stationary point of the derivative of the phase function, that is 0(£ o , r.) = 0 
and must necessarily lie between the two stationary phase points Thus, using a 
local cubic approximation, ip(£',T,) may be written as follows: 


0 ( 4 , 7 -.) ~ ip{lo,T.) + (£' - £ 0 )4>(£ 0 ,t.) + 
(t'-lo) 


(4 - 4) 5 

2 ! 


0(4, 7-.) 


+ 


3! 


0 (4,t,). 


(3.20) 


A phase matching condition is also required at the caustic point Q 0 for the case when 
4,2 — > 4 corresponding to observation points at the smooth caustic; therefore we have 


'0(4,7 -,) — 0(4,2, 7-,), 


(3.21) 
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which is satisfied only if fa(l 0 i t, ) = 0 in Equation (3.11). Using Equation (G.19) this 
caustic condition can be expressed as follows: 

(i2o — s') • f 0 = 0 . (3.22) 

Then, using Equations (3.11) and (3.22), fat 0 ,r,) and its derivatives in (3.20) are 
given by: 

(3.23) 

(3.24) 

(3.25) 


fato, 

T.) 

= S'o + Ro , 


kto, 

T.) 

— fao(Qo) , 


kt 0. 

T ‘) 

= fao(Qo) ~ 

fax 

fa? 

i (to, 

T .) 

— fao(Qo) — 

fax 

fa? 



= m(Q 0 ) , 


4 >mn{Q 0) 

A d m+n fat, 

dl ,m dT' 

fa) 

n 


= 0 , 


fax(Qo)fa?(Qo) 


(3.26) 

(3.27) 


(<o. T *) 


Using (G.18), (G.19), (G.35)-(G.37), and (G.59)-(G.62) the derivative functions fa 0 , 
fa o, fa n fa„ fao, fax , and fa? at Q 0 = (t 0 ,T.) can be expressed in terms of the geometric 
parameters of the surface S at Q 0 , that is 

(?; - Ro) - to, 

( ~ ~ \ ~ 1 — ( a 0 ’ 4 >) 1 — (Ho ' to) 

faQo) (Ro - s’ 0 ) • fio -\ — 1 

a o 

(« i-to) (*’o -n) (Ro ■ to) (Ro ■ To) 

Ro 


^io(Qo) — 

<j>2o{Qo) “ 


R« 


Si 


^30 (Qo) = 
+ 

+ 

faiiQo) — 


m\(n ~ 1 1- ( 5 o- f< 0 2 

« 2 (Qo) [Ro - a 0) ’ n o d • 1 


»o 


Ro 


(3.28) 

(3.29) 

(3.30) 

(3.31) 


(A - K) ■ {"!(<?.)<. + [£ • V«,(Q.)] n„} 

1 


*0 'to 




3K:(C?o)(5o -n 0 ) - 


s' n 


Ro-t 0 

Ro 


/ ^ \ 1 (-Ro * t 0 ) 

3k,(Qo) (-Ro • n 0 ) + — 

1 — 3 (a J, • 


(3.32) 


K l(Qo) ( S 0 ’ n o) 


*0 
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+ 


fiuiQo) — 


+ 


Ro -T 0 
R o 

R 0 •/„ 

iZc 


_ l-S^-^o) 3 

^(Qo) (Re ■ Wo) + ^ 


(3.33) 


1 - 3 (s’ -fo) 21 


K'liQo') ( 5 o ' »o) 




/ ^ \ 1 3 (Ro ' To) 

M<?o) (i^o ■ no) + ^ - 


(3.34) 


where £ 0 , f 0 are the principal surface directions at Q 0 > n 0 is the outward unit normal 
at Q 0 , and k 1<3 {Q 0 ) = l/R.^Qo) are the principal surface curvatures at Q 0 along 
the £ 0 and ? 0 directions, respectively. The caustic point Q 0 is therefore found from 
Equations (3.25) and (3.29)— (3.31), that is 


^2o(Qo)<A>2(^?o) ^n(Qo) — 


(3.35) 


and must satisfy the caustic condition in (3.22). Notice that in the case of plane wave 
incidence with observation in the far-zone of S’, that is s*, Ro — » oo, Q 0 is a stationary 
point on a zero curvature line and satisfies Ki(Q 0 ) = 0 along with (3.22). Also worth 
noting is the fact that when the plane of incidence coincides with one of the principal 
surface directions, the cross- derivative functions <^u(Qo) and either 4>?\{Qo) or ^ 12 (^ 0 ) 
reduce to zero. This is also true for the plane wave far-zone scattering problem where 
all three cross- derivative functions reduce to zero for any arbitrary plane of incidence. 

Next, making the usual linear transformation, the stationary phase integral of 


Equation (3.18) transforms to the following expression: 

11 00” 

E‘(P) « e - J *“ J g{s)e- ih((,,+ ‘ 3 ' 3) ds , (3.36) 

• a 

where 

a = s' 0 + Ro , or (3.37) 

= \ 1M + <) + M + »;)] . or ( 3 - 38 ) 

= •[,+•» ±|(-r. or (3.39) 

= < + (3.40) 
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0 = *o •(?■-£>) 

3 


*, i 


or 

, 2/3 


or 

2/3 


[™(Qo)\ 

( 3 W/3 

= -{4 [W + O-W + OJj . 

= - {| [(•; + *.)-«,+ . 
= ±\A - (»i — *;) -&•(»; -A.) I 


I 


(Q «)J 


+ if Qo 3 S 
- if Qo e s, 


G(s) = G(t,r.) 

dt 

da 


dt 

da ’ 
, 1 


(3.41) 

( 3 - 42 ) 

(3.43) 

(3.44) 

(3.45) 

(3.46) 


M<WJ 

and the proper branches for /?, a a , and dt' /da depend on the sign of m(Q 0 ) given 
in (3.26). Also, “ 00 ” = 00 exp[— jaig(dt/da)\. The asymptotic evaluation of (3.36) is 
then carried out by employing Chester’s expansion [7], and following the same steps 
as in the 2-D case of Chapter 2 the scattered electric field is given by 


E’(P) ~ 

e~ ika 

y/k 

o 0 Ai’(7, ( a ) + jb 0 Ai*(7, C 

.) 

+ 


i.r,) 

' G(t a ,T.) 
jkrp(l a , t,) 

3 (a 0 + Cab 

+ V&(7 + C 

°) 

:). 

E‘(P) ~ 

g-J*a 

y/k 

a 0 Ai(7 ,( a )-jb 



-L 


i,r B ) 

1 

1 

<-o. 

+ 

r* 

C5-J 

o)l 

1 

c 


jki>(t a ,T.) 

(7 + c. 



if m(Q 0 ) > 0, (3.47) 


if m(Q 0 ) < 0, (3.48) 


where Ai( 7 ,£ a ) is the incomplete Airy integral defined in (2.48), 7 = k 2/3 fl, 
C — k 1/3 a a , and a, a 0 , and b 0 are given by the following expressions: 

q = •:+•:- sh - 1 ( 7 c. + c/3) , ° r 

S = Sgn[m( (?„)), 

_ (- 7 )** 


a 0 

bo 


G(t, r.) 


|W,r.) 


± G(t 7 ,T.) 


\ v»( 4 ,r.) 


(3.49) 

(3.50) 

(3.51) 

(3.52) 
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From Equations (3.13) and (3.19), G(£ 1 i 2 ,t,) may be written as follows: 


(3.53) 


(3.54) 

(3.55) 


sin _ n _ /i* cos o[,7 A (Qi*) • -R (<?>.*) 

L r {l l , 2 ,T.)-y r— 5 

V S U Sr i,7\j \<t>rr{l U 2,T.)\ 

where R (Q 12 ) is the dyadic GO reflection coefficient given by 

R{Qi,i) = i2.e , x (Q 1 , 2 )el(Qi, 2 ) + J R h e;|(Qi. 2 )eI(C?i. 2 )> 

.ft, ifc = ^1; fora perfectly conducting boundary. 

From Equation (3.11), t})(£ 1i2 ,t,) is given by 

" ( C _ ^ ^tt(A, 2 J 7 ’»)^tt(^1, 2 J T») 7 "*) 

V’(*l,2 5 7 '.) — 2 ft \ 

rrr \*T ,2 > *J 

Using Equation (G.69), (3.56) can be written as follows: 

W...-0 = TTTT) (t + sM (?- + y-) ’ (3 ' 57) 

Pt\ y l2/ \ 5 1,2 P rl ,r2 / 

where pj r are the radii of curvature of the reflected wavefront given in (G.70) 
or (G.71), and 6[ 2 is the usual angle of incidence (between the incident ray and 
the outward normal of S at Q 12 ). Thus, using (3.53) and (3.57) a 0 and b 0 may be 
written as follows: 


(3.56) 


4l 

y/k e j *(— 7)** 

’i(Qx)- ^(Q.) 

/ />i, p;i 

. *<■ J 

’ “ 2y/f 


V (#»« + *i)(p;i + *:) 


!(<?>)• S «,) 

/ *4 />; 2 




it 

V (p« + 

+ «;) 



(3.58) 


Also, from Equations (3.11), (3.15), (3.16), (G.18), (G.20), (G.82), and (G.83) we 
have that 


i’i.L ,t.) 

G(L,r.) 

Ha-Si 


(a 1 - ? a) • £ = - sin i? a (cos <p' a + cos <p a ) , 


[jit MQ«) • [(”* -^De^ei + (n, 


27f <j> rr (t a ,T.) 



— sin t? a sin tp' a , 


(3.59) 

(3.60) 

(3.61) 

(3.62) 
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n a -s d 

= sin i? a sin ip a , 

(3.63) 

4>rr(ta,T.) 

= sin ’ A (VJi)' 

(3.64) 

1 

Pi 

3^ R a sin 2 ’ 

(3.65) 


where t? a is the angle of diffraction as shown in Figure G.3, <p' a and <p a are the angles 
between l a and the projected incident and diffracted ray directions in the plane per- 
pendicular to ? 0 as shown in Figure G.4, e'V and e’^ d are the unit vectors fixed in the 
incident and diffracted rays as shown in Figure G.3, p d a is the distance between the 
caustic at Q a and the second caustic of the diffracted ray, and R a > 0 is the radius 
of curvature of the edge at Q a . 

Finally, using (3.59)-(3.64) the scattered electric field is given by 
E‘(P) ~ e -j=- a 0 Ai-( 7 ,C a )+i6o^Ai*( 7 ,C) 

+ {^'(Qa)- d„ (Qa) 

+ \ - + j£ 1^*% if m(Q 0 ) > 0, (3.66) 

Vk T + C J J 

E’(P) ~ a 0 Ai( 7 ,C) - jh 0 ^Ai( 7 ,Ca) 

+ ^-^(Qa) • Dnu (Qa) 

| , if m(Q 0 ) < 0, (3.67) 

where D nu (Q a ) is the non-uniform PO dyadic edge diffraction coefficient for oblique 
incidence and is given by 


J_ P - j**i °° + 

Vk [7 + Ca 2 . 




Dnu {Qa) = D.(<p' a ,<p a ' 1 # a )e i ± e d ± + D h (^l, v> a ; t? a ) e e J , (3.68) 


smy> a 


D.AVa, ± ^ b (cOS (fi' a + COS <f a ) ' 


(3.69) 


The reflected, complex ray, and edge diffracted field components in (3.66) and (3.67) 


can then be separated as usual using the properties and complex plane topology of 
the incomplete Airy integral in Appendix C. 


64 



3.1.1 Reflected field solution on the lit side of the caustic 


From (3.66) and (3.67) and using the appropriate form of the incomplete Airy integral 
in Appendix C, the reflected field on the lit side of the caustic is given by 

e t (Pl) ~ e ~ Jkai wi(-7*)] 

- u 2 e- J *° 2 [An 1 /* W a (-7i) - jBtlT 11 * W;(-70] } . ( 3 - 70 ) 

if m(Q 0 ) > 0 , 

e'(p l ) ~ [A«7; ,, w;(-7,)+js<7; ,, *w"(-7.)] 

- u je -'*-[i«7i /, W;(- 7 ,)+jB,7,- , "W"(-7.)]}. (3.71) 

if m(Qo) < 0 , 


where Wi,2 are the Fock-type Airy functions [16], u 1>2 , 7 /, and a 1<2 are the same as 
in (2.84)-(2.86), and A t , B t are given by 


A t 


> = 


A(Qi)-R(Q>) 


s\ 


Ptl Pr 1 


A{Q*)-R{Q>) 


3 2 


(pi 1 + S l)(P r r 1 + *r) 


Ph Pr2 


(3.72) 


(Pt2 + a l){p T r2 + 

The first term in (3.70) and (3.71) is the uniform contribution from Q , where the 
second term in (3.70) and (3.71) is the uniform contribution from Q 2 . 


3.1.2 Complex ray field solution on the dark side of the 
caustic 


From (3.66) and (3.67) and using the appropriate form of the incomplete Airy integral 
in Appendix C, the complex ray field on the dark side of the caustic is given by 

E C (P S ) ~ jy/ne- ik <^[An/*Ai('r.)TBn: 1,i Ai'{‘y.)]u(-C a ) (3.73) 

for m(Q o )>0, 

where Ai(7„) is the ordinary or Miller- type Airy function, 7. and a* a are the same as 
in (2.98) and (2.99), and A„ B. are given by 


U1 

A{Q\)-R{Q\) 

1 Pn Pr\ 

1 B. J 

*1 

V (P^ + ^(P" + •?) 
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± A(Ql)-R{Ql) 

o*c 
*2 

Although the values from both complex reflection points appear in (3.73), in the deep 
dark side of the caustic only the contribution from one of the complex reflection points 
is involved since the contribution from its complex conjugate counterpart violates the 
radiation condition and produces an unbounded result. 


P'AP'r 


(p« + a V)(P r r\ + *; c ) 


(3.74) 


3.1.3 Edge diffracted field solution 

The edge diffracted field is written as follows: 


E\P) = Et(P) + E‘,(P ) , 


(3.75) 


where E* u (P) and E*(P) are the non-uniform and curvature correction terms of the 
edge diffracted field, respectively. From (3.66) and (3.67) and using the appropriate 
expression for the incomplete Airy integral in Appendix C, they are given by the 
following expressions: 


Et(P) ~ £'((?.)• D.. (<?.), 


K(p d . + O 




E!(P) ~ + 

l Ta * Sa . 


+ jB d 


if m(Q 0 ) > 0 , 


(3.76) 


+ (3.77) 


K(P) 


~ e 


- jB d 


-J*(«i + *a) J 

V 


f*’ , “ li '\(7.,(.) 

Ta T Sa . 




7a + 


if m(Q Q ) < 0 , 

1 if Ca > (-7a) 1/2 

2 if Ca < ~(-7a) 1/2 

3 otherwise 


Ca 2 . } ’ 


l = < 


1 


if 7 a < 0, or 


i = < 


1 if Ca > 0 

2 if Ca < 0 


> if 7a > 0, 


(3.78) 


(3.79) 


(3.80) 
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— ♦ — * 

where 7 *, ( a , A d , and B d are given by 


C. = 


A d 

B d 


6k'»i ■ (si - Ro) 
u 0 k l,s 


m 


(<?.) 



2 

ji 

m(Q 0 ) 


U o = < 


,1 

e ; < 


2 v^ 

1 

2 ^ 


+ 1 if Qo 3 5" 

-l if g 0 e S, 

(-7.) 1/4 X 

(-7o)- /4 ^ 

ji'J'A, 

7.- 1,4 A 


if 7 a < 0, 
> if 7a > 0. 


(3.81) 

(3.82) 

(3.83) 

(3.84) 

(3.85) 


Notice that the curvature correction term of the edge diffracted field contains infor- 
mation about the reflection points through the vector quantities A d and B d given 
in (3.84) and (3.85) above. Also, Dnu is given in (3.68) and (3.69), and James’s cor- 
rection factors [20] given in (2.117) and (2.118) can still be employed to improve the 
PO diffraction coefficient away from the optical boundaries. 


3.2 Derivation of the UGO/EUTD solution using 
the method of steepest descent 

In this section, the UGO/EUTD solution for the reflected, complex ray, and edge 
diffracted fields is formulated using the method of steepest descent for the asymptotic 
reduction of the PO radiation integral representation for the scattered field formulated 
in the previous section. We begin with the expression for the scattered field in (3.36), 
that is 

#(P)«e-'*“ /(£,«.;*), (3.86) 

where a is given in (3.37)-(3.40) and J(/9,s a ;fc) is a stationary phase integral in the 
.9-plane given by 

oo 

f(0,s.;k) = J (3.87) 
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with /?, s a , and G(s) given in (3.41)-(3.46). Only the case of m(Q 0 ) > 0 will be 
considered here since the results for m(Q Q ) < 0 can be easily obtained by letting 
/ 3 = 0exp(j2ir/3) and s a = s a exp(jir/3) in (3.87). Then, using the same procedure 
as in section 2.2, that is deforming the original path of integration in (3.87) into 
the appropriate steepest descent paths through the endpoint s a and the two real or 
complex saddle points s 12 = ±(— /3) 1/2 as shown in Figures 2.6 and 2.7, and expanding 
the amplitude function G(s) in Taylor series near each pertinent critical point, (3.87) 
may be written as follows: 


where 


I(/3,S a ]k) ~ -jyfik l/3 G(s 1 )W,(k , ' 3 0)u(s l - s a ) 


n 

I(P, s*] k) 


+ 

+ 


js/nk 1/3 G(s 2 ) ’W 1 (k 2/3 0)u(s 7 - 3 a ) 
k~' /3 G(s*)g:(V /3 l3,k' li 3 a ) + 0(k ■“) if 0 < 0, 
1 if 8 a > 

' 2 if s a < s 7 
3 if otherwise, 

27 r k~ 1/s Q(s 7 ) Ai(k 2/3 0)u(—s a ) 


+ k- l/3 G(s a )gl(k^0,k^s a ) + 0(k~') if 0 > 0, 
| 1 if s a > 0 
I 2 if s a < 0, 


(3.88) 

(3.89) 


(3.90) 

(3.91) 


G(s ia ) = 
G(s a ) = 

G(l, a ,T.) = 

4>(1,,,,t.) = 


G(l l j,T.)(-0) 1/i 

G(l a ,T.) 


±2 


\ i)(l ia ,T.) ’ 


i’(^,r.) 

ITk cosfl; jA(Qjj) • R (QA) 

cos 
4>tt(L 

jkA(Q a ) • D nu (Q a ) 


+ (± + _5_), 

1,2 l 7 "*) \^1,2 Ptl t i2/ \^1,2 P rl,r2 } 


s\ 


< (Pi + * 2 ) ’ 


(3.92) 

(3.93) 

(3.94) 

(3.95) 

(3.96) 
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with D n « (Qa) given in (3.68) and (3.69). The first two terms in (3.88) are the uniform 
contributions from the real reflection points in the lit side of the smooth caustic, where 
the third term in (3.88) is the uniform edge diffracted field contribution in the lit side 
of the smooth caustic. The first term in (3.90) is the uniform contribution from the 
complex reflection point in the dark side of the smooth caustic, where the second 
term in (3.90) is the uniform edge diffracted field contribution in the dark side of the 
smooth caustic. 


3.2.1 UGO expression for the reflected field 

The UGO expression for the reflected field is as follows: 

E'(P l ) 



/ PnPrl 


(Pn+Sl){p r r 1 

+•1) 

+ u 7 E'(Qt') • 77 (Q?)^ 

/ PnP T ri 

p -) k *2 

1 {pn + ^ r ) (Pri 

+ *;) 


(3.97) 

where u, = 1 if Q t € S and zero otherwise, u 2 = 1 if <? 2 G 5 and zero otherwise, and 
p\ T are the radii of curvature of the reflected wavefront given in (G.70) or (G.71). 
The UGO dyadic reflection coefficient, % (<5i, 2 ), is given by 

H (Qi,0 — x e 1 + 7Z. h (Q li2 )e )( e n , (3.98) 

where %., h {Q U v) are given in (2.171)-(2.173). The caustic point Q 0 for the 3-D case 
is found from (3.22), (3.29)-(3.31) and (3.35), usually through a computer search 
procedure. In the case of plane wave incidence with the observer in the far-zone of S , 
Q 0 is simply an inflection or zero-curvature point satisfying /c, (Qo) = 0, where Ki(Q 0 ) 
is a principal surface curvature at Q 0 . When the observer is in the deep lit side of the 
smooth caustic, 7 l Mih -» R. itl = and (3.97) becomes the usual GO expression for 
the reflected field. When the observation point lies on the smooth caustic, we have 
that p T tl n and the spread factors in (3.97) become singular. In this case the 

UGO scalar diffraction coefficients, 7t, h , assume the following form: 


N 


(?;. + »:) (<>;. + «:) 
PtoPr 0 


(3.99) 


where m(Qo) is given in (3.26), and thus the singular spread factors are cancelled and 
the reflected field remains finite. 
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3.2.2 UGO expression for the complex ray field 


The UGO expression for the complex ray field in the dark side of the caustic is as 
follows: 




Ptl,l2 PH, r* 


e-^ 2; m(C?o)<0, 


^ {pit,* + *z) (prU + a r 2 ) 

(3.100) 

where p\ c r are the radii of curvature of the complex ray wavefront and are simply ana- 
lytic continuations in complex space of the usual real ray wavefront radii of curvature, 
Pt, T > given in (G.70) or (G.71). The dyadic UGO complex ray reflection coefficient, 
7l c is given by 

ft‘(<?;,J = ft; e';e” + ft; e , ( 3 . 101 ) 

where R,\ h is given in (2. 179)— (2.181). When the observer is in the deep dark side of 
the smooth caustic, ft' „(<?, J -» R\ h = =Fl and (3.100) reduces to the usual complex 
GO expression for the complex ray field. 

3.2.3 EUTD expression for the edge diffracted field 

The EUTD expression for the edge diffracted field is as follows: 

#(P) ~ 5 W.) i 


pi 




a i {pi + *f) 

where T) {Q a) is the EUTD dyadic edge diffraction coefficient given by 

V(Q a ) = V.(Q a )e i J d ± +V ll {Q a )e i ] .e - 


(3.102) 


V.AQ*) = 


-trU 


2 sin ‘d a \/2irk 


F{l[a,C)- 


(3.103) 

(3.104) 


[cos ( 2 ^) cos(*i±^)J 
-ft(7«» Ca) is the EUTD transition function given in (2.185), ip' a and are the angles 
between the projected incidence and observation directions on a plane perpendicular 
to e„ (unit tangent vector to the edge at Q a ) and t a (= e a x n„, n a is the outward unit 
normal of S at Q a ), i? a is the usual angle of diffraction, and p d a is the distance from 
the caustic at Q a and the second caustic of the diffracted ray given by 


1 1 n e - (?'-?•) 

_|_ c V a at 


s' 


R a sin 2 


(3.105) 
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where n c is the associated unit normal vector to the edge at Q a directed away from 
the center of curvature (in general ^ n a ), and R a > 0 is the radius of curvature of the 
edge at Q a . The parameters of the EUTD transition function 7 a and £ a for the 3-D 
case are given by the following expressions: 


la 


U a 

Ca 


U 0 


sini9 0 (cos <p' 0 + cos y> 0 ) 

\LM,<PaMr ’ 

Sg n[Ld{<P' a ,Va\tia)] , 

| sin t?p(cos <p' 0 4- cosypp) - sin t? a (cos + cosy> a )| 1/2 

Uo \L d (<p' a ,<p a Mr 

-|-1 if Q 0 3 S 

< 

-1 if Q 0 € S, 


(3.106) 

(3.107) 

(3.108) 

(3.109) 


where the angles d 0 , <p ' 0 , and <p 0 associated with Q 0 are defined in the same manner as 
i? a , (p' a , and <p a associated with Q a , and L d (<p' a ,<p a ] t? a ) is the EUTD distance parameter 


given by 


Ld(<P' a , Va\$a) = 


m 20 (Q a ) = 


+ 

m u (Q a ) = 

m 02 (Q a ) = 

+ 

m 3 o(Qo) = 

+ 

+ 


+ 


m 2l{Qa) — 


{ m 30 (Q a )ml 2 (Q a ) - m n (Q a ) 


2 m\(Q a ) 

[3m 21 (Q a'j'WoiiQ a.) Tre i2(Qa) Wl ll(Qo)]} ? 

sin 2 8' 

«i(Q a ) sini? a (sin tp' a + sin ip a ) + 


s’ 


sin 2 8 ta 
~ ’ 
a 

cos 0\ a cos 0 f ra cos 0 la cos 0 T 


s' 


8° 


sin 2 0* 

« 2 (<? a ) sin i? a (sin <p’ a + sin <p a ) + 


s' 


sin 2 8 T 


(cOS0 <o + COSf>J a )«2(Q a ) 

sin T? a (sin <p' a + sin <p a ) • V/ti(Q a )] 


cos 

cos 8 la 

^r 

cos 6 

8' 


sin 2 0 * 

3/ti (Q a ) sin t? a sin (p' a H — 


sin 2 0* a 


3«i (Q 0 ) sin a sin y> a + 


1 — 3 cos 2 0\ 
/^(QJsintf.siny?;, + 1 

S' 


(3.110) 


(3.111) 

(3.112) 


(3.113) 


(3.114) 
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+ 


cos 8 r 


1 — 3 cos 2 8, 


_ (n \ cose ia 

m i2 \Qa) — 


K i{Qa) sini? 0 sin <p a + 


/r\ \ - a • » 1 — 3 cos 2 8' 

fh(Q a ) sm i? a sin ip a H — 


+ 


cos 


«2(Qo)sint? a sin <p a + 


1 — 3 cos 2 8 t 


m c (Q a ) = S in 2 t? a ^i + ij, 


(3.115) 


(3.116) 

(3.117) 


where 8' tara are the angles between the incident ray direction at Q a and T a , f a (principal 
surface directions at Q a ), respectively, 8 lara are the angles between the diffracted 
ray direction at Q a and i a , r a , respectively, and n l<7 (Q a ) are the principal surface 
curvatures at Q a in the direction of l a , f a , respectively. It is also understood that 
when L d ((p' a ,(p a ]‘8 a ) < 0, the EUTD transition function in (3.104) should be replaced 
by its complex conjugate version. 

In the derivation of the results for the 2-D case in Section 2.2.3, it was shown that 
when the observation point is on the deep lit side of the smooth caustic and near an 
optical boundary, the EUTD transition function assumes the following form: 


* *V), 

where F(x ) is the UTD transition function and i j is given by 

Ja + C _ Sil1 1 ?a(cOS <P' a + COS lfi a ) 


(3.118) 


(3.119) 


2 (— 7 a ) 1/4 2 (— 7 a ) 1 / 4 | X d (^, y ? a ; t ? a )| 1/3 

From (3.106) and (3.110)— (3.117), near an optical boundary on the deep lit side of 
the smooth caustic, L d (<p' a ,<p a ; t? a ) assumes the following form: 


L d (<p' a ,<Pa]#a) = 


\/k 


Trt2o(<? a )Tn 02 (Q a ) - m 2 ,((g a ) ' 8/2 

2m e (Q a ) 


Near the RSB, <p a ss 7r — <p' a and using (G.69) we have 

wijo (Q a )rn 0 ,{Qa) - m* n (Q a ) w sin 2 i? a cos 2 ( — - — 'j 

Also, near the RSB m e (Q a ) is given by 

m s (Q a ) fa sin 2 ^ ^ , 


(3.120) 


(3.121) 


(3.122) 
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where from (3.105) 


(3.123) 


J_ _ J_ _ 2 (n a -n e )(?i, -w 0 ) 
Pi si R a sin 2 tf a 

Thus, using (3. 1 19)— (3. 122) , rj 2 near the RSB is given by 


r? 2 w 2A: cos 2 


(Va + VA 

' + OPi*p r Ta Sln * *«" 

\ 2 ) 

. Pl(p r la + S d a)(p T ra + a i) . 


= *X'a(v>i + »».)• (3-124) 


(3.125) 


(3.126) 


Near the ISB, <p a ~ 7 r -f ^ and using (G.38)— (G.40) the quantity m 20 m 02 m 2 , at Q 
assumes the following form: 

m 20 (Q a )m 02 (Q a ) - m'XQ*) ~ sin 2 cos 2 g " ) + 3”) ' 

Also, near the ISB p d a = sj, from (3.105) and m e (Q a ) is given by 

m c (Q a ) ~ sin 2 i? a ^ + 7 ) ' 

Thus, using (3.120), (3.125), and (3.126), rf near the ISB is given by 

where L i r are the UTD distance parameters for spherical wave incidence. Both L' r 
and o(<p' a ± tp a ) are defined in [4]. The results in (3.118), (3.124), and (3.127) demon- 
strate that near the optical boundaries in the deep lit side of the smooth caustic, both 
the EUTD transition function and its argument reduce to the conventional UTD for- 
mulations as expected. 


(3.127) 
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Chapter 4 

Plane Wave Scattering from 
Polynomial Surfaces 


In this chapter, a UGO/EUTD solution for the plane wave scattering and diffraction 
from polynomial surfaces is developed. The derivation of the solution follows a similar 
procedure as the uniform asymptotic analysis in Chapters 2 and 3, with a generalized 
form of Chester’s expansion being employed for the asymptotic reduction of the PO 
radiation integral. The observer is assumed to be in the Fraunhofer region of the 
scatterer, and due to the symmetry introduced by this assumption the alternative 
method of steepest descent employed in Chapters 2 and 3 yields the same expres- 
sions, except for the edge diffracted field solution where the resulting edge diffraction 
coefficients are slightly different. 

The predominant scattering contribution from polynomial type surfaces would 
come from the neighborhood of the zero-curvature points on the boundary [22] which 
are associated with caustics of the GO reflected rays. In addition, the contribu- 
tions from the complex stationary phase points in the PO integral are interpreted 
as diffraction from the zero-curvature points rather than complex ray contributions. 
This zero-curvature diffraction interpretation is indeed equivalent to the complex ray 
interpretation, yet it eliminates the need for a complex extension of the scattering 
boundary and the evaluation of geometric parameters in complex space. 

The PO approximation is also employed in this chapter, thus, only the first order 
scattering mechanisms will be considered, namely reflection, zero- curvature diffrac- 
tion, and first order edge diffraction. The canonical function involved in the expression 
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for the total scattered field from a polynomial surface is the generalized incomplete 
Airy integral [22]. In the special case of a cubic polynomial surface, this generalized 
integral reduces to the usual incomplete Airy integral and the expressions for the scat- 
tering mechanisms are in agreement with the results of Chapters 2 and 3 specialized 
for the case of plane wave incidence with observation in the far-zone. In the case of 
a fourth order polynomial surface, the generalized incomplete Airy integral reduces 
to a special case of the incomplete Pearcy integral. It should be understood that the 
results of the analysis employed in this chapter do not remain uniformly valid for 
observation points in the intermediate or near zone of the scatterer. In the near zone 
problem, the dominant contribution to the PO integral is no longer localized at the 
zero- curvature points and a more general treatment along the lines of Chapters 2 and 
3 must be employed. 

The outline of this chapter is as follows: The geometric properties of polynomial 
surfaces and their relation to the scattering mechanisms are discussed in Section 4.1. 
The UGO/EUTD expressions for the reflected, zero- curvature and edge diffracted field 
contributions from finite two and three dimensional polynomial surfaces are derived 
in Sections 4.2 and 4.3, respectively. 

4.1 Polynomial surfaces and high-frequency scat- 
tering considerations 

Let’s consider a two dimensional polynomial surface. As a plane parametric curve it 
may be expressed as follows: 

f(u) = ui + y(u)j/, a<u<b. (4-1) 

where Y{u ) is an nth order polynomial function given by 

y(u) = P n (u) = a 0 + a,u + a 2 u* + f a n u n . (4.2) 

The distinguishing feature of polynomial defined surfaces are the existence of zero- 
curvature points which are real roots of the second derivative of the surface, that 
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IS 


F"( u m) = a 2 + a,ti m + a 4 t£ + |-a n <, 2 = 0; m = 1,2 • • • ,n - 2 , (4.3) 

a n = n(n - 1 )a„ . (4.4) 

If a < u m < b is a first order real root of (4.3), then the surface near u m is locally 
cubic, that is Y’"'(u rn ) ^ 0 and y(ii) may be approximated as follows: 

Y(u) « Y(u m ) + Y'(u m )(u - u m ) + ^p^(« - u m y . (4.5) 

In general, if u m is an 1th order real root of (4.3) or an /th order zero- curvature point 
with 1 < l < n — 2, we have that 


F ( *>K) = 0; 2<*<Z+1 (4.6) 

y<‘ <’>(<.„) / 0, (4.7) 


and near u m , F(u) may be approximated as follows: 

y ( u ) » r(«„) + r(»„)(ti - u.) + - ».)'« . (4.8) 

It should also be noted that near an odd order zero- curvature point (/ = 1,3,5,* • •) 
the polynomial surface is concave-convex, whereas near an even order zero- curvature 
point (/ = 2,4,6, • • •) the polynomial surface is strictly concave or convex. The unit 
tangent, normal, and surface curvature vectors for the polynomial boundary defined 
in (4.1) are as follows: 


t(u) 

n(u) 

/c(u) 

k(v,) 

h(u) 


r'(tt) i + y'(u)y 

k(u) h(u) ’ 

f(u) x z , 


f "(u)h(u) — r'(u)h'(u) 


|«(u)| = 


h s (u) 

Y"{u) 


= — /c(u)n(u) , 


h>(u) ’ 

|f*'(u)| = |l 4- [y'(u)] 2 | ' = Jacobian. 


(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 


The results of Chapter 2 and specifically Equation (2.14) show that when both 
the source and the observer are in the far-zone of the scattering surface, the zero- 
curvature points are associated with caustics of the GO reflected rays. Figure 4.1 
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Figure 4.1: First order scattering mechanisms associated with a polynomial boundary 
illuminated by a plane wave. 

illustrates a polynomial boundary illuminated by a plane wave. The boundary C 
contains an edge point, Q e , and two zero-curvature points, Q p 1>p2 , with their associated 
GO ray caustics. When the two zero-curvature points are sufficiently far apart, the 
high-frequency scattering from C can be effectively analyzed by treating each caustic 
separately. 

The scattering mechanisms of interest in the analysis that follows are reflection 
that remains uniform as the reflection point Q r — > Q p , first order edge diffraction 
that remains uniform as Q e —* Q p , and zero-curvature diffraction or complex ray 
(evanescent) field contribution that compensates the discontinuity of the reflected 
field across the GO ray caustics. Although additional higher order mechanisms such 
as edge exited surface rays and whispering gallery modes may also exist and con- 
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tribute significantly to the scattered field in certain regions of the coordinate space 
surrounding a polynomial surface, they will not be considered in the present analysis. 


4.2 Uniform asymptotic analysis (2— D case) 


The canonical geometries for the uniform asymptotic analysis are shown in Figure 4.2. 
Figure 4.2a shows a semi-infinite perfectly conducting concave- convex polynomial 
surface that contains an odd order zero-curvature point. Figure 4.2b shows a semi- 
infinite perfectly conducting polynomial surface that contains an even order zero- 
curvature point. The electric field E 3 scattered by the 2-D polynomial boundary C 
illuminated by a plane wave and observed at a point P in the far-zone of the boundary 
can then be expressed in terms of the usual radiation integral over the electric current 
J on C as follows: 

E'{P) ~ ^ / c [ 5 x J x AQ')] H^ 2 \kR)dl ' , (4.14) 

where 


R = | Q'P j « p — t(Q') • 3", (far-zone approximation) 
s — x cos 6 + y sin 8 , 
r(Q') = x x + Y(x)y = r(x) ; a < x < oo , 


(4.15) 

(4.16) 

(4.17) 


y(*) is a polynomial function, p — OP , and s is the observation direction. As 
before, the PO approximation to (4.14) is employed using the GO currents in (2.2), 
where the incident magnetic field at Q' under the plane wave illumination assumption 
is given by 


H\Q') = Z~ l (s' x p) e-W** ; S'-p = 0 , (4.18) 

= Z; l rxE\Q'), (4.19) 

E\Q‘) = , (4.20) 


p is the polarization vector of the incident plane wave, and s' is the direction of 
incidence given by 

s' = — £ cos O' — y sin 6' . (4.21) 
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Figure 4.2: Canonical geometries for the uniform asymptotic analysis: (a) polynomial 
boundary containing an odd order zero-curvature point, (b) polynomial boundary 
containing an even order zero- curvature point. 
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Thus, using the PO approximation and the large argument form of the Hankel func- 
tion given in (2.6) the far-zone scattered electric field is given by 


E-{6',e ( 4 . 22 ) 

VP 

where ilssp was used in the magnitude term of (2.6). The vector quantity e(9',6) is 
the angular field dependence and is given by the following stationary phase integral: 

oo 

e(6\6) = J F(l') d£' , (4.23) 

t a 

where 

H?) = { i x a x [n(<') x S ' x p] } , (4.24) 

m = f(2') ■(»-?•'), (4.25) 

and n(/') is the outward unit normal to the surface at Q‘. The next step in the 
procedure is the uniform asymptotic evaluation of e(6',6) using a generalized form 
of Chester’s expansion [22]. 


4.2.1 Uniform asymptotic evaluation of e(0',0) 

First, the stationary phase integral in (4.23) must be transformed into a canonical 
form. This is done by using the surface approximation near the zero-curvature point 
given in (4.8), or alternatively by expanding the phase function in Taylor series around 
the zero- curvature point. Both methods give the same result, and if l is the order of 
the zero-curvature point, that is x p is an Ith order root of the equation f''(x) = 0, 
the appropriate transformation is as follows: 


where v = / + 2, 


4 >(?) = t ( s ) = a + /3s - f , 


a = 'r(O) = , 

0 = r --( s - r )^ 


(4.26) 


(4.27) 

(4.28) 
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and t p is the unit tangent vector at Q p 

Hence, using (4.26) the stationary phase 


integral in (4.23) transforms to the following canonical form: 



— 

"oo” 

e(0',O) = e J *“ J G(s) ds , 

•o 


(4.29) 


where 





\? (~ ~<\ 61 
s • = [*•'{•-•) Ts 

-t p - (s-s l ) 

• = J 0 

■ti " 

o 

II 

m 

% -3 

(4.30) 


+1 if Q p 3 C 

u a = < 

-1 if Q, e C, 



(4.31) 


<?(») = ?(.*)%, 



(4.32) 


dt _ dt dx _ Uu- 1)!' 
ds dx ds m(x, v) 

L 

V 


(4.33) 


, \ r { M (z p ) ■ (5 - s ') 

m{x,v) = KP J \ ', 

l r (®)l 



(4.34) 

— 

u » (dt\ 

oo = oo exp -jarg 1 — 1 

5 


(4.35) 


and t a is the unit tangent vector at Q a . The proper branches for /3, s a and dt' jds 
depend on the sign of m(x p , v). 

Next, we employ a generalized form of Chester’s expansion [22] for the amplitude 
function in (4.29), that is 

G(s) = £ [o m>0 (a*- 1 + P) m + a mil s (s"- 1 + /3) m 

m—0 

+ — + S m „-, s*- 3 + P) m ] , (4.36) 

and since only the leading terms in the asymptotic expansion of (4.29) will be re- 
tained, (4.36) may be written as follows: 

0(>) = £ •"<- + (»"“ + P) 9 (‘) , ( 4 . 37 ) 

n=0 

where 

»W = £ (-— + ^) m - 1 + a„., < ‘ + 0 )-' 

m= 1 

+ ••• + 2 s v ~ 3 (a* -1 + /3) m-1 ] . (4.38) 
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Due to the symmetry introduced by the plane wave far-zone observation assumption, 
only the vector coefficient a 0 , 0 is non-zero in (4.37), that is 

O 0> , o = — * nd ( 4 - 39 ) 

Oo.n = 0; n = l, — 2, (4.40) 


where d?(s),) is given by 


G(s k ) = F(£„) 


dt 


dt 


'(j, _ i)(_ y g)^e ;2( *- 1)( ‘'“ 2),r/(,/ " l) 

ds 

,= i * 

W) 


(4.41) 

(4.42) 


and s h are the stationary phase points of r(s), that is t'(s*) — 0, and are given by 
s„ = (-/S) z±r . k = 1,2,- ••,«/- 1. (4.43) 

Thus, using (4.37) along with integration by parts e(0',Q) may be written as follows: 


e(6',8 ) ~ e 




a 0 , 0 J - j^g(3 a )e ik ^ + ’^ 

* »a 

and using (4.26) and (4.37) we have that 

“oo” 

e(M) - e iha d 0fi j e ik{ ^‘ v ' v) ds 

UqjO ^(^a) 


+ 0(*T 2 ), (4.44) 


+ e 




+ 0(fc -2 ) 


We then let s = tk~ u \ 7 = /3k’- l '- 1),l/ , and C„ = s a k 1/v in (4.45), thus 

oi ka 


(4.45) 


?(«■.») - S 


+ e 


Jkr{» a ) 


CM 


L >s^cr + c-) (4 - 46) 

Finally, depending on the sign of m(jc p ,i/) the angular field dependence e(0,8) is 
given by 


e jka ^ -=r 


e(M) ~ ^e-' 0 ^( 7 ,Ca^) + ^“ ) 


' jy/k ( 7 + Cr 1 ) 


+ 


(4.47) 
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if m(x s ,u) > 0, 

jFV.) 

H'{la) jVk(~f+Cr l ) 

if m(x p ,i/) < 0, 


g>*« ____ 


(4.48) 


where the 0(& J ) terms were omitted, and the parameters 7, £ a , and e 0 are given by 
the following expressions: 


7 

u P 

C 

w(x, 1/) 
e 0 


u p t p - (s- s’) u;(a:p, 1/) , 

sgn[m(* p ,i/)], 

£,* (a - S’ 1 ) ty(a,i/) - t p ■ (s - s’) u/(;Cp,i/)|‘'~ 1 , 

(1/ - 1)! " 
m(x,i/) ’ 

(- 


g y f„ ' IJEL 
p-i h (<>) vW.)i 


(4.49) 

(4.50) 

(4.51) 

(4.52) 

(4.53) 


The function ^^(7,^,1/) is the generalized incomplete Airy integral [22] defined as 
follows: 

00 

K u ( 7, C, I/) = f e*-'**’” dz ; 1/ = 3, 4, 5, • • • (4.54) 

c 

The second term in both (4.47) and (4.48) is a purely edge diffraction contribu- 
tion where the first term in (4.47) and (4.48) involving the generalized incomplete 
Airy integral contains all three scattering mechanisms, namely reflection, edge, and 
zero-curvature diffraction. These scattering contributions may be separated using the 
properties and complex plane topology of ’K u (‘j,£,u) in Appendix C. Also, the pa- 
rameters a, 7, £ a , and e 0 can be expressed in terms of the local geometric parameters 
on the boundary relevant to each scattering mechanism. 


4.2.2 Reflected field contribution 

There exist two cases. The first case corresponds to a concave- convex polynomial 
boundary that contains an odd order zero-curvature point. In this case v is odd 
in (4.54) and the reflected field consists of two specular contributions that only exist 
in the lit side of the caustic. The second case corresponds to a strictly concave or 
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Figure 4.3: Reflection near an odd order zero-curvature point. 

convex boundary that contains an even order zero-curvature point. In this case v is 
even in (4.54) and the reflected field consists of a single specular contribution that 
exists on both sides of the caustic. The concave-convex case will be considered first. 


A. Reflection from a concave-convex polynomial boundary 

The geometry for the reflected field from a polynomial boundary that contains an 
odd order zero- curvature point is shown in Figure 4.3. Using (4.22), (4.47), (4.48), 
and the appropriate expression for K u (^,C,v) for v odd and 7 < 0 in Appendix C, 
the reflected field in the lit side of the caustic is given by 


E r (d',Q\p) 


{«x e >kQ1 ei J,(- 7 i,«') + e *" 1 


if 

m(x p ,i/) < 0, 

e~ ikp 

E r {6'Ap) 


{u, e, / r *(- 7 . ,») + J r (- 7*, ")} — 


if 

m(x p ,u) > 0 , 


(4.55) 


(4.56) 
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where tt li2 = 1 if Q 12 £ C and zero otherwise, and the function 4 .( 7 , v) is defined as 
follows: 

* foe exp( jir/2v) 

4 ( 7 , 1 /)=/ eHv+^M dz ; V is odd. (4.57) 

J 00 exp(— j3?r/2i/) 

The sign of m(x p ,v) designates which side of the boundary C is illuminated. For 
the situation depicted in Figure 4.3 m(x p , v) < 0 . From (4.24)-(4.26) and (4.53) the 
parameters a li2 , 7 1i2 , and e , i2 may be written as follows: 

<* 1,2 = ^(<?i, 2 ) • (S” - - ?') ± u p k~ l 7iT‘ > (4.58) 

7,2 = {k(^)\[r(Q l ,)-f(Q p )}.(s-r)\]~ :r , (4.59) 

(4.60) 


r - c >f / x/ ~ 1 7 ^t 

- 1 ' 2 ” 6 V 27T 7l ' 2 


p- R (^ 1 , 2)1 \/ 2 I^(Ql, 2 )| cos 8 1 2 , 


where 4 (£?t, 2 ) is the GO dyadic reflection coefficient given in (2.60) and (2.61), and 
p a (C? 12 ), 0' 12 are the surface radius of curvature and angle of incidence at Q 12 , respec- 
tively. Thus, using (4.58)-(4.60) the reflected field from a concave- convex polynomial 
boundary for the case of m(x p ,i /) < 0 , that is p g (Q t ) < 0 and p„(Q?) > 0 , may be 
written as follows: 


E'iO'Ap) 


^E'iQtyn'mJ ^PM) cos 6{ e jkl,< - Ql) ‘ T 


+ «, E\Q,y n (Q,)J^P„(Qi) cos 61 




-jr (4 - 61) 


where E‘(Q l7 ) is the incident plane wave at Q 12 given in (4.20), and ft (Q l:7 ) is the 
UGO dyadic reflection coefficient and is given by 


H {Q ,2) — e ± e ± + 'R'hi.Qij ) e nC 1 1 


(4.62) 


K'AQi,*) = (4.63) 


R., h = =Fl , for a perfectly conducting boundary. 


(4.64) 


Equation (4.61) shows that the complex conjugate version of the UGO reflection 
coefficient should be used when p a {Qi, 7 ) < 0 . When u = 3 , (4.63) reduces to the 
UGO reflection coefficients derived in Chapter 2 and given in (2.171). 


85 



/ 




Figure 4.4: Reflection near an even order zero-curvature point. 

B. Reflection from strictly concave or convex polynomial boundary 
The geometry for the reflected field from a polynomial boundary that contains an 
even order zero-curvature point is shown in Figure 4.4. Using (4.22), (4.47), (4.48), 
and the appropriate expression for 7 , v) for v even in Appendix C, the reflected 
field on either side of the caustic is given by 

e~* kp 

E r (9',e-,p) ~ u r e J * ap e r /;(- 7 r ,z/)— — ; ifm(* p ,j/) < 0 , (4.65) 

VP 

e~* kp 

E T (6',6\p) ~ u T e jkar e r 7 r (— 7 r , u) — ; if m(x p , i/) > 0, (4.66) 

VP 

where u T = 1 if Q r € C and zero otherwise, J r ( 7 ,*/) is the same as in (4.57) with v 
even , and the parameters a P , 7 P , and e P are given by 

«. = ■ r (Q,).(5--}')±u,r l CVr)' l ' A ’ ( 4 - 67 ) 

7 . = {* ( r ^ T )rW.)-r (<?.)] •(*-?■)!}” . (4.68) 


I 
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(4.69) 


. _ j V — 1 u-1 

* _ 

I 

~ e V 2 t 

P ■ R (Qr) 

V 


^k(<?r)|cos 0*. 


Thus, using (4.67)-(4.69) the reflected field from a strictly concave or convex polyno- 
mial boundary for the case of m(x p1 v) < 0 or p g (Q r ) > 0 depicted in Figure 4.4 may 
be written as follows: 


_ / 1 e~ jkp 

E'(e',o ]P ) ~ u, p(Q,y n -j=- . 


(4.70) 


where E’(Q r ) is the incident plane wave at Q r given in (4.20) and H ( Q r ) is the UGO 
dyadic reflection coefficient and is the same as in (4.62)-(4.64) with v even. Again, 
when p„(Qr) < 0 the complex conjugate version of H (Q r ) should be used. 


4.2.3 Zero-curvature diffracted field contribution 

As in the case of the reflected field, we also have two cases for the zero-curvature 
diffracted field. The first case involves diffraction from a an odd order zero-curvature 
point associated with with a concave-convex polynomial boundary. Since the reflected 
field is discontinuous across the caustic in this case, the zero curvature diffracted field 
has to compensate this discontinuity. The second case involves diffraction from an 
even order zero-curvature point associated with a strictly concave or convex polyno- 
mial boundary. In the latter case the reflected field is continuous across the caustic, 
thus the zero-curvature diffracted field only accounts for higher order surface discon- 
tinuities. 


A. Diffraction from an odd order zero-curvature point 

The geometry for the diffracted field from an odd order zero- curvature point is 
shown in Figure 4.5. Using (4.22), (4.47), (4.48), and the appropriate expression 
for f° r v °dd and £ < 0 in Appendix C, the odd order zero-curvature 

diffracted field is given by 

_ e~ ikp 

E p (6\6-,p) ~ u 0 e ,hap 2 p I pt (-i„v)— — , (lit side) (4.71) 

VP 

-» e~ jkp 

E p (6', 0; p) ~ u 0 e ikap Z p I pd { 7 P , v) —— , (dark side) (4.72) 

VP 
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where u 0 = 1 if Q p 6 C and zero otherwise, and the functions 7 p/ (7,i/), I pd ( 7,1/) are 
defined as follows: 


KM = / 

J a 

KM = / 

J c* 


ooexp(-j 3 ir/ 2 »/) 

e ,(l,+ ‘ /l/) <fz ; v is odd, 

00 ex P[j( -*’ + 3 * 721 ')} 

00 exp( jir/ 2 »/) 

e )(y*+ *''/<') fc. v j s odd. 

00 exp[j(ir-ir/ 2 i/)] 


(4.73) 

(4.74) 


From (4.24), (4.25), (4.27), (4.28), and (4.53) along with some algebraic manipula- 
tions, the parameters a p , 7 P , and e p are given by 


where 


& P 


7 P 


r{Q P ) -(? p - «*)> 

k (v ~ l),v | cos <p p + cos <p p \ w(Q p , i /) , 


S K 1 v> ~ 2 

e 3 < 
y/iir 


L v (<P' p ,¥’p;i / ) [p- (~e±e± + ej , e ;)] , 


™{Q p,v) 

m(Q p ,u) 


L pty f ,'Pp\ v ) 


(i/ - 1)1 I- 

m(Q p ,v)\ 

— k {v - 2) (Q P ) n p -(s — s’) , and 
™(Q P 1 U ) COS - — 'j . 


(4.75) 

(4.76) 

(4.77) 


(4.78) 

(4.79) 

(4.80) 


The angles <p p and <p p are defined in Figure 4.5, V , '~ 2) (Q P ) is the u — 2th derivative 
of the surface curvature at Q p , and n p is the outward unit normal at Q p . Thus, us- 
ing (4.75)-(4.77) the zero-curvature diffracted field from a concave- convex polynomial 
boundary may be written as follows: 


E’id’Ap) ~ u 0 E\Q p y V co (Q P )e 3k ^ >P - , (4.81) 

VP 

where E^Qj,) is the incident plane wave at Q p given in (4.20), and V co (Q P ) is the 
EUTD dyadic odd order zero- curvature diffraction coefficient given by 


T>'°(Q P ) 


vzWpi'Ppi 1 ') 


v:°(<p' p , <P P ; v) Wx + T>l 0 (ip' p , <p p ; u) e *e 

hii-lp,”), 


=F — 70=- tpW.W') 




Kd{lpt * 


II ’ 

(lit side) 
(dark side) 


(4.82) 

(4.83) 


Notice that when v — 3, I p i(— 7 P ,3) = 0 and the EUTD expression for the zero- 
curvature diffracted field is equivalent to the complex ray field solution derived in 
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c 


in 

//? 

Caustic 

Figure 4.6: Diffraction from an even order zero-curvature point. 

Chapter 2 specialized for the case of plane wave incidence with the observer in the 
far-zone. 

B. Diffraction from, an even order zero- curvature point 

The geometry for the diffracted field from an even order zero-curvature point is shown 
in Figure 4.6. Using Equations (4.22), (4.47), (4.48), and the appropriate expression 
for 5^(7, Cj *0 f° r v even and ( < 0 in Appendix C, the even order zero- curvature 
diffracted field on either side of the caustic is given by 

E’(d',9-,p) ~ u 0 e^r p j;( 7p ,i/)-— ; ifm(* p ,»/)<0, (4.84) 

VP 

_ e~ ikp 

E p (9',6-,p) ~ u 0 e J * Qp e p J p (7 P ,i/) — — ; if m(* p ,i/) > 0. (4.85) 

VP 

Again, the sign of m(x p , 1/) depends on which side of the polynomial boundary is 
illuminated, with m(x p , 1/) < 0 corresponding to the situation depicted in Figure 4.6. 
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The parameters u 0 , a p , 7 P , and e p are the same as in the odd order zero-curvature 
diffraction case. The function J p ( 7, v) is defined as follows: 

A [oo*xp{j-n f2v) 

J p {j,v) ~ I dz \ v is even. (4.86) 

J oo exp[j( ir — Sw/2i/)] 

Thus, using (4.75)-(4.77) the even order zero- curvature diffracted field may be written 
as follows: 

__ e~ 3kp 

E p (6\ 6 ; p) ~ u 0 £’(C? P )- V“(Q P ) e’""** —=r , (4.87) 

VP 

where E'(Q P ) is the incident plane wave at Q p given in (4.20), and T> CC {Q P ) is the 
EUTD dyadic even order zero- curvature diffraction coefficient given by 


V“(Q P ) 




T> c ;{<p' pi <p P ; v) e ± e ± + v~{tp' p , <p p ; v) e\e\ , 

J p (1p^) Xrn{x p , 
J P {l P , ^) Hrn(x p , 


e jf 


») < o, 

v) > 0. 


(4.88) 

(4.89) 


4.2.4 Edge diffracted field contribution 

The geometry for the edge diffracted field from a polynomial boundary is shown in 
Figures 4.7a and 4.7b. Unlike the reflected and zero-curvature diffracted fields, the 
two methods of analysis employed in Chapter 2, namely Chester’s expansion and the 
method of steepest descent, result in two different expressions for the edge diffracted 
field. However, the resulting diffraction coefficients make use of the same transition 
function and are in agreement to the first order term for observation points near the 
RSB. They also reduce to the exact expression for the half- plane diffraction coefficient 
outside the transition regions. When u = 3, the general transition function reduces 
to the EUTD transition function derived in Chapter 2. The results obtained using 
Chester’s expansion method will be considered first. 


A. Chester’s expansion formulation for the edge diffracted field 
Using (4.22), (4.47), (4.48), and the appropriate expression for K u ( 7,C, l/ ) in Ap- 
pendix C, the edge diffracted field may be written as follows: 

EfV'Ap) ~ {e jkaa e a / a ( 7 a, Ca) ^) 
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. ■,* Lit side 


Caustic 

U <S> • 


(a) 



Figure 4.7: Diffraction from an edge in a polynomial boundary: (a) concave- convex 
case, (b) strictly concave or convex case. 



+ e 


J M*a) 


W) 

k<t>'{ta) 1 i(7 + Cr’) 


+ -7 


EWAP) ~ {e^e 0 / a *( 7a ,C,«/) 

'jF(4) e a 


_l_ 


,-i*p 




; if m(x p ,i/) > 0, (4.90) 


-jfcp 


— — ; if m(x p , i/) < 0. (4.91) 

V? 


*^'(4) i (7 + Cr‘)_ 

From (4.24)-(4.26) and (4.49)-(4.53) along with some algebraic manipulations, the 
parameters a a , 7 *, £ a , and e a are given by 

<*a = f(Q a ) • (s d - s’) - u p k~ l (7„Ca + O) , (4.92) 

7 a = Up A ( ‘ /_ 1 )/ *'(cosyj' + cosy> p )u;(x p ,i/) , (4.93) 


c. 


u, 


k ilu 


(cos ip' a + COS <p a )w(a,l/) - (cOS + COS^)™^,!/)]' 


U ( 7 4 . £*'-• )-lil±L 
pl7a + C ‘ > k<i>‘(£ a ) 


where u p and w(x,v) are given in (4.50) and (4.52), respectively, and 

f{Q = ]j H [?• (-exe x + ejejj)] cos ) • 

The function / a (7, £, ^) is defined as follows: 

A Pooexp[>Vo(*')I 

/«( 7,C,*') = J ( dz ; C > 0, 

^ a ( 7 ) C ) if 1 / is odd 


^1(7? Cj ^) = S 

V’oH = 


if C < 0, 


-/„( — y, -£) if t/ is even, 

tt / 2 u if [sgn(C)] 1 ' -1 (7 + O' -1 ) > 0 

—Ztt/ 2 v otherwise. 

Also, from (4.24) and (4.25) we have that 

HQ = [p- (-sin^ej.ej. + smp.ejejf)] , 

${t a) = cos <p' a + cos <p a . 


(4.94) 

(4.95) 

(4.96) 

(4.97) 

(4.98) 

(4.99) 


(4.100) 

(4.101) 


Thus, using (4.92), (4.95), (4.96), (4.97)-(4.101) and the PO correction factors 
in (2.117) and (2.118), the edge diffracted field may be written as follows: 


£?(M;P) - £’(<?.)• 


d e 


-jkp 


Vp ’ 


(4.102) 
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(4.103) 

(4.104) 


V;(Q a) 

i* 

Cd 


v: 1 (<p'a,<pa)e ± e ± + Tr l (<p' a ,<p a )e , „e', 

-e~ J ' f 1 ^(7*1 

2 \/ 27 rA: cos ^ “ ~ yg ) cos 

{ 7 a if i/ is odd 

sgn(Ca)7» ^ 17 is even 7 

ICJ. 


The EUTD transition function J-( 7 , £, v) is defined as follows: 

^(7,C,^)=i(7 + C- 1 ) /*(7,C,^)e^ c+c ‘' /, ' ) , 


(4.105) 

(4.106) 


(4.107) 


and when v = 3 , it reduces to the expression derived in Chapter 2 and given 
in (2.185). It is also understood that when p g {Q a ) < 0 , the complex conjugate version 
of should be used in (4.104). 

When the observation direction is near an optical boundary and far removed from 
the caustic, that is 7 d + » 0 with 7 d < -1, I a {~1,C,v) can be expressed in terms 

of the Fresnel integral and the transition function reduces to the UTD 

transition function F(x) as was shown in Chapter 2 for the v = 3 case. Therefore, 
near the optical boundaries and away from composite shadow boundaries the EUTD 
diffraction coefficient in (4.104) is in agreement with the regular UTD diffraction 
coefficient. In addition, away from transition regions, that is 7 d + ^ U it can 

be easily shown by employing integration by parts in (4.97) that 7 , £,i/) ~ ► 1 
and (4.104) reduces to the exact expression for the half-plane diffraction coefficient. 


B. Steepest descent formulation for the edge diffracted field 

For the steepest descent analysis we begin with the expression for the scattered field 
in ( 4 . 22 ) and (4.23), that is 

- _ g-jfcp 

E‘{6 ', 6\ p) « e ika /,(/?, a 0 , v\ k) — — , (4.108) 

VP 

where a, (3 , and s a are given in (4.27), (4.28), and (4.30), respectively, and m(x p , v) 

— ♦ 

given in ( 4 . 34 ) is taken as positive in the analysis that follows. I,(/3,s a ,i/]k) is a 
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stationary phase integral in the complex s-plane and is given by 

oo 

/.(&*., v\k) = J G(s) ds . 


(4.109) 


&) has the same complex plane topology as the generalized incomplete Airy 
integral in Appendix C. Thus, by deforming the original path of integration in (4.109) 
into the appropriate steepest descent paths through the integration endpoint as shown 
in Figures C.1-C.3, the edge diffracted field contribution in (4.108) may be written 
as follows: 


EM\e;p)*e^I d ((3,s a ,i,;k) 


o-jkp 


VP ’ 


where 


-♦ , fooexpfji P d (u)] _ 

W, s m , v,k)= j G{s ) e J * (p * + * M ds , 

J Bn 


(4.110) 


(4.111) 




7t/2i/ if /3 + < -1 > 0 
— 37 t/2i/ otherwise 

7T — 7t/2i/ if /? + > 0 

— 7r + 37 t/2i/ otherwise 

— tt + tt/2i/ if /3 + s" -1 < 0 
7r — 37 t/2j/ otherwise 


; for > 0, (4.112) 

* ; for s a < 0 and v odd, (4.113) 
^ ; for s a < 0 and v even. (4.114) 


Next, it is assumed that G(s) is regular and slowly varying along the path of integra- 
tion in (4.111) and thus, for large k the principal contribution to the integral would 
come from the neighborhood of the integration endpoint. Upon expanding G(s) in 
Taylor series around s a and making the substitution s = the edge diffracted 

field may be written as follows: 




EWA p) 




Jka a 


F(t.) 


me.) 

if m(x p ,i/)>0, 

?(t.) 


< 7 . +C - ') 


-jkp 


Vp 


~ — e 


jka a 


W'(t a ) 

if m(x p ,i/)<0, 


( 7 . + C“)/;( 7 .,C.,*') 




sTp 


(4.115) 


(4.116) 
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where a„, 7 *, ( a are the same as in (4.92)-(4.94) and /«(7,C» I/ ) i g defined in (4.97)- 
(4.99). Finally, using (4.100), (4.101) and the PO correction factors in (2.117) 
and (2.118) the steepest descent formulation for the edge diffracted field for the case 
of p„{Q a ) > 0 is given by 


EWAp) ~ 
VHQa) = 
£>:»*) = 




d e 


-jkf> 


E'(Qa)' VliQa) e ^ , 

V'.' 2 {<p' a ,<p a )e 1 e 1 + T>*{<p' a ,<p a )e\fil, 

— e~ j < [ 1 1 

2\/27r k cos 


(4^l) cos(*^)J 


Cdt 5 


(4.117) 

(4.118) 

(4.119) 


where j d , ( d , and v) are the same as in (4.105)-(4.107). When p g (Q a ) < 0 the 

complex conjugate of ^( 7 ,C> I/ ) should be used in (4.119). 

The major difference between the steepest descent diffraction coefficient in (4.119) 
and the diffraction coefficient in (4.104) derived using Chester’s expansion method is 
that in (4.119) the transition function operates on the entire non-uniform coefficient 
as opposed to the Chester expansion formulation in (4.104) where the transition 
function operates only on the reflection term of the non-uniform coefficient. Thus, 
near the RSB and away from the caustic, (4.119) would be in agreement with the 
UTD diffraction coefficient only to the first order term. Away from transition regions, 
however, it reduces to the exact expression for the half-plane diffraction coefficient as 


well. 


4.3 Uniform asymptotic analysis (3— D case) 

The surfaces considered for the 3 -D analysis are assumed to be rotationally symmetric 
or slowly varying in the transverse direction, with the 2 — D polynomial boundaries of 
Figure 4.2 in the previous section being the generator curves representing one of 
the principal coordinate curves of the 3-D surface. The general procedure employed 
in this section is very similar to the one employed in Chapter 3 for the near-zone 
problem. Hence, some of the details will be omitted. 

As usual, the electric field E’ scattered from a 3-D polynomial surface containing 
an edge illuminated by a plane wave and observed at a point P in the far-zone of the 
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surface is expressed in terms of the radiation integral over the electric current J on 
5, that is 

E‘(P) a ^ fj [j x } x J(Q’)] ^ (4.120) 

s 

where 


R = 


Q' is any point on 5, r 0 = 


Q'P j w r a — f(Q') • S', (far-zone approximation) (4.121) 

S = * sin 6„ cos <f„ + y sin Q 0 sin <p„ 4- z cos 0 O , (4.122) 

s' = — x sin ^cos <p' o — y sin 9' 0 sin <p' 0 — z cos^ , (4.123) 

is the distance between the origin of the reference 


OP 


coordinate system and the observation point, r(Q') = OQ’ is the surface vector, 
and S’, S are the directions of incidence and observation, respectively. The PO 
approximation to (4.120) is employed once again using the GO currents in (2.2) with 
the incident magnetic field under the plane wave assumption given in (4.18). Thus, 
the far-zone scattered electric field may be written as follows: 


1-.) ~ e (*>:,* o,<Po) 


-jkr 0 


(4.124) 


where R ~ r a was used for the magnitude term in (4.120). The vector quantity 
e(6' o , <p' o , 6 a , i p D ) is the angular field dependence and is given by the following stationary 
phase integral: 


oo oo 

?(«>:, *.,»>.)« 1 1 F(e',r')e’^ e ^dtdr', 


(4.125) 


la" oo 


where 


F(l\r') = ^ {s x s x [n(l,T') x r x p ]} , 

Ztt 


(4.126) 

(4.127) 


n(/',r') is the outward unit normal to the surface at Q', and t, r' are the principal 
coordinate curves of the 3-D polynomial surface S. The edge is located at (£ q ,t') and 
therefore extends along the principal coordinate curve r' which is assumed to have a 
slowly varying curvature. The next step in the procedure is the uniform asymptotic 
evaluation of e(6' o ,<p' 0 ,6 0 ,tp 0 ). 
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4.3.1 Uniform asymptotic evaluation of e (O' o , (p' o , 0 O , (p a ) 

First, it is assumed that the 3-D polynomial surface contains an isolated zero- 
curvature line that satisfies the following relation: 

C ( T ') = W>r')L tp = °’ (4.128) 

and is arbitrarily close to the edge. Also, the phase function in (4.127) contains at 
least one real and a number of complex stationary phase points that are situated near 
the zero curvature line and satisfy the stationary phase condition given in (3.8). The 
corresponding zero- curvature diffraction point is a stationary point along the zero- 
curvature line, it has coordinates Q p = (l p) r a ) and it can be arbitrarily close to the 
edge diffraction point with coordinates Q a = (l a ,r a ). The zero- curvature diffraction 
point satisfies the following phase condition: 

= r tt {l p ,T,) = 0 and 4 > r = f p • (?-£*) = 0, (4.129) 


where t p is the unit tangent along the zero- curvature line at Q p . The second condition 
in (4.129) is analogous to the law of edge diffraction, and both conditions in (4.129) 
must be satisfied simultaneously at the zero curvature point. If / is the order of the 
zero curvature point, that is Q p = (£ p , r„) is an /th order root of r u (t, t.) = 0, we 
have that 


d n _ 

r 

dl 


(O') 


(/p,T a ) 


= 0 ; 


2 < » < Z + 1, 


(4.130) 


and <f>(t,T') contains l + 1 stationary phase points. When l is even, the surface near 
the zero-curvature point along the t principal coordinate is strictly concave or convex, 
and one of the stationary phase points is real and corresponds to a single specular 
reflection point. When / is odd, the surface near the zero-curvature point along the 
l’ principal coordinate is concave- convex, and two of the stationary phase points are 
real and correspond to a pair of specular reflection points for observation directions 
in the lit side of the caustic. In the dark side of the caustic all stationary phase points 
are complex. 

The 3-D analysis can be effectively reduced to a 2-D situation using the same 
procedure employed in Chapter 3, that is a local quadratic approximation is initially 
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used for the phase function and the integral along the r' principal coordinate is 
evaluated using the stationary phase method. Thus, e(6' o ,<p' o ,0 o ,<p o ) is given by the 
following expression: 


where 


oo 


nd' _ \ [jk sxsx [n(/',{.) x s’ x p] 

‘ ’ J " Vfc ’ 


= r. 


^ T,) , and 


2 *rr(/',T.) 
4MO.) ^ OVf. 


(4.131) 

(4.132) 

(4.133) 

(4.134) 

(4.135) 


The asymptotic reduction of the integral in (4.131) would follow the same steps as the 
asymptotic reduction of the integral in (4.23). Therefore, the following transformation 
is introduced in (4.131): 


W, t.) = V>(0 = a + fit + — , 

1/ 

where v = / + 2, and 


(4.136) 


a = = *(£,)• (?- ?*’) 

dt 

dt 


0 = i P -(s~r) 


(4.137) 

(4.138) 


and £ p is the principal unit tangent at Q p along the t surface coordinate. Hence, 
using (4.136) the field angular dependance is given by 


where 


e(0' o ,<p' o ,0o,<Po) *e iha J </(*) e’ fc( * + ‘ 1 ' / *' ) dt , (4.139) 

ta 



7 •* dt 

dt 

* 

ta = U a 

e. • (• - *■) Tt 

t = t a 

t = 0- 


(4.140) 
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(4.141) 


u a 

<7(0 

dt 

dt 

m(v, u) 


f +1 if Q p 3S 
| —1 if Q P G S, 


G(t,r.) 


dt 
dt ’ 


dt dv _ (u — 1 )! 
dv dt m(v,v ) 
f (y) (v p ) •(?-?') 


1 

1 / 


if'(v)r 


(4.142) 

(4.143) 

(4.144) 


and r(v ) is the parametric generator polynomial boundary of the 3 D surface, 
or a parametric representation of the t principal coordinate curve. The proper 
branches for / 3 , t a and dt/dt depend on the sign of m(v p ,v) and as usual “oo” = 
oo exp(— jarg(d£'/ dt). The asymptotic evaluation of (4.139) is then carried out by 
employing the generalized Chester expansion [22], and in exactly the same manner 
as was done for the 2 — D case in the previous section. Without any further details, 
the field angular dependence is given by the following expression: 


e(01,¥>l,0 o ,¥>°) 




Vk 


i 0 K u ( 1 ,C a ^) + e ik ^- ) 




+ 


j\/5t(7+c«r l J 




eiKiV'o’ 6 »,<Po) ~ •^/=e 0 A'‘(7,C,^) + e' 


kip(t a ,T.) 

; if m(v p , */) > 0 , 

jG(L,r.) 


(4.145) 




jy/k{~f + CY J 


_klj)(l a ,T,) 

; if m(v p , v) < 0 , 


(4.146) 


where K u ( 7 ,^, 1 /) i s the generalized incomplete Airy integral [22] defined in (4.54). 
The parameters 7 , and e * 0 are given by the following expressions: 


7 

u p 




w(v, 1 /) 


U p i p - (? - 8 •) w{v p , V ) , 

sgn[m(v p ,i/)] , 

u a k l ' v \ l a - (?-?’) w(v a ,v)-t p - (?-?‘) w(» M »')| 5ir > 
(*,-!)! ‘ 


m 


( v p ^) 


v - 1 


(_7)a(* Q \ 1 

V- 1 S * \ \4>{lk,T.)\ 


(4.147) 

(4.148) 

(4.149) 

(4.150) 

(4.151) 
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where (4,r,), k = 1,2, ■ • • ,v - 1 axe the stationary phase points of Also, 

using (4.127) and (4.132)-(4.134) we have the following relationships: 


Tp{L,T.) 

= <t>(l a ,T.) = r(Q a )-(s d - s') , 

(4.152) 


= ^(^jo 7 ".) = — ^ K i(,Qk) cos 6 k , 

(4.153) 

G(t k ,T.) 

l~jk cos 8[p- R (Qk) 

» 2X sJ\ir.(C,T.) | 

(4.154) 

4>rr{tk,T.) 

= -2 k 2 (<5*) cos <?;, 

(4.155) 

G(l a ,r.) 

= jkj)- Dnu {Qa)\fp{, 

(4.156) 

pi 

R a sin 2 d a 
n c • (s d — s ‘) ’ 

(4.157) 


where Z)„ u (Q a ) is non-uniform PO dyadic edge diffraction coefficient given in (3.68) 
and (3.69), p d is the distance between the caustic at Q a and the second caustic of the 
edge diffracted ray, R a > 0 is the edge radius of curvature, i? a is the usual angle of 
diffraction, and n c is the unit edge normal. 

Next, the scattering mechanisms in (4.145) and (4.146) are separated using the 
properties and complex plane topology of in Appendix C. 

4.3.2 Reflected field contribution 

There exist two cases for the reflection from a three dimensional polynomial surface 
as well. For the first case we consider a polynomial surface that contains an odd order 
zero-curvature line with v odd in (4.145) and (4.146). The reflected field consists of 
two specular contributions in this case, that only exist on the lit side of the caustic. 
For the second case we consider a polynomial surface that contains an even order 
zero- curvature line with v even in (4.145) and (4.146). The reflected field in this case 
consists of a single specular contribution that exists on both sides of the caustic. 

A. Reflection near an odd order zero-curvature line 

Using (4.124), (4.145), (4.146) and the appropriate expression for with v 

odd and 7 < 0 in Appendix C, the reflected field on the lit side of the caustic is given 
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by 


E r (0' o ,(p' o ,6 o ,yj 0 ;r 0 ) 


{«. ') 



+ 

g ~jkr 0 

u -2 e J *“ 2 e 2 I T ( 7 , ,u)} ; 

T o 

if m(v p ,i/) < 0, (4.158) 

E T (O'oi<Po>0o,<p o -,r o ) 


R e J *°* e, /*(— 7 l , i/) 



+ 

g-jfcro 

u 2 e J *° 2 e 2 / r ( 7 2 ,t/)} ; 

if m(v p ,u) > 0, (4.159) 


where u 1>2 = 1 if Q l 7 € 5 and zero otherwise, the function is defined in (4.57), 

a 12 , 7 i ,2 are the same as in (4.58) and (4.59), and e 12 for the 3-D case is given by 

P' R ( ^ i , 2 ) j 2 » (4.160) 

where R ($, i2 ) is the GO dyadic reflection coefficient given in (3.54) and ( 3 . 55 ), 
and /?i, 2 (Qi, 2 ) are the principal surface radii of curvature at <? lj2 . Thus, using (4.58) 
and (4.160) the reflected field from a polynomial surface that contains an odd order 
zero-curvature line for the case of Ri(Q { ) < 0 , and R^Qi) > 0 may be written as 
follows: 


e,- = e-’f < 


c l,2 


E T {B' o ,ip' o ,d o ,ip 0 -,r 0 ) ~ 


~ e { (Qi )• n'iQ^y/MQiWQi) eJkf(Q> ),r 


L 2 


+ y &&,)• n (Q,)yfKmMQ7)e^ ) >T 


-jkr Q 


(4.161) 


where E i {Q l 7 ) is the incident plane wave at Q l 7 given in (4.20), R (Q 12 ) is the UGO 
dyadic reflection coefficient given by 


n (Qi,0 = n.(Q lt2 ) e\e\ + H h (Q ia ) e\e[ , (4.162) 

and R,. <h {Qia) are the same as in (4.63) and (4.64). It is also understood in (4.161) 
that if Rt {Q i, 2 )R 2 {Qij) < 0, ^ il 1 (Q li2 )fi 2 (Q li2 ) = J^|^i(Qi, 2 )^2 2 (Qi,i)|. 


B. Reflection near an even order zero-curvature line 

Using (4.124), (4.145), (4.146) and the appropriate expression for K u { 7 ,^, 1 /) with v 
even in Appendix C, the reflected field on either side of the caustic is given by 


E T (0' o >'p' o >0o, <Poir 0 ) ~ u r e ihar e r /;(- 7 r ,i/) 
E r {0' o ,<p' o ,6 o ,ip o ;r o ) ~ u T e jkar e r / r (~ 7 r , u) ■ 


o—jkr 0 

T 0 

,-jkr 0 


; if m(Q p ) < 0, (4.163) 
; if m(Q p ) > 0, (4.164) 
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and for the 3-D case are given by 


,1 / !j=l 
e } i 


where 


KQp,") 


pi = 


sin i? p | cos <p' p -|- cos ip p \ , 

(4.169) 

> <Pp, tip-, v) [p • (-KK + c ;,e ;)] ^ , 

(4.170) 

( u “ !)• h 
m (Qp, u ) 

(4.171) 

-« ( r" s) (Qp)V(a-s*) , 

(4.172) 

w(Q p ,u) sm i? p cos ( p 2 j , 

(4.173) 

R p sin 2 i? p 
n p • (is ” - s’) ’ 

(4.174) 


n p is the outward unit normal at Q p , d p is the angle between the tangent along the 
zero-curvature line at Q p and the direction of incidence (analogous to the angle of 
diffraction,) tp p , <p p are the angles of incidence and observation, respectively, projected 
on the plane perpendicular to the tangent along the zero-curvature line at Q p , p p is 
the distance between the caustic at Q p and the second caustic of the zero-curvature 
diffracted ray, R p > 0 is the radius of curvature of the zero-curvature line, and n p 
is the unit normal of the zero-curvature line at Q p . Also, the v — 2th derivative of 
the surface curvature in (4.172) is along the the l' principal coordinate line. Thus, 
using (4.75), (4.169), and (4.170) the odd order zero-curvature diffracted field from a 
polynomial surface may be written as follows: 

- = p-jkro 

r.) ~ u. £■(<?,)• , (4.175) 


~ <*«£•(«.)■ V(Q,ypi , (4.175) 

r a 

where E‘(Q P ) is the incident plane wave at Q p given in (4.20), and V co (Q P ) is the 
EUTD dyadic odd order zero- curvature diffraction coefficient for oblique incidence 
given by 


V‘°(Q p ) = V:°( l p' p ,' Pp ,4 p] v)e i ± e’’ ± + , 

„ v ... . f ( 


(4.176) 


T»co ( i In . v__ / P< (-7 P ,I/), (lit side) 

— T . q nr— ^pyPp’i ^p? ^p? V) < 

sm t? p v 27 T IJLw). (dark side) 

(4.177) 
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where u r = 1 if Q r G 5 and zero otherwise, J r ( 7,1/) is the same as in (4.57) with u 
even, a r , 7 r are the same as in (4.67) and (4.68) and e r for the 3-D case is given by 


Cr 


= < 


\v- 1 

2tt 


7r 


2(^-0 


?• («.)] • (4-165) 


Thus, using (4.67) and (4.165) the reflected field from a polynomial surface containing 
an even order zero- curvature line for the case of Ri(Q r ) > 0 may be written as follows: 


£’(«>:, ~ y E‘(Q r y K (<?,) *' 




(4.166) 


where ^’(Q r ) is the incident plane wave at Q r given in (4.20) and % (Q r ) is the UGO 
dyadic reflection coefficient and is the same as in (4.162), (4.63) and (4.64) with v 
even. When Ri(Q r ) < 0 the complex conjugate version of jl (Q r ) should be used. 


4.3.3 Zero-curvature diffracted field contribution 

As in the case of the reflected field, we also have two cases for the zero-curvature 
diffracted field. The first case involves diffraction from a zero-curvature point associ- 
ated with an odd order zero- curvature line in a polynomial surface. The reflected field 
is discontinuous across the caustic in this case, therefore the zero-curvature diffracted 
field has to compensate this discontinuity. The second case involves diffraction from 
a zero-curvature point associated with an even order zero-curvature line in a poly- 
nomial surface. In the latter case no discontinuities have to be compensated by the 
zero-curvature diffracted field since the reflected field is continuous across the caustic. 


A. Diffraction from an odd order zero- curvature point 

Using (4.124), (4.145), (4.146), and the appropriate expression for K v) f° r v 
odd and £ < 0 in Appendix C, the odd order zero-curvature diffracted field is given 


by 


E p (0' o ,<p‘ o ,e o ,<p o -, ro ) 

£^>:A,*>o;r 0 ) 


oS kr o 


u 0 e ikap e p I pl (- 7 P , v) , (fit side) 


r a 

jkr Q 


u 0 e Jkap e p / pd (7„, v) , (dark side) 


(4.167) 

(4.168) 


where u 0 = 1 if Q p € S and zero otherwise, the functions I pt ( 7, »/), I pd ( 7, v) are defined 
in (4.73) and (4.74), respectively, a p is the same as in (4.75), and the parameters Ip 
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Notice that when v = 3, I p t( — 7 P ,3) = 0 and the EUTD expression for the zero- 
curvature diffracted field is equivalent to the complex ray field solution derived in 
Chapter 3 specialized for the case of plane wave incidence with the observer in the 
far-zone. 


B. Diffraction from an even order zero-curvature point 

Using (4.124), (4.145), (4.146), and the appropriate expression for K u ( 7,C>^) for v 
even and £ < 0 in Appendix C, the even order zero-curvature diffracted field on either 
side of the caustic is given by 


E P {0o><Po>8o,<Po-,r o ) 

Er(9' 0 ,<p' 0 ,e 0 ,ip 0 -,r 0 ) 


g-J" r o 

~ «0 e ikap Cp J P { lp,v) ; if ™(<3 p> *0 < 

To 

g j ^ ro 

~ u 0 e ikap e p J p (7p, v) ; if rn(Q p , v) > 0. 

r a 


(4.178) 

(4.179) 


The parameters u 0 , a p , 7 P , and e p are the same as in the odd order zero-curvature 
diffraction case, and the function J p ( 7,1/) is defined in (4.86). Thus, using (4.75), 
(4.169), and (4.170) the even order zero-curvature diffracted field may be written as 
follows: 


E p (8' o ,<p' o ,6 0 ,<p oiO ~ tic E\Q p y V ce (Q P )J P p e 




o-i kr o 


(4.180) 


where E^Qp) is the incident plane wave at Q p given in (4.20), and V C ‘(Q P ) is the 
EUTD dyadic even order zero-curvature diffraction coefficient for oblique incidence 
given by 


V"(Q P ) = V\\<p' p ,v P ,ti p \v)e\e p ± + V c ;(v' p ,<p p ,d p \v)e\f\, 


(4.181) 


'fk# 




J p*(7p>» / ), rn(Qp^) < 0 

m(Q p ,u)>0. 

(4.182) 


4.3.4 Edge diffracted field contribution 

As in the 2-D case of the previous section, the two methods of analysis employed in 
Chapter 3, namely Chester’s expansion and the method of steepest descent, result in 
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two different expressions for the edge diffraction coefficient that are in agreement to 
the first order term, however, for observation points near the RSB. They also reduce 
to the exact expression for the half-plane diffraction coefficient for oblique incidence 
outside the transition regions. Let’s consider Chester’s expansion method first. 


A. Chester’s expansion formulation for the edge diffracted field 

Using (4.124), (4.145), (4.146), and the appropriate expression for K u ( 7,C, I/ ) for in 

Appendix C, the edge diffracted field may be written as follows: 


Et(0:,<p' o ,e B ,ip o ;r o ) ~ le^e a I a ( la ,C a ,^ + e ik ^^ 


+ 


J(7+V')]} 




klp(l a ,T.) 


9 -jkr a 


; if m(v p , u) > 0, (4.183) 


Et(9' o ,<p’ o ,0 o ,<p o] r o ) ~ \ei ka ‘e <i r a { la ,C,v) + e ik *^ 


jG{l a ,T.) 


krl>(l a ,T.) 

TeTfCro] } ] if < °’ (4 - 184) 


where a 0 is the same as in (4.92). From (4.147)-(4.151), along with some algebraic 
manipulations, the parameters 7 0 , £ 0 , and e a for the 3-D case are given by 


7a = U p k (l/ l)/v sin i? p (cos <p' + cos <p p )w(v p ,i /) , 


(4.185) 


C« = u a k’ /v sini? 0 (cos<^ + cos ip a )w(v a ,v) — sim? p (cosyjp + cosy? p )ti;(v p ,i/)| t ' 1 , 

(4.186) 

(4.187) 




W(ta,T.) 

where u p and w(v,u) are given in (4.148) and (4.150), respectively, 

f(L,T.) « ^ [p- (-elei + e jej)] y/p* cos , (4.188) 


and the function J„(7,C> I/ ) 1S defined in (4.97)-(4.99). Thus, using (4.152), (4.156), 
(4.157), (4.92), (4.187), (4.188), and the PO correction factors in (2.117) and (2.118), 
the edge diffracted field may be written as follows: 






a e 


~)kr 0 


(4.189) 
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where V\(Q a ) is the EUTD dyadic edge diffraction coefficient for oblique incidence 
derived using the Chester expansion method and is given by 


V\{Q a ) = + V*(<p , a ,<p a -,d a )e\£\, (4.190) 


— e 


- 2 sin „\/2jrft 




COS 


•7 r (7<oCd,*') 


cos 


(4.191) 


7 d , ( d are the same as in (4.105), (4.106), and is defined in (4.107). It is 

understood that when R, (Q a ) < 0, the complex conjugate version of J-(j, £, v) should 
be used in (4.191). 


B. Steepest descent formulation for the edge diffracted field 

For the steepest descent analysis one begins with the expression for the scattered field 
in (4.124) and (4.139), that is 


E‘(9' 0 , y 9 e , <p a ; r„) « e Jka 1.(0, t a , u; k ) 


e~ jkr ° 


(4.192) 


where a, j3, and t a are given in (4.137), (4.138), and (4.140), respectively, and 
l,(0,t a ,u]k ) is a stationary phase integral in the complex t-plane given by 


1.(0, t a , w,k) = J Q(t) e jk(f>t+t ‘' / “ ) ds . (4.193) 

to 


/,(/?, < 0 , v\ k) has the same complex plane topology as the stationary phase integral 
in (4.109) for the 2-D case. Therefore, the edge diffracted field contribution in (4.192) 
is identified in the same fashion as in the 2-D case, and without further details it 
may be written as follows: 


El(9’ o ,<p' o ,9 0 ,<p 0 -,r 0 ) 


EWo ,<P' 0 ,9 0 ,<p 0 -,r 0 ) 


~ e 


,jka a 


G(l„r,) 

kj)(l a ,T,) 

if m(v p , v) > 0 , 

„ I± 

if m(v p , I/) < 0 , 


(7a +C ‘) L(la,(a,v) 


D -jkr Q 


(7a +C ') ^’(7a,Ca,^) 


-jkr Q 


(4.194) 


(4.195) 


where a a is given in (4.92), 7 a , £ a are the same as in (4.185) and (4.186), and I a ( 7 , £, u) 
is defined in (4.97)-(4.99). Finally, using (4.156) and the PO correction factors 
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in (2.117) and (2.118), the steepest descent formulation for the edge diffracted field 
for the case of Ri{Q a ) > 0 is given by 




n d p i^(Qa)’* d 
r a 


0 -i kr o 


(4.196) 


where p d a is given in (4.157), and X>‘(Q a ) is the EUTD dyadic edge diffraction coeffi- 
cient for oblique incidence derived using the method of steepest descent and is given 
by 


X K(Q a ) = 1 + ^>.i^)^;, (4-197) 


-e-*i 


2 sin 4 a y/tok 


COS 


(^) cos(^)J 




(4.198) 


where 7*, ( d , and ^ r (7,C) I/ ) are the same as in (4.105)-(4.107). When R^Qa) < 0 the 
complex conjugate of •7 r (7,£,i') should be used in (4.198). 



Chapter 5 


Numerical Results and Discussion 


In this chapter, numerical results for the scattered fields from cubic and fourth order 
polynomial strips illuminated by a plane wave are presented and discussed. The ac- 
curacy and limitations of the UGO/EUTD solution derived in this report are assessed 
by comparison with an independent method of moments (MM) solution. Results ob- 
tained using the classic GO/UTD solution are also shown to illustrate its failure near 
caustics and composite shadow boundaries and the need for the new solution. The 
total scattered field is plotted in terms of echo width which is defined as follows: 


L e (6,6') = lim 

p— ► oo 



(5.1) 


Both monostatic and bistatic results for the plane wave scattering from cubic and 
fourth order polynomial strips are presented and discussed in Sections 5.1 and 5.2, 
respectively. 


5.1 Scattering results for cubic polynomial strips 

In this section, results for the plane wave scattering from cubic polynomial strips 
are presented and discussed. The method denoted as UGO/EUTD1 in the results 
that follow uses the edge diffraction coefficient in (4.104) derived using the Chester 
expansion. The method denoted as UGO/EUTD2 uses the edge diffraction coefficient 
in (4.119) derived using the method of steepest descent. The scattering geometry and 
the relevant parameters are illustrated in in Figure 5.1. Before proceeding with the 
results for the total scattered field, the various field mechanisms are examined and 
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Figure 5.2: Scattered field contributions (TM case) from a cubic polynomial strip 
with a 0 = 2.0A, a, = 0.5, a 2 = 0.1 A -1 , a, = 0.1A -2 , a = — 1.5A, b = 1.5A, and angle of 
incidence O' = —45°. 


the EUTD1, EUTD2, and UTD edge diffraction coefficients are compared. The first 
example considered is a cubic polynomial strip with a 0 = 2.0A, a, = 0.5, a 2 = 0.1A" 1 , 
a 3 = 0.1A~ 2 , a = — 1.5A, and b = 1.5A. In this case both edges are far removed from 
the zero-curvature point, and thus the two RSBs and the caustic of the reflected rays 
are clearly distinct. Figure 5.2 shows a plot of the magnitude of the various field 
contributions to the total scattered field due to a TM polarized plane wave incident 
at an angle O' = —45°, computed using the UGO/EUTD1 method. Notice that 
the two reflected field components exhibit a combined total of three discontinuities. 
The first discontinuity occurs across the caustic of the reflected rays (0 fa —85°) 
and is compensated by the zero-curvature diffracted field that only exists in the 
dark side of the caustic. The second discontinuity occurs at the RSB associated 
with the edge at Q a (6 fa —52°) and is compensated by the edge diffracted field 
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Figure 5.3: Edge diffracted field magnitude (TM case) from a semi-infinite cubic 
polynomial boundary with a 0 = 2.0A, a, = 0.5, a 2 = 0.1A -1 , a 3 = 0.1A -2 , a = — 0.33A, 
and angle of incidence 8' = —45°. 

from Q a . Similarly the third discontinuity occurs at the RSB associated with the 
edge Q b (8 m -23°) and is compensated by the edge diffracted field from Q b . Also 
notice that the two edge diffraction terms become singular near the ISB ( 6 = 135°); 
however, they combine to give a finite result. The second example considered is 
semi-infinite cubic polynomial boundary with flo = 2.0A, a, = 0.5, a? — 0.1A , 
a, = 0.1A -2 , a = -0.33A. In this case the edge at Q a coincides with the zero- 
curvature point, and thus a composite shadow boundary is formed. Figure 5.3 shows 
a plot of the magnitude of the edge diffracted field due to a TM polarized plane 
wave incident at an angle 8' = — 45°, computed using the EUTD1, EUTD2, and 
UTD diffraction coefficients. Notice that UTD produces a singular result near the 
composite shadow boundary. The EUTD1 and EUTD2 results are not in complete 
agreement near the composite shadow boundary, however they both remain finite. All 
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Figure 5.4: Edge diffracted field magnitude (TM case) from a semi-infinite cubic 
polynomial boundary with a 0 = 2.0A, dj = 0.5, a 2 = 0.1A -1 , a 3 = 0.0A~ 2 , a = — 0.33A, 
and angle of incidence O' = —45°. 

three diffraction coefficients seem to be in good agreement away from the composite 
shadow boundary. The third example considered is semi-infinite parabolic boundary 
which is a special case of a cubic polynomial boundary with a 0 = 2.0 A, a, = 0.5, 
a 2 = 0.1A -1 , a 3 = 0.0A -2 , a = — 0.33A. In this case no caustic or composite shadow 
boundary exists and UTD is expected to give accurate results for the edge diffracted 
field. Figure 5.4 shows a plot of the magnitude of the edge diffracted field due to 
a TM polarized plane wave incident at an angle O' = —45° computed using the 
EUTD1, EUTD2, and UTD diffraction coefficients. All three diffraction coefficients 
are in good agreement, however, the EUTD1 diffraction coefficient seems to be in 
better agreement with UTD near the RSB. Figure 5.5 shows the bistatic echo width 
of a cubic polynomial strip with a 0 = 2.0A, a l = 0.5, a 2 = 0.1A -1 , a s = 0.1A -2 , 
a = — 0.33A, and b = 1.5A illuminated by a TM polarized plane wave incident at an 
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Figure 5.5: Bistatic echo width (TM case) of a cubic polynomial strip with u 0 = 2.0A, 
a, = 0.5, a 2 = 0.1A -1 , a, = 0.1 A -2 , a = — 0.33A, b — 1.5A, and angle of incidence 
8' = -45°. 

angle 8' = —45°. The classic GO/UTD solution produces a singular result near the 
composite shadow boundary whereas the UGO/EUTD1 and UGO/EUTD2 remain 
finite. Notice that the UGO/EUTD1 solution is in better agreement with MM near 
the composite shadow boundary whereas the UGO/EUTD2 solution seems to be more 
accurate away from the composite shadow boundary. The differences, however, are 
quite minor (< 1 dB-A) and therefore the two EUTD diffraction coefficients should 
be considered equivalent for all practical purposes. For the subsequent results only 
the UGO/EUTD1 results will be shown. 

Figures 5.6 and 5.7 show plots of the bistatic echo width of a cubic polynomial 
strip with a 0 — 2.0A, a, = 0.5, a 2 = 0.1A -1 , a, = 0.1A -2 , a = — 1.5A, b = 1.5A, and 
angle of incidence 8' = —45°, for the TM and TE polarization cases, respectively. 
The UGO/EUTD result shows excellent agreement with MM for the TM polarization 
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Figure 5.6: Bistatic echo width (TM case) of a cubic polynomial strip with a 0 = 2.0A, 
Oj = 0.5, a 2 = 0.1A -1 , o 3 = 0.1A -2 , a = — 1.5A, b = 1.5A, and angle of incidence 
d' = -45°. 


115 




Figure 5.7: Bistatic echo width (TE case) of a cubic polynomial strip with a 0 = 2.0A, 
a t = 0.5, a 2 = 0.1A -1 , a, = 0.1A -2 , a = -1.5A, 6 = 1.5A, and angle of incidence 
B' = -45°. 
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case. For the TE polarization case, the higher order mechanisms which are not 
included in the UGO/EUTD solution are quite significant and therefore the agreement 
is reasonable only near the specular and forward directions where the first order 
mechanisms are dominant. For observation directions near grazing the higher order 
mechanisms such as edge excited creeping waves and whispering gallery modes excited 
in the concave part of the cubic strip must be included to obtain accurate results for 
the TE polarization case. The classic GO/UTD solution gives an erroneous result 
near the caustic ( 6 85°), and on both its lit and dark sides. This is expected since 

the classic UTD solution does not include zero-curvature diffraction information, and 
the GO expression for the reflected field gives singular results near caustics. Away 
from the caustic, the UGO/EUTD and GO/UTD results are in good agreement for 
both polarizations. Figures 5.8 and 5.9 show plots of the monostatic echo width of 
the same cubic polynomial strip for the TM and TE polarization cases, respectively. 

The monostatic UGO/EUTD result shows excellent agreement with MM for the 
TM polarization case. For the TE polarization case, the agreement is good only near 
broadside where the first order mechanisms are dominant. The classic GO/UTD 
solution gives an erroneous result near the caustics (0 ~ —67° and 6 ~ 113°), and on 
both the lit and dark sides. Away from the caustics, the UGO/EUTD and GO/UTD 
results are in good agreement for both polarizations. 

Figure 5.10 shows a plot of the bistatic echo width of a cubic polynomial strip 
with a 0 = 2.0A, a i = 0.5, a 2 - 0.1A -1 , a 3 = 0.1A -2 , a = — 0.33A, b — 1.5A, and angle 
of incidence 6' = -45°, for the TM polarization case. For this geometry, the edge at 
Q a coincides with the zero- curvature point, thus the RSB and the caustic coalesce to 
form a composite shadow boundary. The GO/UTD result exhibits a singularity at 
the RSB associated with the edge at Q a (0 ~ —85°), whereas the UGO/EUTD result 
remains finite and is in excellent agreement with MM. Figure 5.11 shows a plot of 
the monostatic echo width of the same geometry. Again, the non- uniformity of the 
GO/UTD solution is clearly evident, whereas the UGO/EUTD solution remains valid 
across the composite shadow boundary and uniformly reduces to the the GO/UTD 
result away from the transition region. 
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Figure 5.8: Monostatic echo width (TM case) of a cubic polynomial strip with a, 
2.0A, a, = 0.5, a 2 — 0.1A -1 , a s = 0.1A -2 , a = — 1.5A, and b = 1.5A. 
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Figure 5.9: Monostatic echo width (TE case) of a cubic polynomial strip with a 0 = 
2.0A, a t = 0.5, a 2 = 0.1A -1 , a, = 0.1A“% a = — 1.5A, and b = 1.5A. 
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Figure 5.10: Bistatic echo width (TM case) of a cubic polynomial strip with a 0 = 2.0A, 
a, = 0.5, a 2 = 0.1 A -1 , a s = 0.1A" 2 , a = -0.33A, b = 1.5A, and angle of incidence 
6' = -45°. 
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Figure 5.11: Monostatic echo width (TM case) of a cubic polynomial strip with 
= 2.0A. a, = 0.5, a 2 = 0.1A~ 1 , a 3 = 0.1A" 2 , a = -0.33A and b = 1.5A. 




Figure 5.12: Bistatic echo width (TM case) of a cubic polynomial strip with a 0 = 2.0A, 
o, — 0.5, a 2 = 0.1A l , a, = 0.0A -2 , a = — 0.33A, b = 1.5A and angle of incidence 

O' — 


The remaining four results in this section illustrate the flexibility and universal 
nature of the UGO/EUTD solution. Figures 5.12 and 5.13 show plots of the bistatic 
and monostatic echo width, respectively, of a cubic polynomial strip with a 0 = 2.0A, 

a ‘ = °‘ 5 ’ a ’ = °- 1A-1 > = 0.0A- 2 , a = -0.33A and b = 1.5A. This geometry 

of course corresponds to parabolic screen and the zero- curvature point theoretically 
moves to negative infinity. It is well known that the classic GO/UTD solution gives 
accurate results for these type of geometries. The UGO/EUTD solution remains valid 
and is in excellent agreement with both GO/UTD and MM. 

Figures 5.14 and 5.15 show plots of the bistatic and monostatic echo width, respec- 
tively, of a cubic polynomial strip with o 0 = 2.0A, a, = 0.5, a* = 0.0A-', a, = 0.0A" 2 , 
a = -0.33A and b = 1.5A. This case corresponds to a flat strip since both the cubic 
and quadratic coefficients are set equal to zero. The UGO/EUTD solution remains 
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Figure 5.14: Bistatic echo width (TM case) of a cubic polynomial strip with a 0 = 2.0A, 
a, = 0.5, a 2 = 0.0A -1 , a 3 = 0.0A -2 , a = — 0.33A, b - 1.5A and angle of incidence 











Figure 5.16: Geometry and relevant parameters for the scattering from a fourth order 
polynomial strip. 

valid for this special case also, a fact that demonstrates its flexibility for treating 
general surfaces that are highly curved, slightly curved, or completely flat. 

5.2 Scattering results for fourth order polyno- 
mial strips 

In this section, results for the plane wave scattering from fourth order polynomial 
strips are presented and discussed. The scattering geometry and the relevant param- 
eters are illustrated in Figure 5.16. The results in this section will be limited for the 
case where the two first order zero-curvature points, Q pl and Q p2} are sufficiently far 
apart. When Q pl and Q p2 are close together, the two caustics form a cusp and the 
UGO/EUTD solution would fail. The UGO/EUTD solution does apply when Q pl 


126 




Figure 5.17: Bistatic echo width (TM Case) of a fourth order polynomial strip with 
a 0 = 1 . 0 A, a, = 0 . 5 , a 2 = - 0 . 4 A- 1 , a 3 = O.lA-% a 4 = 0 . 1 A~ S , a = -1.5A, b = 1 . 0 A, 
and angle of incidence 6' = —45°. 


and Q p 7 coalesce to form a second order zero- curvature point, however, no results for 
this case are available at this time. 

Figures 5.17 and 5.18 show plots for the bistatic echo width of a fourth order 
polynomial strip with a 0 = 1 . 0 A, a { — 0.5, a 2 = — 0.4A ’,03 = 0 . 1 A , a 4 = 0 . 1 A , a 
- 1 . 5 A, and b = 1 . 0 A, illuminated by a TM and TE polarized plane wave, respectively, 
incident at an angle 6' = -45°. The UGO/EUTD result shows good agreement with 
MM for the TM polarization case. For the TE polarization case, the agreement is 
good only near the specular and forward directions where the first order mechanisms 
are dominant. The classic GO/UTD solution fails dramatically near the caustics 
( 1 0 « -112° and 6 « -35°). Figures 5.19 and 5.20 show plots for the monostatic echo 
width of the same geometry for the TM and TE polarization, respectively. Again, 
the UGO/EUTD result shows good agreement with MM for the TM polarization case, 
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Figure 5.18: Bistatic echo width (TE Case) of a fourth order polynomial strip with 
a 0 = 1.0A, a l = 0.5, a 7 = -0.4A- 1 , a, = 0.1A" 2 , a, = 0.1A-*, a = -1.5A, 6 = 1.0A, 
and angle of incidence 6' = —45°. 
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Figure 5.19: Monostatic echo width (TM Case) of a fourth order polynomial strip 
with a 0 — 1.0A, a! = 0.5, a 2 = — 0.4A -1 , a s = 0.1A % a 4 = 0.1A 3 , a = — 1.5A, 
b = 1.0A. 
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Figure 5.20: Monostatic echo width (TE Case) of a fourth order polynomial strip 
with a 0 = 1.0A, a, = 0.5, a 2 = -0.4A" 1 , a 3 = O.IA" 2 , a 4 = 0.1A"*, a = -1.5A, and 
b = 1.0A. 
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Figure 5.21: Bistatic echo width (TM Case) of a fourth order polynomial strip with 
a 0 = 1.0A, a , = 0.5, a 2 = -0.4A" 1 , a, = O.IA'% a 4 = 0.1A" 3 , a = -1.1A, b = 1.0A, 
and angle of incidence 9' = —45°. 

whereas the agreement is only good near broadside for the TE case. Also, notice that 
when the concave part of the strip is illuminated and whispering gallery modes are 
excited, the agreement is worse than the case where the convex part of the strip is 
illuminated. 

Figures 5.21 and 5.22 show plots for the bistatic and monostatic echo width, 
respectively, of a fourth order polynomial strip with a 0 = 1.0A, a, = 0.5, a 2 = — 0.4A -1 , 
a 3 = 0.1A -2 , a t = 0.1A -3 , a = — 1.1A, and b = 1.0A. For this geometry, the edge at Q a 
coincides with the zero- curvature point, Q pi , forming a composite shadow boundary. 
The UGO/EUTD solution shows good agreement with MM, whereas the GO/UTD 
solution fails both near the composite shadow boundary and the caustic associated 
with Q p2 . 


131 





Angle Q (degrees) 

Figure 5.22: Monostatic echo width (TM Case) of a fourth order polynomial strip 
with a 0 = 1.0A, a, = 0.5, a 2 = — 0.4A -1 , a s = 0.1A -2 , a 4 = 0.1A _S , a = — 1.1 A, 
b = 1.0A. 





Figure 5.23: Bistatic echo width (TM Case) of a fourth order polynomial strip with 
o 0 = 1.0A, a , = 0.5, a 2 = -0.4A" 1 , a 3 = 0.1A" 2 , a< = 0.1A' 3 , a = -1.1A, b = 0.6A, 
and angle of incidence 6' = —45°. 


Finally, Figures 5.23 and 5.24 show plots for the bistatic and monostatic echo 
width, respectively, of a fourth order polynomial strip with a 0 = 1.0A, a, = 0.5, 
a 2 = -0.4A" 1 , a 3 = 0.1A- 2 , a 4 = 0.1A"*, a = -1.1A, and b = 0.6A. In this 
case, both edges coincide with the zero- curvature points forming a pair of composite 
shadow boundaries. The UGO/EUTD solution remains uniform across both compos- 
ite shadow boundaries and shows good agreement with MM. 
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Figure 5.24: Monostatic echo width (TM Case) of a fourth order polynomial strip 
with a 0 = 1.0A, a, = 0.5, a 2 = -0.4A" 1 , a, = 0.1A" 1 , a 4 = O.lA-% a = -1.1A, 
b = 0.6A. 



Chapter 6 

Summary and Conclusions 


A uniform geometrical optics (UGO) and an extended uniform geometrical theory of 
diffraction (EUTD) were developed in this report for evaluating high frequency elec- 
tromagnetic (EM) fields within transition regions associated with a two and three di- 
mensional smooth caustic of reflected rays and a composite shadow boundary formed 
by the caustic termination or the confluence of the caustic with the reflection shadow 
boundary (RSB). 

The UGO is a uniform version of the classic GO. It retains the simple ray optical 
expressions of classic GO and employs a new set of uniform reflection coefficients 
involving the Fock type Airy functions that provide for a finite reflected field near a 
smooth caustic on the lit side. Table 6.1 summarizes the UGO reflection coefficients 
for the lit side of a smooth caustic. The UGO also encompasses a uniform version 
of the complex GO ray field that exists on the dark side of the smooth caustic. The 
UGO complex ray reflection coefficients are summarized in Table 6.2 and involve the 
ordinary Airy function. They provide for a finite field near the smooth caustic on 
the dark side. The EUTD is an extension of the classic uniform geometrical theory 
of diffraction (UTD) and accounts for the non-ray optical behavior of the UGO re- 
flected field near caustics by using a two-variable transition function involving the 
incomplete Airy functions in the expressions for the edge diffraction coefficients. It 
also uniformly recovers the classic UTD behavior of the edge diffracted field outside 
the composite shadow boundary transition region. The EUTD edge diffraction coef- 
ficients are summarized in Table 6.3. The UGO/EUTD expressions derived in this 
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Table 6.1: UGO reflection coefficients for the lit side of a smooth caustic (near-zone, 
arbitrary source illumination). 


Reflection 

Coefficients 

Tl. ih = Te-'f 7 ; /4 W,(— 7*) 

soft # TM 

^ for a < > acoustic surface or < > polarization 

hard TE 

/ / 

for a PEC surface. 

Parameter 

7< = [§*lto + ^)-(« , ‘ + J r )l] V * 

Canonical 

Function 

W,( 7 ) = f 

" L 23 

Fock-type Airy function [16] 

Notes 

7i\ h should be used when p r J(p r t + a r ) < 0. 


Table 6.2: UGO complex ray reflection coefficients for the dark side of a smooth 
caustic (near-zone, arbitrary source illumination). 


Reflection 

Coefficients 

n.,H = T2v^ 7 :'< e (i/3h * /2 Ai( 7 .) 

f soft 1 . ( TM 1 

for a < > acoustic surface or < > polarization 

{ hard J [ TE J 

for a PEC surface. 

Parameter 

7. = 

j*l(«o +R 0 )-(s i ‘ + s r ‘)\ 

2/3 

Canonical 

Function 

Ai( 7 ) ee ^ f e j(7 * + *’ /s) dz 
J L 7l 

Ordinary (Miller-type) Airy function [17] 
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Table 6.3: EUTD edge diffraction coefficients valid near a composite shadow boundary 
(near-zone, oblique incidence, arbitrary source illumination). 


Diffraction 

Coefficients 


-e-'f 


T) — — c * 

*' h 2sint 




K7..C) 


soft 

^ I 

► acoustic surface or < 

' TM 1 

hard J 


1 TE J 


=p for a I 
for a PEC surface 


Transition 

Function 


^(7.0 = i(7 + c , )g-(7,cy ( ’ < ' lt ’ nl 


g ; (7,£), i = 1,2,3 are the incomplete Airy functions [9] 

for 7 < 0, or 


5,V7>S/> * — 01 c 111 c II 

1 iff >(-7)‘« 
i = < 2 if C < — (~7) 1/2 * 


1 = < 


3 otherwise 

1 if c > 0 

2 if C < 0 


> for 7 > 0 


Arguments 


- - _ s< rn (L \L2/» 8in M c08 4l+ c08 .fol 
7a s g n O<U fc |L d (^, ¥ j a ;«> <1 )| 1 / 3 

> 1 ! /a Igini9o(co8i^n+co8<^o)-sintfa(co8yn+co8ya)| 1 / 2 

“ U “ |id(^.Va;rila)| 1 / 6 

+ 1 if Qo3S 
-1 ifQ 0 G5 


u a = < 


Notes 


i d is the EUTD distance parameter given in (3.110)— (3.117) 
^*(7 a> Ca) should be used when i d < 0 
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Table 6.4: UGO reflection coefficients (far-zone, plane wave illumination). 


Reflection 

Coefficients 

f soft 1 . [ TM | 

=F for a < acoustic surface or < polarization 

[ hard [ TE 

for a PEC surface. 

Parameters 

7 . = {* (a) rw.) - ?w,)i ■ (i - 

v — l -f 2, / is the order of the z-c point, Q p 

Canonical 

Function 

i»ooexp(jir/2i/) 

/ r ( 7 , v) = / e j(7 ' +I /u) dz 

* oc exp( — j3 

Notes 

7Z] h should be used when p g (Q r ) < 0. 


report are approximate ones because the surface current used in the spatial domain 
radiation integral representation for the fields was found using the PO approximation. 
Therefore, the UGO/EUTD expressions account for only first order mechanisms and 
remain valid within regions where grazing fields do not exist. 

The UGO/EUTD analysis was employed to investigate the far-zone RCS problem 
of plane wave scattering from two and three dimensional polynomial defined surfaces, 
and uniform reflection, zero-curvature, and edge diffraction coefficients were derived. 
These new UGO/EUTD coefficients involve a generalized form of the complete and 
incomplete Airy integrals that serve as canonical functions. The UGO reflection 
coefficients and the EUTD zero-curvature and edge diffraction coefficients for far-zone 
observation and plane wave illumination are summarized in Tables 6.4, 6.5, and 6.6, 
respectively. 

Numerical results for the scattering and diffraction from cubic and fourth order 
polynomial strips were also shown and the UGO/EUTD solution was validated by 
comparison to an independent MM solution. The UGO/EUTD solution was also 
compared with the classic GO/UTD solution. The failure of the classic techniques 
near caustics and composite shadow boundaries was clearly demonstrated and it was 
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Table 6.5: EUTD zero-curvature diffraction coefficients (far-zone, plane wave illumi- 
nation, oblique incidence). 


Diffraction 

Coefficients 

■».. j ,, u e .1 U-l,,"), (lit side) 

‘• h sin rfpV^27r P’ )l 

{ IpA-l P ,v), (dark side) 

soft ) . ( TM ) 

for a < > acoustic surface or < > polarization 

hard y TE J 

for a PEC surface. 

Parameters 

7 P = fc 11 ^u>((5 P ,t')sini? p |cos( l oJ ) -|-cos^ p | 

I p = w(Q p ,v) sini? p cos ( P 2 v ’ p ) 

*0 = " 2) (<3 P )(sin <fi' p + sin <p p ) 

v = l + 2, / is the order of the z-c point, Q p 

Canonical 

Functions 

poo exp( — j3ir/2v) 

I pt (-f,i/)= j e^" ,xJtzV lv) dz ; (*/ odd) 

oo exp[j( — w + 3ir/2»/)] 
p oo exp( j?r /2i/) 

= f e H ‘ ,x+z '' /l ' ) dz; (u odd) 

J oo «xp(j( n — 7r/2i/)] 

^oo exp(i?r/2u') 

J p ( 7 , i ') = 1 e j(7X+I ; (i/ even) 

J oo exp[j( ?r— 3ir/2t/)] 

Notes 

Jp( 7 P ,i/) should be used when rfi(Q p ,i/) < 0. 
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Table 6.6: EUTD edge diffraction coefficients (far-zone, plane wave illumination, 
oblique incidence). 


Diffraction 

Coefficients 


Notes 


T>' — — ~ e J -t r- zn 

M ' h 2 sin t9 a \/z7rS „( \ ' ( v>!i+ya 



soft . TM 

FF for a < > acoustic surface or < 

hard TE 

\ / V 

for a PEC surface. 


■7 r (7.C,* / ) = j( 7 + C 1 Ua(7»C 5 * / )e J< ' TC+< ‘ //,/) 


polarization 


/ ooexpfj^ot^)) 

e j(7I+ * ' v) dz 

f if [sgn(C )]— ‘(t + C"-‘) > 0 


e i( - 1, + x '' lv) dz 


V’o(^) = 


-37t/2i/ otherwise 


if i/ is odd 


C. = 1C 


sgn(C) if v is even, 

7„ = sgn[m(i; p ,i/)]A:V- sinti p (co6 <p‘ p + cos <p p )u>(v p ,v) 
C« = u a k l | sin i? p (cos ip p + cos <p p )uj(v p , v) 

- sin i? a (cos <p' a + cos y? a )u;(v 0 , v) | ^ 

f +1 if Qo 3 S 


— 1 if Qo € 5 


u>— 1)! |£ 


w ( v ^)=\^)\ V » m(u,i/)= r 

i/ = l + 2, / is the order of the z-c point, Q p 


<^*(7<o Co *0 should be used when J?,(Q a ) < 0 











shown that the UGO/EUTD results remained valid and uniformly reduced to the 
classic results away from the transition regions. 

It is believed that the work accomplished in this report is an important and timely 
first step in the development of a useful engineering solution for the interaction of 
high frequency electromagnetic waves with CAD generated surface components. Such 
components are often designed using polynomial or spline surface patches that are 
then put together to create the desired models for complex structures. It is therefore 
advantageous to use the CAD geometry information directly for RCS or antenna 
analysis. Topics for future research include the extension of the EUTD edge diffraction 
coefficients to handle wedge type geometries, and also the development of a diffraction 
coefficient for a vertex in a polynomial or spline surface patch. 
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Appendix A 

The Ordinary and Fock-Type Airy 
Functions 


Both the ordinary (Miller-type) and Fock-type Airy functions satisfy Airy’s differen- 
tial Equation [23], that is 

y"{<r) - <ry(<r) = 0 . (A.l) 


The ordinary Airy functions Ai(<r) and Bi(cr) are two independent solutions of (A.l), 
and in integral form they are defined as follows: 



(A.2) 

(A.3) 


where the contours of integration Z, 21 , I 2S and Z 3I are shown in Figure A.l. 

Their power series form which may be used for small argument approximations is 
given by [21] 


where 


Ai(cr) = cj(a) - <hg(<r), 

Bi(<r) = y/Z [c, /(<t) + Cj g{<r )\ , 


fi*) 

9(<r) 


c , = 


<h - 


, 1 s 1 • 4 1-4-7 , 

1 + 3! <r +_ 6T°' 9! a +•" 

2 , 2 - 5 , 2 - 5-8 

‘ ,+ 4! <T + 1 r" + ToT‘ r +- 
Ai( 0 ) = Bi(0)/\/3 = 3- 2/ 7r(2/3), 

— Ai'(O) = Bi / (0)/v / 3 = 3- 1/s /r(l/3). 


(A.4) 

(A.5) 


(A.6) 

(A.7) 

(A.8) 

(A.9) 


142 




Figure A.l: Contours of integration for the ordinary Airy functions. 
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The asymptotic forms of the ordinary Airy functions and their derivatives are as 
follows: 


1 . a > 1 : 


Ai(o-) 

AiV) 

BiM 

BiV) 


--1/4 


2 


-.1/4 

^ _ _ _ e -(»/»)«* /a 


2^ 


-.-1/4 


<r 1/4 


-.1/4 


(A.10) 
(All) 
(A.12) 
(A. 13) 


2. <r <C — 1: 


~ ^r 8 in (r’" + i)’ 


// V 0-1/4 /2 ... 7T\ 

**<-'> ~ -^ cos (r + j) ’ 


r-'/ 4 


, v ^' 1/4 /2 ... 7T\ 

~ + 4 )’ 


~ 7 ? 8in (r 7 ' , ’ + i) 


(A.14) 
(A. 15) 
(A. 16) 
(A.17) 


Complete asymptotic expansions for the ordinary Airy functions may be found in [21]. 

The Fock-type Airy functions W 12 (<r) [16] are also independent solutions of (A.l), 
and in integral form they are defined as follows: 


W,(<r) = 4 =/ e j(<ri+lJ/,) dz , 

V 7 T w L23 

W 2 (cr) = 4= / e*“ + *’ /s) dz = W;(rr), 

V^T ^^13 

In terms of the ordinary Airy functions they are given by 

W li2 (<r) = ypR [Bi(o-) ± j Ai(<r)] . 


(A.18) 

(A.19) 


(A.20) 


Of particular interest are the asymptotic forms of W 12 and their derivatives for a < 
-1. Using (A.14)-(A.17) and (A.20) they are given by 


W,(— <r) ~ e J fo-- 1/4 e j(2/s)ffS/2 , 


(A.21) 
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(A. 22) 
(A. 23) 
(A.24) 


W',(-<r) ~ -e _i V /4 > 
W 2 (-<r) ~ e-*io -~ lh e 
W'(-ff) ~ -e j V /4 e 


° 5 

-j(2/3)or 3 / 2 
-j(2/3)<r S / 2 _ 
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Appendix B 

Numerical Evaluation of the 
Incomplete Airy Functions 


In this appendix, a method for the accurate and efficient evaluation of the incomplete 
Airy functions is developed. The incomplete Airy functions are functions of two vari- 
ables and they satisfy the parabolic partial differential equation applied by Fock [24] 
to the study of fields near the surface of a smooth convex scattering body, that is 


‘ d 2 
pp 


0- 3 


d 

dt 


&(0,O = o ; 


i = 1,2,3. 


(B.l) 


In integral form, the solutions of (B.l) are defined as follows: 


a /*ooexp(j>,) 

g ,(0,0 = J ( e 3{e *+ z /s) dz ; i= 1,2,3, (B.2) 

where the upper limit lies within one of the three sectors of the complex 2 -plane in 
which the integral converges, that is 


2 (* “ l)f - ^ 5 *=1,2,3. (B.3) 

The contours of integration for the incomplete Airy functions are shown in Figure B.l. 

Only the function g, (/?,£) will be considered in this appendix since the other two 
functions, namely g 2 and g s , can be obtained from g, and the ordinary and Fock-type 
Airy functions, that is 

&(^0 = g.(/9,0-2irAi(/?), (B.4) 

and 


&(0.0 = MO - iV*W s (j3) . (B.5) 
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z-plane 


Figure B.l: Contours of integration for the incomplete Airy functions. 
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The arguments (3 and £ will be taken as real, since in most practical applications real 
f3 and £ are of primary interest. However, this is not a requirement for the analysis 
that follows and the resulting formulae are valid for arbitrary values of (3 and £. In 
addition, £ will be restricted to positive values since for negative values of £, g, may 
be obtained using the following expression: 

g,(/3, -£) = 27t Ai(/?) - g*(/3,0 , ( B - 6 ) 

with (*) denoting the complex conjugate operation. 

B.l Derivation of the series solution 

In order to obtain a series solution for g,(/3, £), we begin with the parabolic differential 
equation in (B.l) and assume two independent solutions of the form 

y,(/J,£) = f>n(W, ( B - 7 ) 

n~0 

and y,09,{) = n(()/3" ■ (B.8) 

n=0 

Substituting (B.7) and (B.8) in (B.l) we obtain the following expressions: 

( B ' 9 > 

n=2 n=0 n=0 

and ^n-lWfr’-EMfr'-iEWfr s 0. (B.10) 

n = 2 n = 0 n ~° 

For Equations (B.9) and (B.10) to be satisfied, the sum of coefficients of like powers 
of (3 must be zero for any value of £. To obtain the first independent solution, y,, we 
let o,(f) = 0, and setting the sum of coefficients of like powers in (B.9) equal to zero 
we obtain 


«i(0 = 

r-> 

> 

o 

(B.ll) 

a 2 (0 = 

£««). 

(B.12) 

and a„(0 = 

fln-3(0 + J a n- 2(0 ^ 0 

( i\ 5 71 - 0 * 

n(n — 1) 

(B.13) 
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Similarly, for the second independent solution, y 2 , we let MO = 
sum of coefficients of like powers in (B.10) equal to zero we obtain 

0 and setting the 

> 

© 

II 

w 

(B.14) 

6,(0 = o, vf, 

and M0 = tn ~ l( l )+, M-. "S3. 

n(ra — 1) 

(B.15) 

(B.16) 

Thus, y, and y 2 are given by: 


y.(0.O = M0 + £°n(O0 n . 

n-2 

(B.17) 

and y,(/?,0 = b,{(V + £ MfV»“ , 

(B.18) 


n= 2 


where a„(0 and &„(£) for n > 2 can be expressed in terms of a 0 (£)> MO and t ^ ieir 
derivatives, respectively, via Equations (B.11)-(B.16). Now, g,(0,O must be equal 
to the sum of the two independent solutions, that is 

g,(0,O = «o(0 + MO0 + X^[®n(0 + MO10” > (B.19) 

n = 2 

and ^zg,(0,O = + + ( B - 20 ) 

Op n= 2 

Although the coefficients a n (0 and 6 n (0 can be combined into a single coefficient, 
keeping them separate greatly simplifies their evaluation in a computer code. 

Finally, it remains to find expressions for a 0 (0 and MO in order to complete the 
solution. This can be accomplished by applying the proper boundary conditions at 
(3 = 0 using the integral form of g,(/?, 0 given in (B.2), that is 

oo 

«.(«) = < B - 21 > 

e 

and 

oo 

= f jze j ' 3/i dz. (B.22) 

0=0 e 

The functions a 0 (0 and MO can be expressed in terms of the incomplete Gamma 
function [21] as follows: 

o 0 (0 = e w ® 3 -2/ * T(l/3, — jO/3) > (B.23) 

6,(0 = 3 -l/ * T (2/3, — i0/3) , (B.24) 
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6.(0 = Jjg.(0>o 



where 


T{x,y) = J t'~ l e~‘dt, R{x)>0. (B.25) 

V 

Using the series solution of the incomplete Gamma function [21], a 0 (£) and &,(£) are 
computed using the following expressions: 

N(() (iP/W 

«.(«) « «*"* 3- ,n T(l/3) - f S ^ , (B.26) 

and 4,(0 « 3-‘" r(2/3) - if £ ■ (B.27) 

(3n + 2)n! 

The number of terms in the series, N(£), is given by the empirical formula: 

N(() = 2 + , (B.28) 

and results in a truncation error of less than 10~ 6 . The derivatives of a 0 (£) and &,(£) 
are given by: 

<K() = M -"(() + »X*‘”(0+ n>3, (B.29) 

&!*’({) = i(t-l)a'*-'>(0+>{a'‘>(f), »>1, (B.30) 

where 

<00 = , (B.31) 

and a' 2) (£) = -j?e i(i/s . (B.32) 

The constants d k and e k are obtained using the following recursive relationships: 

dk ~ -f e*.., , (B.33) 

and e k = e*_,+2, k >3, (B.34) 

with d, = e 2 = 0. 

When £ — ♦ -oo, our solution should reduce to the series solution for the ordinary 
Airy function, Ai(/3), given in Appendix A. In this case we have 

a 0 (£-> — oo) = 2ir Ai(0) = 2.23070703, (B.35) 

bi(C -► -oo) = 27rAi'(0) = -1.62621025, (B.36) 
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and using (B.11)-(B.18) we obtain 

g,(/M - -°°) = 2*Ai(0)(l+i/?> + ^,8‘+i^/3* + ---) 

+ 2xAi'(0)(/3+|^ + ^'+m^r + -) 

= 27r Ai(/3 ) , (B.37) 

which is a necessary condition for the validity of (B.19). 

B.2 Derivation of asymptotic formulae 

In this section, three asymptotic formulae for the incomplete Airy function, <7i (/?,£)> 
involving several terms are derived and serve as large argument forms. We begin by 
introducing the large parameter fi in the integral form of g,(/?,f)i examine the 
following integral: 

oo 

/ 0 (<7, 7 ; ft) = / e J ‘ n( ~ + * 3/3) dz . (B.38) 

7 

The objective is to obtain asymptotic expansions for I 0 as fi — > oo for various dispo- 
sitions of the saddle points and endpoint. Then, the asymptotic formulae for gi(/3,() 
may be obtained using the following expression: 

*(0,0 = n i/3 4(* = /W- 2/3 ,7 = (B.39) 

Since we are interested in real /3 and £ with £ > 0, the analysis that follows will be 
restricted to real a and 7, with 7 > 0. 

Asymptotic Formula for {3 + 1 

This case corresponds to the endpoint being far removed from the possibly neighboring 
saddle points of the integrand in (B.38), or 7 > M 1/J - Although the sign of <r 
is irrelevant in this case, the saddle points 2^2 = ±(— are taken as real for 
simplicity. Also, the original integration path P 0 in Figure B.2 is deformed into the 
steepest descent path, P a d p , leading away from the endpoint at 2 = 7 and into the 
sector 0 < arg(z) < tt/ 3 of the complex 2 -plane. The asymptotic evaluation of / 0 
can then be performed by means of Watson’s lemma or alternatively by repeated 
integration by parts as follows: 
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I ( a . ~.m = [ — e in ^ z3/3) dz = 3 

Jjn((T+z*)dz 




n(o- + 7 2 ) 


+ 


/ 5 


2z 


e >n ("+* a /»dz. 


(BAO) 

jn(<T + z>y'' 

T 

Integration by parts can then be applied to the second integral in (B.40), and repeat- 
ing the same procedure a few more times yields the following expression for J 0 (<r,7, 0). 

27 


I.(*,7;0) * 


(IM 

\ n [<7 + 


i 

+ & 

i 

+ & 


2 j 


7 
12J7 


1 


[(o- + 7 2 ) 4 (o’ + 7 2 ) 5 

— 40j 840j7 2 


+ 


L(°'+7 i n 

407 


1 ' 

fi 7 [(o- + 7 2 ) 8 (o - + 7 2 ) 7 J 

1680j7 


1207 s 


>jY ]\ 

7 2 ) 9 Jj 


■> « 1 v pm- 6 ). (B.41) 

[(0- + 7 2 ) 7 (0- + 7 2 ) 8 (o- + 7 2 ' 

Now, using Equations (B.39) and (B.41), the asymptotic formula for gl (/3,0 when 

£ » |/3| 1/2 is given by 


gl (/3,0 ~ e i(W+<4/S) 

40£ 

+ 


l/? + f’ + O0 + £’)* T (P + { , )‘ (0 + ?)' 

120^ - 40j 840 j^ 2 1680j^ ] 


2£ 


+ 


2j 


12jf 


(/3 + e) 6 (/3 + e) 7 03+£ 2 ) 8 (/* + £ 2 ) 

The derivative of g,(/3,0 is given by 


(B.42) 






[0 +> (/3+^) 2 + (^+^) s (^+er 

280£ + 120jT 


2jf 


8£ 


8j — 12f 100 

\S "r 


(/3 + £ 2 ) 5 (/3 + e) 6 (^ + ^) 7 

1680^ - 280j , 672 0jf - 1680f _ 15120jP 

+ - - T~n H 


(/3 + £ 2 ) 8 


v+ey 


(/3 + £ 2 ) 10 


(B.43) 


Asymptotic formula for (3 — 1 

This case corresponds to real and widely separated saddle points (<r < -1) in the 
integrand of Equation (B.38), with the endpoint 7 arbitrarily close to the saddle 
point Z\ = (-o-) 1/2 , as shown in Figure B.3. For an asymptotic evaluation of (B.38) 
that holds uniformly as the endpoint 7 approaches the saddle point z M we make the 
following transformation [9]: 

q(z) = <rz + z s /3 = q(z ,) + a 2 = + * 2 = t(j) , (B.44) 
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Z — pic 


Figure B.3: Contour of integration of the ii 

7 ~ (-*y /2 . 
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with arg(s) restricted so that J 0 (cr, 7; fi) converges as & — ► 00. Hence, employ- 
ing (B.44) in (B.38) we have 

00 

U<r,r,0) = e-'*"'-'’"/ it, (B.45) 

c 

with the upper limit taken in the sector 0 < arg(a) < 7r/2 of the complex 3 -plane. 
The quantities £ and f(s) are given by 


c = ± 


" t + 773 + 




7 < ( — °’) 1/2 > 


(B.46) 


and 


f /.x _ _ llf) 2^ 


ds q'(z) <r -f z 2 
Equation (B.45) can be written as follows [25]: 


(B.47) 


n) = j /(0) /e-’ * + ^ / [ 


M - /( 0 ) 




(B.48) 


and using integration by parts in the second integral we get 

[/to - m 


U’-r-V) = - Yju [ 


o ;nc 2 


(B.49) 


where 

9(s) = »/■(») - m + m 


(B.50) 


In a similar way, the second integral in Equation (B.49) can be written as follows: 


where 


00 

J g(s)e jn,:i ds 
< 

= *(0)/e'"- d.- 2ja 

'9(C)-9(oy 

c \ 

e *c 7 



~ 2 ja] h ^ e °' d *’ 



(B.51) 


*9'(*) - g(*) + 9 ( 0 ) 


(B.52) 
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Thus, using (B.51), Equation (B.49) becomes 

[/(0) - ^j(O)] / e'"-’ d. 

,v„f 1 [ /(C) -/(0) 1 

\ m [ ( J 




my [ c 


g(C)-g(o) 




(2 jSl) f 

The same procedure is repeated once more for the last integral in (B.53), that is 

m k{0) }J em ' da 

mo - <=(0) 1 j 


(B.53) 


oo 

J h(s)e ]n,i ds ~ 


+ 


M o) 

l 






+ e 


+ o(n-). 


J 1 

1 

'o 

j 

\ 2 jQ 

L c J 


(B.54) 


where 


k{s) = 


sh’(s) — h(s) + fc(0) 


(B.55) 


and finally combining (B.53) and (B.54), the asymptotic expansion of / 0 (o-, 7; Q) when 
(T <C — 1 is given by 




(2ja) 


■Mo) 


OO 

J e jn ' ! ds + e in( 


»-y+-T 


J 1 

/“N 

O 

1 

\ 2jQ 

L c J 


+ 


+ 


MO -g(°)l 


fh(a-h(o) 


my 1 c 

*( 0 - 


my 


-fc(Q) jj 


+ o(n- s ). 


(B.56) 


(2jfi) 4 [ 

In order to express the functions /, <7, h, and k in (B.56) in terms of a and 7, we 
need to derive an expression for f(s) and its derivatives when s « 0. Since f(s) is 
regular near s = 0, it can be expanded in a Taylor series around s = 0 by means of 
Lagrange’s theorem. This is done using a procedure introduced by Erdelyi [26] that 
yields 


f(s) = 1 T(3n/2 1)( l) n 1 

J( ' (-*) 1 * h{n~ l)!r(n/2) [H-ay^r 1 

*L f (^ = 1 V r(3n/2 - !)(-!)— 1 

ds^ S) (— <r) l/4 (n-k- l)!T(n/2) [3(-<r)>' 4 ] 


n-k-l 


(B.57) 

(B.58) 
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Hence, using (B.47), (B.50), (B.52), (B.55), (B.57), and (B.58) we have 
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and 


(B.59) 

(B.60) 

(B.61) 

(B.62) 

(B.63) 

(B.64) 

(B.65) 

(B.66) 

(B.67) 

(B.68) 

(B.69) 

(B.70) 


Also, the Fresnel integral in (B.56) may be expressed in terms of the UTD transition 
function F(x ) [4], that is 


00 p 

y e"”’ is = fi- 1 ' 1 v/re'"' -u(-i)) - 


= n- ,,, A(i )) 1 


(B-71) 


where T) = D‘ /2 £. The general properties of the UTD transition function and de- 
tails on its computation are provided in Appendix D. Now, using (B.59)-(B.71), the 
asymptotic expansion of 7 0 (o-,7;D) when a <C — 1 is given by 

e -j$n(-<o s /2 

*»(*, 75 Q) - ^ (1^/4 A (0 + eJn< ^ / s ) ^ K 7^5 «) + 0( D - S ) , ( B .72) 
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where 
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with a, = 0.20833333, s 2 = 0.33420139, and s, = 1.02581260. 

The asymptotic formulae for g, (/?,£) and its derivative when < -1 a. 
obtained using (B.39)) and (B.72)-(B.74), that is 
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(B.73) 

(B.74) 
•e then 
(B.75) 

(B.76) 

(B.77) 

(B.78) 
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(B.79) 


(B.80) 

(B.81) 


and 
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(B.82) 


It can be easily shown using the large argument form of the UTD transition 
function (see Appendix D) that when £ > (-/?) 1/2 or -q » 1, Equation (B.75) appro- 
priately reduces to the first six terms in (B.42). Also, (B.79) and (B.80) remain finite 
near the caustic when £ — * (— P) l/ ^ and rj — > 0, however, they become numerically 
unstable and an alternative formulation should be used. Applying L’Hospital’s rule 
in (B.56), and then using (B.39) and (B.57)-(B.70) we have 
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(B.84) 


where 

/o = 1/3, /, = 1.25, / 2 = 4/3, 

g 0 = 0.29629630, < 7 , = 3.00781250, g 2 = 5.83333333, 

h 0 = 0.69135802, h, = 4.74804688, h 2 = 13.3333333, 

Jb 0 = 2.63374486, k , = 6.48405151, fc, = 23.8333333. 

Also, A (77) and its derivative near the caustic are given by 


A(t} « 0) ~ 
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;-(-0) l ' J l l/2 ’ 
-0) 1 / 7 . 


(B.85) 

(B. 86 ) 


Asymptotic formula for P » 1 

This case corresponds to imaginary and widely separated saddle points (<r 1) in 

the integrand of Equation (B.38), with the endpoint 7 arbitrarily close to the saddle 
point Z\ = ja ip , as shown in Figure B.4. For an asymptotic evaluation of (B.38) that 
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z-plane 


Figure B.4: Contour of integration of the incomplete Airy integral when a 1 and 
7 ~ ja 1 ! 2 . 

holds uniformly as the endpoint 7 approaches the saddle point z , , we make the same 
transformation as in (B.44), that is 

q(z) = <rz + z s / 3 = + s 2 = + * 2 = ^(*)» (B.87) 

with arg(s) restricted so that 7 0 (<r,7;fl) converges as j -» 00. Hence, employ- 
ing (B.87) in (B.38) we have 

00 

I 0 {<r, 7 ; ft) = J f(s)e itu ' ds , (B.88) 

c 

with the upper limit taken in the sector 0 < arg(s) < 7r/2 of the complex a-plane. 
The quantity ( is given by 

( = dt[<T7 + y/3-Jp''] ,, ’i *(7);0, (B.89) 
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and f(s) is the same as in (B.47). Then using the same procedure as in the case of 
<r < -1, / 0 (<r, 7; ft) when cr > 1 is given by 
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+ 

eiiX-T+T*/.)^ 7, x; fi) + 0(n-> ) , 

(B.90) 

A(x) = 

[nx) e -;,*l 

2 jx \ 

* 

; x = n ,n (, 

(B.91) 


F(x) is the Fresnel transition function of complex argument (see Appendix E), and 
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The asymptotic formulae for g,(/3,£) and its derivative when f3 < -1 are then 
obtained using (B.39) and (B.90)— (B.93), that is 
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where 
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Figure B.5: Four different sets of formulae are used for the computation of the in- 
complete Airy function, one for each region in the Figure, and a fifth set that is used 
in the immediate vicinity of the caustic. 

Using the large argument form of the Fresnel transition function given in Appendix 
E, it can be easily shown that when /? + £’» 1 Equation (B.94) appropriately reduces 
to the first six terms in (B.42). 

B.3 Computational aspects and error assessment 

For an efficient and accurate computation of the incomplete Airy function g,(/3, £) and 
its derivative, the argument space is divided into four regions as shown in Figure B.5. 
Four different sets of formulae are used, one for each region in Figure B.5, and a fifth 
set of formulae that is used in the immediate vicinity of the caustic (/? 4- ^ 0. In 

region I, Equations (B.42) and (B.43) are used, in region II, Equations (B.75) - (B.82) 
are used, in region III Equations (B.94 )-(B. 101) are used, and in region IV, the series 
solution given in (B.19) and (B.20) is used. In the immediate vicinity of the caustic 
and specifically when tj < 0.1 in (B.82), Equations (B.75)— (B. 78) and (B.82)-(B.86) 
are used. 
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An empirical expression for the number of terms used in the series solution is 
given by: 

N(p,() = 8|0| + 4, when £ < 2.0 , (B.102) 

= 8|0| + 4 + 3|0|(£-2.O), when £ > 2.0 , (B.103) 

and results in a truncation error of less than 10"®. 

Figure B.6 shows the percent amplitude error of the first asymptotic result (0 + 
£ 2 1) relative to the series solution along the boundary between regions I and IV. 

The results are plotted vs. the parameter 0, with £ = (12 — 20) 1/2 . The asymptotic 
result shows excellent agreement with the series solution, exhibiting a maximum error 
of 0.12%. Figure B.7 shows the percent amplitude error of the second asymptotic 
result (0 <C — 1) relative to the series solution along the boundary between regions II 
and IV. The results for this case are plotted vs. the parameter £, with 0 = —4. Again 
the asymptotic result shows excellent agreement with the series solution, exhibiting 
a maximum error of only 0.075%. 

Figures B.8 and B.9 show plots of the incomplete Airy function g!(0,£) in the 
first and second quadrants of the 0£-plane, respectively. Figures B.10 and B.ll show 
plots of the derivative of the incomplete Airy function with respect to 0 in the first 
and second quadrants of the 0£-plane, respectively. 

A FORTRAN code for the computation of the incomplete Airy functions based 
on the formulae derived in this appendix is available from the author. 
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" 0.1 


g>( 0' () 


Figure B.6: Percent amplitude error of the first asymptotic result (/? + £ 2 ) for 
the incomplete Airy function (solid line) and its derivative (broken line) along the 
boundary between regions I and IV. Results are plotted vs. the parameter (3 with 
£ = (12 - 2/3) 1/2 . 











Figure B.8: Plots of the real and imaginary parts of the incomplete Airy function, 
8 i(A£)j i* 1 the quadrant of the /3^-plane. 
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ary parts of the derivative with respect to {3 
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Figure B.10: Plots of the real and imaginary parts of the incomplete Airy function, 
gXPiOi * n the second quadrant of the /3£-plane. 
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Figure B.ll: Plots of the real and imaginary parts of the derivative with respect to 
(3 of the incomplete Airy function, g,(/9,£)> in the second quadrant of the /?£-plane. 
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Appendix C 

The Generalized Incomplete Airy 
Integral 


The generalized incomplete Airy integral serves as a canonical function for the uniform 
asymptotic approximation of a class of integrals characterized by multiple symmet- 
rically situated stationary phase points that arbitrarily close to one another or to an 
integration endpoint [22]. It is defined as follows: 


-K«(7>C>*') = J e silt * t '' M dz ; i/ = 3,4,5---, (C.l) 

c 

where 7 and £ are considered real, and the path of integration is any path that goes 
to infinity in the sector 0 < arg(z) < w/i/. When v = 3, K u ( 7,C>3) in (C.l) is the 
definition for the incomplete Airy integral [9]. When v = 4, K u ( 7,C,4) is a special 
case of the incomplete Pearcy integral with the quadratic phase coefficient being equal 
to zero. 

The complex plane topology of the generalized incomplete Airy integral when 
£ > 0 is depicted in Figure C.l. Depending on the endpoint location relative to the 
saddle point may be written as follows: 

1. 7 + C~‘ < 0, C > 0= 


K.(l,M = j' 

2 . 7 + C " 1 > 0, C > 0: 


00 exp(— j 3 w/ 2 t/) 
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J 00 exp( — }&k/ 2 v) 
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/ ooexp(jir/2v) v 

e ,h... M dz 


(0.3) 
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Figure C.l: Complex plane topology of the generalized incomplete Airy integral for 

C > o 
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Figure C.2: Complex plane topology of the generalized incomplete Airy integral for 
£ < 0 and v odd. 


The complex plane topology of the generalized incomplete Airy integral when 
£ < 0 and v odd or even is depicted in Figures C.2 and C.3, respectively. Depending 
on the endpoint location relative to the saddle point K u ( 7,C>* / ) ma y be written as 
follows: 

3. 7 + £ t/ ~ 1 < 0, £ < 0, v is odd: 


^-(7. Ci*') = £ 


oo exp[— j(ir — Sw/2i/] 


foo exp(jir/2i/) 
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" oo exp[j(w-*72»/)] 
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5. 7 + C ‘ < 0) C < 0; y_ > s even: 


poo exp{— j(w — w/2»/] 
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e }(i*+' v Mdz . (C.7) 


6. 7 4- C~ l > 0, C < 0> v is even: 


K,(l,(,‘’) = l 


oo e*p[j(fr— 3ir/2v] 


+ f 

J <x 


oo exp{jn/2v) 
oo eKp[j(w—3ir/2i/)] 


C.l The incomplete Airy integral 


A special case of particular interest for K u ( 7 ,Cj i/ = 3) in (C.l) is the incomplete 
Airy integral which serves as a canonical function for the uniform asymptotic ap- 
proximation of a class of integrals characterized by two stationary phase points that 
arbitrarily close to one another or to an integration endpoint [27]. It is defined as 
follows: 

oo 

Ai(7,C) = J l3) dz , (C-8) 

( 

where 7 and ( are considered real, and the path of integration is any path that goes 
to infinity in the sector 0 < arg(z) < 7r/3. 

The complex plane topology of the incomplete Airy integral when the saddle points 
Zia = ±(-7) 1/2 are real (7 < 0) is depicted in Figure C.4. There exist three cases 
depending on the location of the endpoint £ relative to the saddle points z li2 , and 
Ai(7,C) may be written as follows: 

1. C < -(-7) 1/2 : 

ffifr.C) = •*-*•»*+ / 

if * ^23 *'^'31 

= g 2 (7>0 -iv / ir[W 1 (7) - W,(7)], 7 < 0, (C.9) 


2. -(-7) t/a < C < (-7) 171 : 

Al(7,0 = «*■+'»* + f e^> 3 Uz 

J<i ■'J'Sl 

s g.(7,0 + iv5w,( 7 ), 7 < 0 (C.10) 

3. ( > (~7) ,ff : 

/ oo expfjVi ) gR 

e*"*’ /,) dz = g l (7,0 > 7<0. (C.ll) 

The complex plane topology of the incomplete Airy integral when the saddle points 
Z] 3 = ±(_ 7 )‘/ 2 are complex (7 > 0) is depicted in Figure C.5. There exist two cases 


176 




Figure C.4: Complex plane topology of the incomplete Airy integral for 7 < 0. 
depending on the sign of the endpoint C, and Ai( 7 ,C) may be written as follows: 

4. C < 0: 


eo ”« j ' h) e XT+' 3 Mdz+ f e*** m d* 

Jl-j 


Ai <^> = l 

= 8i(7»() + 2xAi(7), 7>0 


( 0 - 12 ) 


5. ( > 0: 


AI(7,0 -jT" W,, ^»*= fc (7.0, 7>0- 


(0-13) 
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Appendix D 

The UTD Transition Function 


The UTD transition function utd is defined as follows: 

oo 

F(x) = 2 jyfi e jz J e~ ir * dr , (D.l) 

where x is cosidered real, with yfx = \y/x\ if x > 0, and y/x = j\y/x\ if x < 0. Also, 
the path of integration in the Fresnel integral is any path that goes to infinity in the 
fourth quadrant. |sc | ) may be obtained using the following relation: 


f(-|*l) = f(H)- (D.2) 

When x < 6.0, F(x) is computed using its power series form, that is 

F(x) ~ v'ttx e j(l+,r/4) sign(x) - 2 jxe ix £ . (D.3) 

„=o *-) n ' 

The number of terms in the series, A(x), is given by the empirical formula 


N{x) = lOy/x , 


(DA) 


and results in a truncation error of less than 10 - \ The large argument form of the 
Fresnel transition function (x > 6.0) is given by: 


where 


f(») ~ £ 


(- 1 )" 2 " ( 0 . 
(2jn)“ 


(D.5) 



and 2 m = (2a + l)(2a + 3) • • • (2a + 2m - 1) . 


(D.6) 
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Appendix E 

The Fresnel Transition Function 
of Complex Argument 


In this appendix, some properties and formulae for the Fresnel transition function 
with complex argument are presented. The Fresnel transition function is defined as 
follows: „ 

OO 

F(() = 2jf£e J{2 / e _> ' 2 ds (E.l) 

€ 

where the path of integration in the Fresnel integral is any path that goes to infinity 
in the fourth quadrant. When ( is real, F(£) = F(C), where F(x) is the UTD 
transition function [4]. The complex plane topology of the Fresnel integral is shown 
in Figure E.l. Also, F(-£) may be obtained using the following relation: 

F(-0 = F{Q - . (E.2) 


The power series form of F(() used in small argument approximations is given by 



- 2}?'“’ E 

n= 0 


tier . 

(2 n + l)n! ’ 


\t\ <!• 


(E.3) 


The asymptotic expansion of F(() used in large argument approximations is given 
by: 

m ~ a S «(-’>) + 1 + E(-ir 1 : 3 ~ (2 ~ ; lfl»i. ( E -“) 

rj = ( E - 5 ) 
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Figure E.l: Complex plane topology of the Fresnel integral with complex argument 
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Appendix F 

Geometric Interpretation of the 
Phase Derivatives in the 2-D PO 
Radiation Integral 


In this appendix, the phase derivatives in the PO integral formulation for the near- 
zone EM scattering from a smooth two dimensional curved boundary are expressed 
in terms of the local geometric parameters of the boundary C as shown in Figure F.l. 
The phase function <f>(l) in the 2-D PO radiation integral is given by 

m = «*(<> +«(<), (F.i) 

where 

[»'(/)]’ = { [x'CA.) + • af/,,)}’ + { , (f.2) 

K(l) = { [«(/„) - ^<)] - + { [^(<„) - , (F.3) 

t 0 is an arbitrary reference point on the boundary C, l = t 0 + St is the arclength, and 
f(£) is the natural representation of C. 

1. First derivative, 4>{l) : 

From Equation (F.l) the first derivative of the phase function is given by 

4>(t) = s i (t) + R(l). (F.4) 

Taking the first derivative of both sides of Equation (F.2) we have 

+ {[?•(/,) + >=(<)] •<&)} FW-JW,)- (F.5) 
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Using the fact that f (£) = t(l) from the Frenet Equations [28] we have that 


{[?<(<.)+?(<)]•((<„)} 
* (<) - ,‘(l) 

(F-6) 


and letting 8 -» 0, that is l 0 -» t we have that 



<«*> 

II 

(F.7) 


and in a similar fashion 


— 

«(<) = -£(/) •«(/). 

(F.8) 


Thus, using Equations (F.7) and (F.8), the first derivative of the phase function is 


given by 

W) =[*'(<) -£(<)]•*(<)• 

(F.9) 

— 

2. Second derivative, : 

From Equation (F.l) the second derivative of the phase function is given by 



< i(t) = 3\l) + R(l) . 

(F.10) 


Taking the derivative of both sides of Equation (F.5) we get 



[i'Wf + aWW = [^)'^o)] 2 + [^)-^o)] 2 


- 

+ { [s'(Q + f (/)] n(/o)} kl) ■ n(Q 

+ {[i' , (4 + r(/)]n(/„)}r(/).r(£ 0 ), 

(P-11) 

— 

and using the fact that ?(/) = Z g {l) = -u(1)k{1) from the Frenet Equations [28], 

— 

where <c e (/) denotes the normal surface curvature, and again letting 8 -+ 0, 
l 0 — > l along with Equation (F.7) we have that 

that is 

— 

r , i [*'(<) •*(<)]’ 

* V) = -**{*) [* (0 • n ( / )] + ,.•(/) a i(t ) 

(F.12) 

— 

In a similar fashion, R{£) is given by 


— 

r ^ i 1 1^(0 * *(0] 

R(l) = *.(<) [*(<) ■ *(/)] + m m . 

(P-13) 

— 
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and thus using Equations (F.12) and (F.13), the second derivative of the phase func- 
tion is given by 


[£(<)-■? '«)]■»(<) [r(<). «(<)]’ [*«•»«)]’ 

1 ’ p,w m + m 

where p g (£) = 1/k 3 (£) is the surface radius of curvature at £. 

3. Third derivative, <f> (£) : 

The third derivative of the phase function is given by 


(F.14) 




(F.15) 


Taking the derivatives of both sides of Equation (F.ll) we get 


•■(/) + «*'(*) s’(*) 


-3 [((<)• »(4)] [«(<)• 8(4 )]«,(<) 

{ [* '(4) + ?(<)] • 8(4} [«(<) ■ 6(4)] n,(t) 

{ [.‘(4) + f (<)] • 8(4)} [»(/) ■ 8(4)] *.(<) 

3 [<(<)• <(4)] [S(/) • t(4)] «„(<) 

{ [.-(4) + ?(/)] • 1(4)} [«(<) ■ <(4)] K,(t) 

{ [i"(4) + f (<)] ■ i( 4)} [»(/) • 1(4)] k.(l) , (F.16) 


and using the fact that k(l) = t(l)K a {l) from the Frenet Equations [28], and again 
letting 8 ► 0, that is l 0 * l along with Equations (F.7) and (F.12) we have that 


>\i) [r« • «(<)] k,(t) - [»'(<) • i(i)] <(<) 

3 P‘(<) • 8(/)] [.'(/) • (’(/)]*,(<) 

m 

3 [?'(<) .n(<)] ! [}•(<)•((«)] 

[»•(<)]’ 

In a similar fashion, R (/) is given by 

R(i) = [£(<) •«(/)] *,(<) + [*(<)• <(/)]»’,(<) 

3 [*(<) ■ 8(<)] [«(/) • ((<)] «,(/) 

m 

3 [«(<)■ «(<)]’[«(/) •(“(/)! 

R?{1) 


(P.17) 


(F.18) 
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and thus using Equations (F.17) and (F.18), the third derivative of the phase function 
is given by 

[j •'(*) - R{1) 1 • *(*) 

<f> (^) = [-R(^) - * *(^)] ■ pi ^ 


+ 


+ 


3 p'(<)-8(<)] [?•(<) •<(<)] 

*'(<) 

3 [«(<) •»(<)] [flW-fW] 

m 


1 V(l) ■ n(<)l 


P,(l) «'(*) 

1 + R(l)-n(l) 


[pm m 


(F.19) 


Finally, using the following relations: 

?'(/) • n(/) = -sin0' 

j'(^)-r(£) = -cos^' 
R(£) • n(l ) = sin 9 
R{l)‘t(t) = cos 6 


(F.20) 

(F.21) 

(F.22) 

(F.23) 


Equations (F.9), (F.14), and (F.19) can be expressed in terms of the angles 6 and 6' 
defined in Figure F.l, that is 


(f>{£) = — (cos 0 + cos O ') , 

M/ N sin 0 + sin 0' sin 2 O' sin 2 0 

m = W) + ^W + " W)' 

v , cos 6 + cos 9' 
<f> (^) = (sin 9 + sin 0 ')k # (^) + 


(F.24) 

(F.25) 


+ 


+ 


3 sin 9' cos 9' 

W) 

3 sin 9 cos 9 


1 


+ 


m 

sin 9‘ 1 


m 


1 sin 9 

+ 


I pm 


(F.26) 


When l is a reflection point as shown in Figure F.2, Equations (F.24)-(F.26) 
reduce to the more familiar expressions 


fa) 

fa) 

P r 


(I 1 - ? r ) • ? r = 0 , ( law of reflection) (F.27) 

cos 2 9' ( — + — ] oc (reflected ray spread factor) -2 , (F.28) 

V* r P T ) 

— -| = (reflected ray caustic distance) -1 , (F.29) 

s' pg(Q r) cos 9 * 
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SOURCE 


RECEIVER 



-9 i 
90 ° + 6 
90 ° + 9 
e x i 


+ 

Figure F.2: Reflection from a 2— D curved boundary 


8 


+ (lr) 


(X 


2 COS 9'K g {Q T ) 


3 cos 2 8' sin 8' 

a 1 


1 V 

Pg{Qr) cos 8’ a * 


3 cos 2 8 r sin 6 r 

a r 

(reflected ray 


1 r 

p g (Q r ) cos 8 t + s r J 

spread factor near a caustic) 


(F.30) 
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Appendix G 

Geometric Interpretation of the 
Phase Derivatives in the 3-D PO 
Radiation Integral 


In this appendix, the phase derivatives in the PO integral formulation for the near- 
zone EM scattering from a smooth three dimensional curved surface are expressed in 
terms of the local geometric parameters of the surface S as shown in Figure G.l. The 
phase function <f>(£ , r) in the 3-D PO radiation integral is given by 

<t>(t,T) = S i (l,T) + R(l,T), (G.l) 

where 

[s'(/,r)] = {[sj + ri^,r)] -4} 2 + +f(/,r)] -r 0 } 2 

+ {[*0 +^»r)] -n 0 } 2 , (G.2) 

R 2 {£,t) = {[^o -r(£,r)] -4} + {[flo - f^,r)] • f 0 } 2 

+ {[^o -r(^r)] -n 0 } 2 , (G.3) 

4, To, and n 0 are the principal unit tangents and unit normal vector at the reference 
point (l 0 , t 0 ) forming a local surface orthogonal coordinate system, l = l 0 + 61 and 
t = r 0 + er are the arclengths from (^ 0 , r 0 ) to (l, r ) along the principal coordinate 
curves of 5 which are always orthogonal, and the surface vector is therefore 

a natural representation of S for which the first and second fundamental forms are 
given by the following expressions [29]: 

I(d£,dr ) = de + dr\ (G.4) 
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Figure G.l: Geometric parameter definitions for the phase derivatives in the 3-D PO 
radiation integral 
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II(d£,dr ) = -«,(£, r)d^ - /e 2 (/,r)dr 2 , (G.5) 

where /ci i2 are the principal curvatures along the principal surface directions 1 and f, 
respectively. Equations (G.4) and (G.5) imply the following relations: 

r t (i, T )-r t (l,T) = 1, (G.6) 

n(l,T) ‘f T (i,r) = 0, (G.7) 

f r (£,T)-fr(£,T) = 1, (G.8) 

ru{l,r) = — «(£, t ) k , (1,t) , (G.9) 

hr(i,T) = 0, (G.10) 

r(^» t) = t ) ■ (G.ll) 

1. First derivatives, <fo(l,r) and <ft r (l,r): 

From Equation (G.l) the first derivatives of the phase function are given by 

M^ t ) = s\(l,r) + Ri{1,t) , (G.12) 

4> T {l,r) = s' t {1,t) + R t (1,t). (G.13) 


Taking the first derivative with respect to l of both sides of Equation (G.2) we have 


s'(^,r)i)(^,r) = {[s‘ + ^> r )] *£>}rir(/,r) -l 0 

+ { [*'o + t)] • Tb} r t {l, t ) • % 

+ {[^o + ^> r )] *»o}r/(/,T) -n 0 . (G.14) 

Using the fact that f/(i,r) = 1 from (G.6) and (G.7) and letting 6,e -» 0, that is 
(4 >,To) — » (^,r), Equation (G.14) becomes 

= ?*(/, r)-/, (G.15) 

and in a similar fashion J7/(^, r) is given by 

Ri{l, r) = -R(l, t) • 1. (G.16) 

Furthermore, using the fact that = f from (G.7) and (G.8) we have 

K(t> T ) = *'(*,r) -f and Rt{1,t) = R(1,t) ■ f , (G.17) 
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thus, the first derivatives of the phase function , r ) and T ) are given by 


M^t) = [s‘(*,r)-£(*,r)] -I, (G.18) 

Ml,r) = [r(/,r)-^(/,r)] -f. (G.19) 

In terms of the angles d' t>r and 6 i<r defined in Figure G.l, the first derivatives of 
the phase function may be written as follows: 

= -(cos#* + co60/) , (G.20) 

<f> T {t,r) = — (cos6' r + cos0 r ) . (G.21) 

2. Second derivatives, ^ T (^,r), and ^> TT (l,r): 

From Equation (G.l) the second derivatives of the phase function are given by 

fc(*,r) = W,r) + R u (l,r), (G.22) 

= s' lr (t,T) + R tr (l,T) , (G.23) 

<M*,T) = sUl,T) + Rrr(l,T). (G.24) 


Taking the derivative with respect to l of both sides of Equation (G.14) we get 

[ij(£, r )] 2 + = [r/(/,r) • £ 0 ] + [*/(*> r ) * ? o] 

+ r) • n 0 ] 2 + { [s o + r(l, r)] • £>} r) • /„] 

+ { [»i + r(t, t)] • Tb} [r«(^,r) • f 0 ] 

+ {[r‘+r^,r)].n 0 }[f^,r)-no] . (G.25) 

Letting £,e -» 0, that is (1 0 ,t 0 ) -» (*,t) and using (G.9), Equation (G.25) becomes 

[i‘(£,r)] 2 + s\1,t)? u {1,t) = 1 - k,(*,t) [s‘(*,r)-n] . (G.26) 

Directly from (G.26) one can also infer the following expression: 

[i‘ T (*, r)] 2 + s\l, t) = 1 - k 2 (1, r) [i“(*, r) • »] . (G.27) 

Taking the derivative with respect to r of both sides of Equation (G.14), then letting 
8,e — ► 0 and using Equation (G.10) we also get 

ij(*, T)a\.(i, r) + s'(l , r)Sj r (/, t) = 0 . (G.28) 
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Then, using Equations (G.15), (G.17), and (G.26)-(G.28) we obtain the following 
expressions: 


. r i 1 — *(A r ) ‘ 

5 «(At) = -*i(*.r) [a'(i,r)-n] + — ±- 

8 ^4, r ) 

. . [f(i,T)-e]{s%r)-f} 

lT ’ 1 s‘(t,r) 

KM,r) = -«,(<, r) [?'(<, t) • n] + 1 ~ l? ’ (< ’ r) ' f] ‘ 


(G.29) 

(G.30) 

(G.31) 


Similarly, the quantities Ru, i?/ r , and R TT are given by 


~l2 


RMr) = g] + 1 

R{£,t) 

Rtt(£,t) = k 7 (£, t) [^(£,r) • raj + — " " 


(G.32) 

(G.33) 

(G.34) 


Finally, using (G.29)— (G.34) the second derivatives of the phase function are given 

by 


~l2 


- v , 1 - [y'(/, r)./l 

4 > u ( 1 , t ) = «i(At) [^(At) - s’^r)] -n + — J 

S \l,T) 


+ 




R(l, r) 

, . [J i (*,r)-/]P-(<,T).?1 [R(l,r)-e][R(e,r)-r] 

^ ' T R(l,r) 


(G.35) 
, (G.36) 


+ 


1 - [fl(l,r) -f]' 

R(£,t) 


‘(At) 


(G.37) 


In terms of the angles 0' n , 6 n , d^ r , and d/ T defined in Figure G.l, the second 
derivatives of the phase function may be written as follows: 

MAt) = (cos 8' n + cosfl ra )/c,(l,r) + ^ — , (G.38) 

*'(A t) ^R(At) 

7 .. . COS COS 0' COS cos 0 r 

m ' t) = ^#r + ^7T- < G39 > 
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(G.40) 


= (cos#; 


3. Third derivatives, 4> m (l,r), <^ r (l,r), <ft rT/ (l,r), and # rrT (1 ,t): 

From Equation (G.l) the third derivatives of the phase function are given by 


&UI (^> r ) — 

s iu(l> T )+ Rut (*> r )> 

(G.41) 

<f>Uri^ T ) = 

S Ut{^i T )^~ Rut (^i r )? 

(G.42) 

l>rr l(l, T ) = 

3 rr/(^ T )~^~ Rrrl r ) ? 

(G.43) 

4 >TTr(t> T ) = 

3 rrr(^i 7 ") - ^ - Rttt (A r )' 

(G.44) 


Taking the derivative with respect to t of both sides of Equation (G.25) we get 

3i)(£,r)s ‘ u (8, r) + a‘(/,r) a iu(t> T ) = 3 [^(^> r ) *^°] [^w(^> T ) *^o] 

+ 3 [f*(*, r) • f 0 ] t) ■ f 0 ] + 3 [r £ (/, r) • n 0 ] r) • n 0 ] 

+ {ft + f&T)].^} 

+ { o + ^> r )] ‘ r 0 } 

+ {[i’o + ^r)] -n 0 } 


riu {f-,r)- l 0 

rui (A r ) ’ T 0 
rui {l, t) • % 


(G.45) 


Next, from Equations (G.6)-(G.ll) one can show that rui (^> r ) is given by 

f£a(*,r) = -*;(/, r)f- [/• ■ Vis^r)] », (G.46) 

and letting 8, e — ► 0 Equation (G.45) becomes 




-K\(t,T) [r(<,r)-<] 
?.V«,(/,T)[r(/,r)-S] 


(G.47) 


Directly from (G.47) one can also infer the following expression: 

3i^(/, r )S j. T (/, t) + a'(l, r)~a ' TTT (l, r) = -/c^,r) [fl'(/,r) • f ] 

- f-Vm(^r)[r(/,r).n] . 


(G.48) 
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If 


Next, taking the derivative of Equation (G.28) with respect to l and r we obtain the 
following expressions: 


23 \{1,t)s' 1t (1,t) + s' t {1,t)s' u {1,t) + «« r (/, r ) = 0 , (G.49) 

2 it(Ar)iUAr) + i^,r)it T (^r) + - < (£ f r);^.0 = 0- (G.50) 

Then, using (G.15), (G.17), (G.29)-(G.31), and (G.47)-(G.50) the quantities s 

3 ttIi 411(1 3 rrr are g^en by 

S 'tu{t-> T ) = ~ K i(^ T ) [ a '(^» r ) ‘t] -t- Vk,(/,t) [?'(/, r) -n] 

s‘(^, r) • I 


+ 


a «r(*.r) — 


s‘(£,r) 

s«(^,r) 
s'^r) ■ l 


S TTi( i , r ) = ~ | M*, r) [g*(l, r) • n] - 

= -«;(/, r)[r'(/,r).f]-f.V^(/,r)[r'(/,T).n] 
s'( 1 ,t) -t 


+ -;.(V) i 3 * 2 ^’ r ) • "] - 


s*(^,r) 


1-3 [j^r) 

• <T i 

s l (l,r) 

J 

1 - • 

f] 2 ) 

3 '(t, T ) 

I 

r)[a’(£,r)-n] 


1 - [«'(*, r) • 

r-]M 


s*(£,r) 


> , (G.52) 
, (G.53) 

, (G.54) 


Similarly, the quantities Rut, Rut, Rrrt , and R TTT are given by 


+ 


RUt ( l,T ) = 
Rrrt t ) — 


R(1,t ) • f 

W 

R{t, 


3«i(/, t ) [.R(£,t) -n] + 
Ki(£,t) |#(/,r) -n] + 




1 - [«(/, T-) - /] 

R(t,r) 


R(1,t) 


B&f { T) T) " ^ + ^ 

Rrrr(l,T) = *’(<, T> [S(<, r) ■ f ] + f • V*,(/, t) [%,t) • n] 


2 ^ 


, (G.55) 
} , (G.56) 
, (G.57) 


+ 
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Finally, using Equations (G.41)-(G.44), and (G.51)-(G.58) the third derivatives of 
the phase function are given by 

itu&T) = [^, r )-? Ur )]. {«?(/, r )/+ p . V «,(/, T )] n } 


S'U,T)-l I . \ ~1 1 ~ [*‘(^ T ) '*] 

+ ,•’(/, r ) | 3 ' Cl (/, r ) t a (^ T )'”] s‘(£,t) 

_ , 1- ffl(/,r).F| 2 ' 

3*i(/,r) [R(t,r)-n] + L ^ r) - J ~ 

t i 

S'd t) • t f r l 1 - 3 [?'(/, r) ■?] 


R(t,r) 


, (G.59) 


+ 




.rg,. ^ l-3[fi(/,r)-f] 2 

+ J2(/,r) |' 5,(/,r )[ Wr )' n J J2(/,r) 

i rTT (/, r ) = [£(/, t) - r (/, r)] ■ {«,(/, t)t + [f • Vk 2 (/, r)] n} 

1 -[?*(*, r)-ff 


1 2 


1 - 3[?‘(*,r) ■ f] 


, (G.61) 




+ 


l 3 *^ [S( '’ T) ' sl + ‘"fc !" 1 


, (G.62) 


In terms of the angles 9' n , 9 n , 9' tr , and 9 t , T defined in Figure G.l, the third deriva- 
tives of the phase function may be written as follows: 

<j>ui(l, T ) = (cos0J + cos0/)kJ(/,t) + (cosflj, + cos0„) [/• Vk,(/,t)] 


cos 9', [„ a , /s \ . sin* 9' t 
+ ~~7~a r 3 cos 9 k 1 (1, t) + ... v 

s‘(i,r) [ 


+ 


cos 9 1 
R(1,t) 
cos 9' 


3 cos 9„k 1 (£,t) + 


sin* 9i 


m,r) 
1—3 cos 1 9' t 


(G.63) 


tUD l/„ . . . A. U V/VO 

^ (< ’ T) = 7{iyj 


+ 


cos 9 t 

AM 


cos0 n /c,(/,r) + 


1 — 3 cos 2 6 1 


R(l,r) 


(G.64) 
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COS 6 1 


8'{l,T) [ 

cos 6 1 
R{1,t) [ 


COS 0' n K 2 (t, t) + 
cos 8 n k 2 (1,t) + 


1—3 cos 2 0' T 


*% T ) 

1 — 3 cos 2 0 T 


4>lrr{ t -> T ) = 

+ 

4 > ttt ( 1 ,t) = (cos^ + cos 0 r )K 2 (£,T) + (cos^ + cos 0 „) [f • V/c 2 (£,t)] 




(G.65) 


+ 

+ 


cos 6' 


sin 2 O' 


LUD . . . om w- 

-^[3 cos «„*,(<• r)+__j 


cos 6 r 
R(t,r) 


3 cos 6 n K 2 (t, r) + 


sin 2 0 r 


(G.66) 


fl(*,r)J ■ 

When (£,r) is a reflection point Q r = (£ r ,T r ), as shown in Figure G.2, one can 
use the derivatives of the phase function in the PO radiation integral to derive some 
familiar expressions involved in the reflected field from a 3— D surface. The first 
derivatives of the phase function given in (G.18) and (G.19) evaluated at Q r represent 
the law of reflection, that is 


f>t(lr,T r ) = (s i -S r )-l= 0 

(, K{£r,T r ) = (3“ - ? r ) • ? r = 0, 

which can be written in a more familiar form as follows: 

(?' — s r ) • n r = 0 
(?' — s r ) x n r = 0 . 


(G.67) 


(G.68) 


The second derivatives of the phase function given in (G.38)— (G.40) evaluated at Q r 
with (G.67) being employed are related to the spread factor of the reflected wavefront 
through the following expression: 

M£ r ,T r )4> rr {i r ,T r )-f>l T {l r ,T r ) = cos 3 O' ^ (G.69) 

oc (reflected ray spread factor) -2 , 


where O' is the usual angle of incidence, and p\ T are the radii of curvature of the 
reflected wavefront due to spherical wave illumination and are given by 



± 


_1_ 1 
s' cos O' 


[«l (Qr) sin 2 0 T + K 2 {Qr) Sm 2 0/] 



4*i(Qr)*2(Qr) 5 


(G.70) 
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n r = OUTWARD UNIT NORMAL AT Q r 
t r = UNIT TANGENT AT Q r 
£ r = UNIT BINORMAL AT Q r 
t r ,T r = PRINCIPAL DIRECTIONS AT Q r 
e-j- = gi x s 1 
e ( r = ef x s r 


INTERSECTION OF THE 

PLANE OF INCIDENCE WITH S 


Figure G.2: Reflection from a 3-D curved surface 
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where the angles 6^ r are defined in Figure G.2 and K\^{Qr) are the principal surface 
curvatures at Q T in the directions of l r and f r , respectively. An alternative expression 
for p r t is as follows: 


1 = 1 f /«g(Qr) 

p T t 3' COS 0 ' 


1± 


N 


1 - 


4«i(Q r )K 2 (Q r )cos 2 0* 

pKg(Qr) 


f 1 I J at 

f = 1 + ~ \ cos 0 > 


(G.71) 


(G.72) 


Kg(Qr) 

where n g (Q r ) is the normal surface curvature along the plane of incidence, Kt(Qr) is 
the normal surface curvature along the transverse plane defined by the unit normal 
vector n r and the unit binormal vector b r , and they can be found using Euler’s 
theorem [29], that is 


K g(Qr) = «i(Qr) cos 2 a + K 2 (Q r )sin 2 a , (G.73) 

K t(Qr) = «i(£? r )sin 2 a + « 2 (C?r)cos 2 a , (G.74) 


where a is the angle between t r and l r as shown in Figure G.2. The third derivatives of 
the phase function given in (G.63)-(G.66) evaluated at Q r with (G.67) being employed 
are related to the spread factor of the reflected wavefront when the observation point 
is near a caustic (p r f = — s r ) through the following expression: 


m(Q,) = - t tr[lr ' Tr) 


^ ^Ur (^r j T r ) 


^rr(4 ? ) 

r ? T r )<h rr 


^rr(frj ^”r) 

oc (reflected ray spread factor near a caustic) _s , 


(G.75) 


where 


$U1 {^ry T r) 

= 

2 cos 0' 

[tr ■ v« x (g r )] 



+ 

cos 9' t 

s' 

3 cos 0*/c 1 (Q r ) + 

sin 2 6' t 

s' 


+ 

cos Of 
8' 

3 COS 8 r K 1 (Q r ) + 

sin 2 0/ 
s r 

^r(^rj T r ) 

= 

— cos 6/ cos 0, ( i + — ^ 
\s' S T ) 

1 . 

^rr(4 ? ) 

= 

2 cos 0* 

K 2 (£? r ) + sin 2 0 r I 

n i 

{S' s r 


(G.76) 
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(G.77) 

(G.78) 




Figure G.3: Diffraction from a curved edge in 3-D surface (side view) 
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cos 6 1 
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1 — 3 cos 2 B' t 
a* 

1 — 3 cos 2 0( 
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1—3 cos 2 8’ r 
a ’ 

1 — 3 cos 2 0 T 

b t 


(G.79) 


(G.80) 


When ( 1,t ) is an edge diffraction point Q a = (/ aj Ta) as shown in Figure G.3, 
one can use the derivatives of the phase function to derive some familiar expressions 
involved in the diffracted field from a curved edge in a 3-D surface. It is assumed here 
for convenience that the edge extends along the surface coordinate curve r. The first 
derivative of the phase function along r given in (G.19) evaluated at Q a represents 
Keller’s law of edge diffraction, that is 
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(G.81) 


ir (4,’i) = (i; - ?a) To =0, 

where f a is the associated unit tangent to the edge at Q a . The second derivative of the 
phase function along r given in (G.40) evaluated at Q a with (G.81) being employed 
is related to the spread factor of the diffracted wavefront, that is 

<f>rr(la,T a ) = am* ^ 

oc (diffracted ray spread factor) -2 , 


(G.82) 


where is the angle of diffraction defined in Figure G.2, and p'i is the distance from 
the caustic at Q a and the second caustic of the diffracted ray and is given by 


J_ _ l_ (»j ~ 

Pi + R° sin 2 i? a 


(G.83) 


where n a is the associated unit normal vector to the edge at Q a directed away from 
the center of curvature, and R a > 0 is the radius of curvature of the edge at Q a . 
The third derivatives of the phase function given in (G.63)— (G.66) evaluated at Q a 
with (G.81) being employed are related to the distance parameter for the EUTD 
transition function, that is 

(p/.T (^(1 , T a ) 


Ld{Qa) — <f>lU (ta, T a) ~ 


3 < Pllr (^aj T “) 


^tt(^O) ) 

a 5 Ta ) (^a? T a) 

^Tr(^a j Ta) 


(G.84) 


where 


&UI (^ai T a) — (cos 0/a + COS 

+ sint? a (sin <p' a + sin<p a ) [£. • /c,(Q a )] 
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cos 0 


la 


*a 

COS 0 


3ki (Q a ) sin i? a sin <p' a + 
3«1 ( Qa ) sin t? a sin <p a + 


sin 2 0'tn 


sin 2 0 ta 


.. tt _ S COS 0'fo cos 0' Ta _ COS 0fax COS 0 r 

*PlT\£a)1~ a J ' 


si 


(G.85) 

(G.86) 


$ Ur (^i r a) — 


cos 


0' 


S a 


«i(Qa)sin i? a sin <p' a + 


1 — 3 cos 2 0, 


la 


si 
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Figure G.4: Diffraction from a curved edge in a 3-D surface (top view) 



(G.87) 


(G.88) 


The angles <p' a and <p a are the angles between the projected directions of incidence 
and observation (on a plane perpendicular to ?„) and l a as shown in Figure G.4. 


202 



Bibliography 


[1] M. Klein and I. W. Kay, Electromagnetic Theory and Geometrical Optics. Inter- 
science Publishers, 1965. 

[2] H. Ikuno and L. B. Felsen, “Real and complex rays for scattering from a target 
with inflection points,” Radio Science, vol. 22, No. 6, pp. 952-958, Nov. 1987. 

[3] H. Ikuno and L. B. Felsen, “Complex ray interpretation of reflection from 
concave-convex surfaces,” IEEE Trans. Antennas Propagat., vol. 36, pp. 1260- 
1271, Sept. 1988. 

[4] R. G. Kouyoumjian and P. H. Pathak, “A uniform geometrical theory of diffrac- 
tion for an edge in a perfectly conducting surface,” Proc. of IEEE, vol. 62, pp. 
1448-1461, Nov. 1974. 

[5] J. B. Keller, “Geometrical theory of diffraction,” J. Opt. Soc. Am., vol. 52, pp. 
116-130, 1962. 

[6] R. F. Harrington, Time Harmonic Electromagnetic Fields. New York: McGraw- 
Hill, p. 127, 1961. 

[7] C. Chester, B. Friedman, and F. Ursell, “An extension of the method of steepest 
descent,” Proc. Cambridge Philosoph. Soc., vol. 53, pp. 599-611, 1957. 

[8] L. B. Felsen and N. Marcuvitz, Radiation and Scattering of Waves. Englewood 
CUffs, NJ: Prentice Hall, 1973. 

[9] L. Levey and L. B. Felsen, “On incomplete Airy functions and their appUcation 
to diffraction problems,” Radio Science, vol. 4, No. 10, pp. 959—969, Oct. 1969. 

[10] M. H. Rahnavard and W. V. T. Rusch, “Surface-curvature-induced microwave 
shadows,” IEEE Trans. Antennas Propagat., vol. 30, pp. 83-88, Jan. 1982. 

[11] N. C. Albertsen, P. BalUng, and N. E. Jensen, “Caustics and caustic corrections 
to the field diffracted by a curved edge,” IEEE Trans. Antennas Propagat., vol. 
35, pp. 297-303, May 1977. 

[12] V. P. Maslov, Perturbation Theory and Asymptotic Method. Paris, France: 
Dounod, 1972 (in French). 


203 



[13] R. W. Ziolkowski and G. A. Deschamps, “Asymptotic evaluation of high fre- 
quency fields near a caustic: An introduction to Maslov’s method,” Radio Sci- 
ence, vol. 17, pp. 1181-1191, 1984. 

[14] P. H. Pathak and M. C. Liang, “On a uniform asymptotic solution valid across 
smooth caustics of rays reflected by smoothly indented boundaries,” IEEE Trans. 
Antennas Propagat., Vol. 38, No. 8, pp. 1192-1203, Aug. 1990. 

[15] R. M. Lewis and J. Boersma, “Uniform asymptotic theory of edge diffraction,” 
J. Math. Phys., Vol. 10, pp. 2291-2305, 1968. 

[16] N. A. Logan and K. S. Yee, Electromagnetic Waves, edited by R. E. Langer, 
University of Wisconsin Press, Madison, 1962 

[17] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions. New 
York: Dover, 1972, p. 446. 

[18] H. Ihuno, “Calculation of far-scattered fields by the method of stationary phase,” 
IEEE Trans. Antennas Propagat., Vol. AP-27, pp. 199-202, Mar. 1979. 

[19] P. Ya. Ufitsev, “The method of fringe waves in the physical theory of diffraction,” 
Sovyetskoye Radio, Moscow, 1962. Now translated and available from the US Air 
Force Foreign Technology Division, Wright Patterson AFB, Dayton, OH, USA. 

[20] G. L. James, Geometrical Theory of Diffraction for Electromagnetic Waves. Lon- 
don, UK: Peter Peregrinus Ltd., 1976, p. 150. 

[21] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions. New 
York: Dover, 1972, p. 260. 

[22] E. D. Constantinides and R. J. Marhefka, “On a class of complete and incomplete 
generalized Airy functions and their application to high-frequency scattering and 
diffraction problems,” paper presented at the 1992 IEEE AP-S/URSI Sympo- 
sium, Chicago, IL. 

[23] G. B. Airy, “On the intensity of light in the neighborhood of a caustic,” Trans. 
Cambridge Phil. Soc., Vol. 6, p. 379, 1838. 

[24] V. A. Fock, Electromagnetic Diffraction and Propagation Problems. Oxford, Eng- 
land: Pergamon, 1965. 

[25] R. M. Lewis, “Asymptotic theory of transients,” Electromagnetic Wave Theory, 
edited by J. Brown. Oxford, England: Pergamon, 1967, pp. 864-869. 

[26] A. Erdelyi, Asymptotic Expansions. New York: Dover, 1956, pp. 41-43. 

[27] N. Bleistein, “Uniform asymptotic expansions of integrals with many nearby 
stationary points and algebraic singularities,” J. Math. Mech., Vol. 17, pp. 533- 
559, 1967. 


204 



[28] M. M. Lipschutz, Differential Geometry, Schaum’s Outline Series in Mathemat- 
ics, McGraw-Hill, 1969. 

[29] D. J. Struick, Lectures on Classical Differential Geometry, Addison- Wesley, 
1950. 


205 




